lusolve

Parâmetros

b

matriz de reais completa

A

matriz quadrada de reais esparsa e invertível

hand

manipulador para fatores de lu esparsos previamente computados (saída de lufact)

x

matriz de reais completa

Descrição

x=lusolve(hand,b) resolve o sistema linear esparso A*x = b.

[hand,rk]=lufact(A) é a saída de lufact.

x=lusolve(A,b) resolve o sistema linear esparso A*x = b

Exemplos

non_zeros=[1,2,3,4];rows_cols=[1,1;2,2;3,3;4,4];
sp=sparse(rows_cols,non_zeros);
[h,rk]=lufact(sp);x=lusolve(h,[1;1;1;1]);ludel(h)
rk,sp*x

non_zeros=[1,2,3,4];rows_cols=[1,1;2,2;3,3;4,4];
sp=sparse(rows_cols,non_zeros);
x=lusolve(sp,-ones(4,1));
sp*x

Ver Também

• sparse
• lufact
• slash
• backslash