Conjugate Gradient (2) with preconditioning

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 ) ;

...

• conjugate gradient (problem sizes 481 and 936) with preconditioning (triangular matrix)

Diapositive précédente

Diapositive suivante

Revenir à la première diapositive

Afficher la version graphique