pdf 224167 bytes - Cimat

Gradiente biconjugado precondicionado con SAI
Hay que hacer cuatro multiplicaciones matrix-vector: p k A, A p k, r k M, M r k. Para hacerlo
eficiente, hay que almacenar A y M dos veces cada una usando CSR y CSC.
x0, coordenada inicial
r0 b− A x0, gradiente inicial
T
r 0 r 0,
pseudo-gradiente inicial
p0 M r0, dirección de descenso inicial
p 0 r 0 M, pseudo-dirección de descenso
, tolerancia
k 0
mientras ∥rk∥ α k ← ̃r k M r k
̃p k A p k
x k1 x k k p k
r k1 r k− A p k
r k1 r k− p k A
k r k1 M r k1
r k M r k
p k1 M r k1 k1 p k
p k1 r k1 M k1 p k
k k1
11/10/11 14/27