# luget

疎行列LU分解の展開

## 呼び出し手順

    [P,L,U,Q]=luget(hand)

## パラメータ

hand

  ハンドル, `lufact`の出力

P

  疎交換行列

L

  疎行列, lower triangular if `hand`が 正則行列から得られた場合は上三角

U

  正方正則上三角疎行列(主対角項は1)

Q

  疎交換行列

## 説明

疎行列 `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`を列圧縮します.

## 例

    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)

## 参照

- sparse
- lusolve
- luget
- clean