Solution pour matrice tridiagonale par blocs

Bonjour,
Si quelqu'un pourra m'aider à écrire un code de résolution d'une équation linéaire Ax=b avec A une matrice tridiagonale par blocs telque
A=(A1 A2 0 ... 0 )
(A2 A1 A2 ... 0 )
(... ... ... )
(... A2 A1)

A1 une matrice tri-diagonale symétriques et A2 une matrice diagonale.
Je l'ai déjà fait avec factorisation LU 1D. Mais, pas en 2D.
Merci d'avance.

Qu'as-tu tenté jusqu'à présent ?

Bonjour,
En fait, j'ai utilisé l'algorithme de factorisation LU pour résoudre l'équation en 1D. En passant à 2D, j'ai utilisé la transformation l=i+j*N (NxN: dimension de ma matrice) et j'ai pensé à 2 solutions: soit considérer la matrice complète (creuse) ou tridiagonale à blocs. Par ma question, je voudrai étendre mon premier code à 2D (c-à-d opter pour la deuxième solution).
Merci.