10.05.2013 Views

pdf 224167 bytes - Cimat

pdf 224167 bytes - Cimat

pdf 224167 bytes - Cimat

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!