Bonjour à tous,

Voilà mon problème :

Alors je réalise un projet calcul numérique basé sur la méthode de la relaxation, une amélioration de la méthode de Jacobi (ou Gauss Seidel je ne sais jamais), une méthode itérative permettant de résoudre un systême linéaire.

Je programme ceci en Fortran 77, (pour info), et j'ai un problème. Pour utiliser la méthode de la relaxation, il faut calculer ce qu'on appelle le w optimal. Le w optimal se calcule à partir d'une formule "w=2/(1+sqrt(1-ra*ra))" où ra est le rayon spectral d'une matrice ne contenant que des 0 sur la diagonale.

Pour calculer le rayon spectral j'ai tenté la méthode de la puissance, dont mes tests montrent qu'elle marche sur une matrice qui ne contient pas que des 0 sur la diagonale, mais malheuresement pas sur la matrice qu'il me faut.

Aussi, je viens demander votre aide, au cas où vous connaitriez une façon de calculer le rayon spectral d'une matrice (donc les couts de calculs ne sont pas en n factoriel), ou alors si vous êtes des spécialistes de la méthode de la relaxation (SOR en anglais)