Nico a écrit :

Michel Talon a écrit :

       SGBSV - compute the solution to a real system of linear equations
A * X
       = B, where A is a band matrix of order N with KL  subdiagonals
and  KU
       superdiagonals, and X and B are N-by-NRHS matrices

etc. etc.

Ok, je regarde... j'ai trouvé aussi une biblio qui s'appelle SuperLU.

En fait me demande si ce sont des méthodes adaptées à etre insérées au sein d'une simu ou plutot si ce sont des méthodes destinées a résoudre un systeme linéaire une fois pour toute ?

Salut!
J'opterais plutot pour du Cholesky (si ta matrice est bien sym. def. positive) qui préserve le profil de ta matrice avec une renumérotation afin de rendre ta matrice la plus diagonale possible (genre Cuthill MacKee) - et evidemment faire le tout en stockage creux. Si ton second membre change tu pourras conserver ta décomposition telle quelle et ne faire que la descente-remontée pour la résolution du système (s'applique aussi au LU - mais il ne conserve pas le profil donc ton stockage deviendra dense).

Peut-on résoudre ce probleme a l'aide de méthodes spectrales ?

pour les méthodes style GMRES (ou autres méthode de Krylov), il me semble qu'ils nécéssitent le second membre, donc tu ne pourras pas réutiliser les calculs des itérations précédentes. A confirmer.

--
Nicolas Bonneel
http://www.bonneel.com