|
|
◆ Zgecon()
| Sub Zgecon |
( |
Norm As |
String, |
|
|
N As |
Long, |
|
|
A() As |
Complex, |
|
|
ANorm As |
Double, |
|
|
RCond As |
Double, |
|
|
Info As |
Long |
|
) |
| |
行列の条件数 (複素行列)
- 目的
- 本ルーチンは複素行列Aの(1ノルムまたは無限ノルムによる)条件数の逆数を推定する. 計算にはZgetrf, Zgesv または Zgesvxにより求めたLU分解が使用される.
norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する. RCond = 1 / (norm(A) * norm(inv(A)))
- 引数
-
| [in] | Norm | 1ノルム条件数を求めるか, 無限ノルム条件数を求めるかを指定.
= "1"または"O": 1ノルム.
= "I": 無限ノルム. |
| [in] | N | 行列Aの行および列数. (N >= 0) (N = 0 の場合, RCond = 1 を返す) |
| [in] | A() | 配列 A(LA1 - 1, LA2 - 1) (LA1 >= N, LA2 >= N)
Zgetrf, Zgesv または Zgesvxにより計算された分解 A = P*L*U のLおよびU. |
| [in] | ANorm | Norm = "1"または"O": 分解前の行列Aの1ノルム. (ANorm >= 0)
Norm = "I": 分解前の行列Aの無限ノルム. (ANorm >= 0) |
| [out] | RCond | 行列Aの条件数の逆数. (RCond = 1/(norm(A) * norm(inv(A)))) |
| [out] | Info | = 0: 正常終了.
= -1: パラメータ Norm の誤り. (Norm <> "1", "O"および"I")
= -2: パラメータ N の誤り. (N < 0)
= -3: パラメータ A() の誤り.
= -4: パラメータ ANorm の誤り. (ANorm < 0) |
- 出典
- LAPACK
- 使用例
- Zgesvの使用例を参照せよ.
|
1.9.6