<< luget Sparse Decompositions spchol >>

Scilab Help >> Sparses Matrix > Sparse Decompositions > lusolve

lusolve

疎な線形システムの解を得る

呼び出し手順

x=lusolve(hand,b)

x=lusolve(A,b)

引数

b

通常の実数行列

A

可逆な実数正方疎行列

hand

疎行列LU分解を計算した際のハンドル(lufactの出力)

x

通常の実数行列

説明

x=lusolve(hand,b) は疎な線形システム A*x = bを解きます.

[hand,rk]=lufact(A) は lufact の出力です.

x=lusolve(A,b)は,疎な線形システム A*x = bを解きます.

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

参照


Report an issue
<< luget Sparse Decompositions spchol >>