Gradient conjugué avec préconditionnement

for iter = 1:max_it

z = M \ r ;

rho = (r'*z) ;

if ( iter > 1 ),

beta = rho/rho_1;

p = z + beta*p;

else

p = z;

end

q = A*p;

alpha = rho / (p'*q ) ;

...

• algorithme du gradient conjugué (tailles de problème 481 et 936) avec préconditionnement (matrice triangulaire)

Diapositive précédente

Diapositive suivante

Revenir à la première diapositive

Afficher la version graphique