◆ Zptcon()

行列の条件数 (正定値エルミート3重対角行列)

目的: 本ルーチンはZpttrfにより計算された, 分解 A = L*D*L^H または A = U^H*D*U を用いて正定値エルミート3重対角行列の(1-ノルムによる)条件数の逆数を求める.

norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する.
RCond = 1 / (norm(A) * norm(inv(A)))

引数

[in]	N	行列Aの行および列数. (N >= 0) (N = 0 の場合, RCond = 1 を返す)
[in]	D()	配列 D(LD - 1) (LD >= N) Zpttrfにより計算された, Aの分解結果の対角行列DのN個の対角要素.
[in]	E()	配列 E(LE - 1) (LE >= N - 1) Zpttrfにより計算された, Aの分解結果の2重対角行列UまたはLのN-1個の副対角要素.
[in]	ANorm	分解前の行列Aの1ノルム. (ANorm >= 0)
[out]	RCond	行列Aの条件数の逆数. 次のように計算する. RCond = 1/(ANorm * Ainvnm) ここで, Ainvnmは本ルーチンで求められたinv(A)の1-ノルムの値である.
[out]	Info	= 0: 正常終了. = -1: パラメータ N の誤り. (N < 0) = -2: パラメータ D() の誤り. = -3: パラメータ E() の誤り. = -4: パラメータ ANorm の誤り. (ANorm < 0)