|
|
◆ Zppcon()
| Sub Zppcon |
( |
Uplo As |
String, |
|
|
N As |
Long, |
|
|
Ap() As |
Complex, |
|
|
ANorm As |
Double, |
|
|
RCond As |
Double, |
|
|
Info As |
Long |
|
) |
| |
行列の条件数 (正定値エルミート行列) (圧縮形式)
- 目的
- 本ルーチンはZpptrfにより計算されたコレスキー分解 A = U^H*U または A = L*L^H を用いて圧縮形式の正定値エルミート行列の(1-ノルムによる)条件数の推定を行う.
norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する. RCond = 1 / (norm(A) * norm(inv(A)))
- 引数
-
| [in] | Uplo | 分解形のUあるいはLのどちらが格納されているかを指定.
= "U": コレスキー分解 A = U^T*U の上三角行列 U.
= "L": コレスキー分解 A = L*L^T の下三角行列 L. |
| [in] | N | 行列Aの行および列数. (N >= 0) (N = 0 の場合, RCond = 1 を返す) |
| [in] | Ap() | 配列 Ap(LAp - 1) (LAp >= N(N + 1)/2)
Zpptrfにより計算されたコレスキー分解 A = U^H*U または A = L*L^H の三角行列UまたはL(圧縮形式). |
| [in] | ANorm | 対称行列Aの1ノルム(または無限ノルム). (ANorm >= 0) |
| [out] | RCond | 行列Aの条件数の逆数. 次のように計算する.
RCond = 1/(ANorm * Ainvnm)
ここで, Ainvnmは本ルーチン内で求められたinv(A)の1-ノルムの推定値である. |
| [out] | Info | = 0: 正常終了.
= -1: パラメータ Uplo の誤り. (Uplo <> "U"および"L")
= -2: パラメータ N の誤り. (N < 0)
= -3: パラメータ Ap() の誤り.
= -4: パラメータ ANorm の誤り. (ANorm < 0) |
- 出典
- LAPACK
- 使用例
- Zppsvの使用例を参照せよ.
|
1.9.6