Bonjour à toutes et à tous,
Je souhaite inverser une matrice bande de très grande dimension (du style (m*m)**2 où m peut valoir 10^2) ce qui prend avec un algorithme classique du style décomposition LU un temps phénoménal, ce qui est problématique car inclus dans une boucle itérative. La matrice est quelconque (donc non symétrique)
Avant tout, je tiens à préciser, puisque c'est visiblement la source de beaucoup de confusions sur le forum, que je souhaite bien INVERSER la matrice et non résoudre directement un système (l'inversion fait partie d'un problème dont les détails sont je pense inutiles)
L'ennuie est que les libraires que j'ai pu trouver dans LAPACK par exemple ne proposent que des RESOLUTIONS de systèmes et non des inversions matricielles. Je cherche donc un programme capable d'inverser une une telle matrice bande.
Merci d'avance pour vos propositions.
Mathieu
Partager