pdf 224167 bytes - Cimat

El algoritmo para calcular las entradas de ̃G es el siguiente:
Sea S L la estructura de G
para i ← 1…n
para ∀ (i , j )∈S L
resolver ( A ̃G i) j =δ i j
fin_para
fin_para
Inversa aproximada rala factorizada
Resolver ( A ̃G i) j =δ i j significa que, si m=η( ̃G i) el número de entradas no cero del renglón i de
̃G, entonces hay que resolver un sistema SPD local pequeño de tamaño m×m.
( A ̃G i) j =δ i j, (i , j)∈S L,
Ejemplo. Si la estructura del renglón 9 de ̃G contiene 4, 7, 8, 9, entonces, para calcular los
valores de ̃G 9 hay que resolver el sistema:
(a 4 4 a4 7 a4 8 a4 9 ̃g
99)( a7 4 a7 7 a7 8 a79 a8 4 a8 7 a8 8 a89 a9 4 a9 7 a9 8 a
9 4
̃g 9 7 0
̃g 0 98
̃g 99)=(0
1).
11/10/11 18/27