Fortran 90 for PCG

iteration_loop: DO

Ö

u = MATMUL(A,p)

up = DOT_PRODUCT(r,r)

alpha = up/DOT_PRODUCT(p,u)

xnew = x + p*alpha

r = r ñ u*alpha

beta = DOT_PRODUCT(r,r)/up

p = r + p*beta

Ö

END DO iteration_loop

Previous slide

Next slide

Back to first slide

View graphic version