Bonjour,

Je cherche à résoudre un système linaire de forme Ax=b avec A une matrice creuse, en utilisant la décomposition de Cholesky.
J'ai déjà codé une subroutine qui me donne une matrice creuse avec un stockage CSR ainsi qu'une subroutine qui effectue la décomposition de Cholesky A=L*tL mais avec une matrice pleine.
Je me demandais donc comment obtenir la matrice creuse L à partir de la matrice creuse A.

Merci d'avance.