| 1 |
\newpage |
| 2 |
\begin{subappendices} |
| 3 |
\section{Algorithme $LU$ version Crout pour les systèmes tridiagonals} |
| 4 |
\label{LUtri} |
| 5 |
\begin{algorithm}[ht] |
| 6 |
\SetAlgoLined |
| 7 |
%\LinesNumbered |
| 8 |
\KwData{\\$A$ la matrice tridiagonale d'un système $AX=B$ avec les termes diagonaux qui sont déjà pivot} |
| 9 |
\KwResult{Valeur du vecteur $X$} |
| 10 |
\tcc{Décomposition de $A$} |
| 11 |
\For{$i=1$ to $2$} |
| 12 |
{ |
| 13 |
$L(i,1)=a(i,1)$ |
| 14 |
} |
| 15 |
$u(1,2)=\frac{a(1,2)}{a(1,1)}$\\ |
| 16 |
\For{$i=2$ to $n-1$} |
| 17 |
{ |
| 18 |
$l(i,i)=a(i,i)-l(i,i-1)*u(i-1,i)$\\ |
| 19 |
$l(i+1,i)=a(i+1,i)$\\ |
| 20 |
$u(i,i+1)=\frac{a(i,i+1)}{l(i,i)}$\\ |
| 21 |
} |
| 22 |
$l(n,n)=a(n,n)-l(n,n-1)u(n-1,n)$\\ |
| 23 |
\tcc{Descente triangulaire} |
| 24 |
$y(1)=\frac{b(1)}{l(1,1)}$\\ |
| 25 |
\For{$i=2$ to $n-1$} |
| 26 |
{ |
| 27 |
$y(i)=\frac{b(i)-l(i,i-1)y(i-1)}{l(i,i)}$ |
| 28 |
} |
| 29 |
\tcc{Remontée triangulaire} |
| 30 |
$x(n)=y(n)$\\ |
| 31 |
\For{$i=n-1$ to $1$} |
| 32 |
{ |
| 33 |
$x(i)=y(i)-u(i,i+1)x(i+1)$ |
| 34 |
} |
| 35 |
|
| 36 |
|
| 37 |
\caption{Algorithme de décomposition de Crout pour une matrice A tridiagonales} |
| 38 |
\end{algorithm} |
| 39 |
\end{subappendices} |
