luget

Parâmetros

hand

manipulador, saída de lufact

P

matriz de permutação esparsa

L

matriz esparsa, triangular infeiror se hand é obtida de uma matriz não-singular

U

matriz quadrada triangular superior esparsa não-snigular preenchida com 1 ao longo da diagonal principal

Q

matriz de permutação esparsa

Descrição

[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.

Exemplos

a=rand(5,2)*rand(2,5);A=sparse(a);
[hand,rk]=lufact(A);[P,L,U,Q]=luget(hand);
full(L), P*L*U*Q-A
clean(P*L*U*Q-A)
ludel(hand)

Ver Também

• sparse
• lusolve
• luget
• clean