Stockage de la matrice

Nous devons stocker uniquement les coefficients de la matrice situés dans la demi-bande inférieure, soit le $ (M+1)\times N$-tableau suivant.

\begin{figure*}\centerline{
\psfig{file=5/IMAGES/stock_lap.eps,height=3.5truecm}
}
\end{figure*}

La bibliothèque LAPACK étant écrite en FORTRAN, ce tableau est rangé en mémoire dans un vecteur de longueur $ (M+1)\times N$ par colonne. Si nous remplissons la matrice depuis un programme écrit en langage C, et puisqu'en langage C les tableaux sont rangés par ligne, il nous faudra stocker la matrice dans le $ N\times (M+1)$-tableau transposé du tableau ci-dessus.

EPFL-IACS-ASN