Salut,
je suis en train d'implementer une decomposition LU et Gauss-Seidel. Mon probleme est de gerer ces algorithme lorsqu il y a des 0 sur la diagonale. Pour Gauss-Seidel, je pense qu'il n y a pas vraiment de solution, il faut tout simplement utilise rautre chose, mais quoi ? (pivot gauss...mais probleme de stabilite ?)
Pour la decomposition LU, par contre je crois qu on peut intervertir des colonnes lorsqu on rencontre un 0, mais je ne sais pas encore comment on gere ca.
Mais ca je pense que ca marche lorsqu on rencontre seulement quelques 0, mais je voulais appliquer cet algorithme pour trouver les valeurs propres de la matrice compagnon d'un polynome, voir le lien suivant : (http://fr.wikipedia.org/wiki/Matrice_compagnon)
et cette matrice ne comporte que des 0 sur la diagonale (enfin presque sauf sur la derniere colonne).
Je voudrais donc avoir votre avis sur la maniere dont on peut gerer ce genre de probleme.
Merci.
Partager