|
|
◆ Zgtcon()
| Sub Zgtcon |
( |
Norm As |
String, |
|
|
N As |
Long, |
|
|
Dl() As |
Complex, |
|
|
D() As |
Complex, |
|
|
Du() As |
Complex, |
|
|
Du2() As |
Complex, |
|
|
IPiv() As |
Long, |
|
|
ANorm As |
Double, |
|
|
RCond As |
Double, |
|
|
Info As |
Long |
|
) |
| |
行列の条件数 (複素3重対角行列)
- 目的
- 本ルーチンは3重対角行列Aの条件数の逆数を推定する. 計算にはZgttrfにより求められた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] | Dl() | 配列 Dl(LDl - 1) (LDl >= N - 1)
Zgttrfにより計算された A のLU分解結果の行列 L を定義するN-1個の乗数. |
| [in] | D() | 配列 D(LD - 1) (LD >= N)
A のLU分解結果の上三角行列 U の対角要素. |
| [in] | Du() | 配列 Du(LDu - 1) (LDu >= N - 1)
U のN-1個の上副対角要素. |
| [in] | Du2() | 配列 Du2(LDu2 - 1) (LDu2 >= N - 2)
U のN-2個の第2上副対角要素. |
| [in] | IPiv() | 配列 IPiv(LIPiv - 1) (LIPiv >= N)
ピボットインデックス. 1 <= i <= N について, 行列の第i行は第IPiv(i-1)行と交換されたことを表す. IPiv(i-1)は常にiまたはi+1である. IPiv(i-1) = i は行の交換が不要であったことを示す. |
| [in] | ANorm | Norm = "1"または"O": 分解前の行列Aの1ノルム (ANorm >= 0)
Norm = "I": 分解前の行列Aの無限ノルム (ANorm >= 0) |
| [out] | RCond | 行列Aの条件数の逆数. 次のように計算する.
RCond = 1/(ANorm * Ainvnm)
ここで, Ainvnmは本ルーチンで求められたinv(A)の1-ノルムの推定値である. |
| [out] | Info | = 0: 正常終了.
= -1: パラメータ Norm の誤り. (Norm <> "1", "O"および"I")
= -2: パラメータ N の誤り. (N < 0)
= -3: パラメータ Dl() の誤り.
= -4: パラメータ D() の誤り.
= -5: パラメータ Du() の誤り.
= -6: パラメータ Du2() の誤り.
= -7: パラメータ IPiv() の誤り.
= -8: パラメータ ANorm の誤り. (ANorm < 0) |
- 出典
- LAPACK
- 使用例
- Zgtsvの使用例を参照せよ.
|
1.9.6