XLPack 7.0
XLPack 数値計算ライブラリ (C API) リファレンスマニュアル
読み取り中…
検索中…
一致する文字列を見つけられません

◆ zpttrf()

void zpttrf ( int  n,
double  d[],
doublecomplex  e[],
int *  info 
)

係数行列のLDLH分解 (正定値エルミート3重対角行列)

目的
本ルーチンは正定値エルミート3重対角行列AのL*D*L^H分解を計算する. この分解はまた A = U^H*D*U の形式とみなすこともできる.
引数
[in]n行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in,out]d[]配列 d[ld] (ld >= n)
[in] 正定値エルミート3重対角行列Aのn個の対角要素.
[out] AのL*D*L^H分解の対角行列Dのn個の対角要素.
[in,out]e[]配列 e[le] (le >= n - 1)
[in] 正定値エルミート3重対角行列Aのn-1個の副対角要素.
[out] 分解 A = L*D*L^H の下2重対角行列Lのn-1個の下副対角要素. これは, 分解 A = U^H*D*U の上2重対角行列Uの上副対角要素とみなすこともできる. (LおよびUの対角要素は1)
[out]info= 0: 正常終了
= -1: 入力パラメータ n の誤り (n < 0)
= i > 0: i×i首座小行列が正定値でない. i < nの場合, 分解を完了できなかった. i = n の場合, 分解は完了したが d[n-1] <= 0 である.
出典
LAPACK