疎行列LU分解の展開
[P,L,U,Q]=luget(hand)
ハンドル, lufact
の出力
疎交換行列
疎行列, lower triangular if hand
が
正則行列から得られた場合は上三角
正方正則上三角疎行列(主対角項は1)
疎交換行列
疎行列 A
に関して
コマンド[hand,rk]=lufact(A)
により得られた
hand
を指定すると,[P,L,U,Q]=luget(hand)
は,P*L*U*Q=A
となるような4つの疎行列を返します.
行列Aはフルランクである必要はありませんが正方行列である必要があります (Aは疎行列と仮定されるため,正方化するために必要に応じて0を追加することができます).
A
が正則でない場合,
行列L
は(rk
個の
独立した非ゼロ列について)列圧縮されます
正則な疎行列 Q'*inv(U)
は
A
を列圧縮します.