pdf 224167 bytes - Cimat
pdf 224167 bytes - Cimat
pdf 224167 bytes - Cimat
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
El algoritmo es el siguiente<br />
x 0, coordenada inicial<br />
r 0 b− A x 0, gradiente inicial<br />
p 0 M r 0, descenso inicial<br />
k 0<br />
mientras ∥r k∥<br />
αk ← r T<br />
k M r k<br />
T<br />
p k A pk<br />
x k1 x k k p k<br />
r k1 r k− k A p k<br />
k r T<br />
k1 M r k1<br />
T<br />
r k M r k<br />
pk1 M r k1 k pk k k1<br />
M r k se calcula con dos multiplicaciones matriz-triangular*vector,<br />
T<br />
M r k =G l<br />
(G l r k ).<br />
No hay que resolver un sistema de ecuaciones en cada paso.<br />
Inversa aproximada rala factorizada<br />
11/10/11 20/27