extração dos fatores LU esparsos
[P,L,U,Q]=luget(hand)
manipulador, saída de lufact
matriz de permutação esparsa
matriz esparsa, triangular infeiror se hand
é obtida de uma matriz não-singular
matriz quadrada triangular superior esparsa não-snigular preenchida com 1 ao longo da diagonal principal
matriz de permutação esparsa
[P,L,U,Q]=luget(hand)
com hand
obtida pelo comando [hand,rk]=lufact(A)
com
A
uma matriz esparsa retorna quatro matrizes esparsas
tais que P*L*U*Q=A
.
A matriz A não precisa ser de posto cheio, mas deve ser quadrada (desde que A é assumida esparsa, pode-se adicionar 0, se necessário, para quadrar A).
Se A
é singular, a matriz L
é
de colunas comprimidas (com rk
colunas independentes
não-nulas): a matriz não-singular esparsa Q'*inv(U)
comprime em colunas A
.