1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465 |
- \documentclass{article}
- \usepackage[pdftex,active,tightpage]{preview}
- \setlength\PreviewBorder{2mm}
- \usepackage[utf8]{inputenc} % this is needed for umlauts
- \usepackage[ngerman]{babel} % this is needed for umlauts
- \usepackage[T1]{fontenc} % this is needed for correct output of umlauts in pdf
- \usepackage{amssymb,amsmath,amsfonts} % nice math rendering
- \usepackage{braket} % needed for \Set
- \usepackage{algorithm,algpseudocode}
- \begin{document}
- \begin{preview}
- \begin{algorithm}[H]
- \begin{algorithmic}
- \Function{Gauss}{$A \in \mathbb{R}^{n \times {n+1}}$}
- \For{($i=1$; $\;i \leq n$; $\;i$++)}
- \State // Search for maximum in this column
- \State $maxEl = |A_{i,i}|$
- \State $maxRow = i$
- \For{($k=i+1$; $\;k \leq n$; $\;k$++)}
- \If{$|A_{k,i}| > maxEl$}
- \State $maxEl = A_{k,i}$
- \State $maxRow = k$
- \EndIf
- \EndFor
- \\
- \State // Swap maximum row with current row
- \For{($k=i$; $\;k \leq n$; $\;k$++)}
- \State $tmp = A_{maxRow,k}$
- \State $A_{maxRow,k} = A_{i,k}$
- \State $A_{i,k} = tmp$
- \EndFor
- \\
- \State // Make all rows below this one 0 in current column
- \For{($k=i+1$; $\;k \leq n$; $\;k$++)}
- \State $c = -\frac{A_{k,i}}{A_{i,i}}$
- \For{($j=i$; $\;j \leq n$; $\;j$++)}
- \If{$i == j$}
- \State $A_{k,j} = 0$
- \Else
- \State $A_{k,j} += c \cdot A_{i,j}$
- \EndIf
- \EndFor
- \EndFor
- \EndFor
- \\
- \State // Solve equation for an upper triangular matrix
- \State $x = \Set{0} \in \mathbb{R^n}$
- \For{($i=n$; $\;i \geq 1$; $\;i$-\--)}
- \State $x_i = \frac{A_{i,n+1}}{A_{i,i}}$
- \For{($k=i-1$; $\;k \geq 1$; $\;k$-\--)}
- \State $A_{k,n+1} -= A_{k,i} \cdot x_{i}$
- \EndFor
- \EndFor
- \\
- \State \Return $x$
- \EndFunction
- \end{algorithmic}
- \caption{Gaussian elimination}
- \label{alg:seq1}
- \end{algorithm}
- \end{preview}
- \end{document}
|