|
|
◆ zgecon()
| void zgecon |
( |
char |
norm, |
|
|
int |
n, |
|
|
int |
lda, |
|
|
doublecomplex |
a[], |
|
|
double |
anorm, |
|
|
double * |
rcond, |
|
|
doublecomplex |
work[], |
|
|
double |
rwork[], |
|
|
int * |
info |
|
) |
| |
行列の条件数 (複素行列)
- 目的
- 本ルーチンは複素行列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] | lda | 二次元配列a[][]の整合寸法. (lda >= max(1, n)) |
| [in] | a[][] | 配列 a[la][lda] (la >= 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] | work[] | 配列 work[lwork] (lwork >= 2*n)
作業領域. |
| [out] | rwork[] | 配列 rwork[lrwork] (lrwork >= 2*n)
作業領域. |
| [out] | info | = 0: 正常終了
= -1: 入力パラメータ norm の誤り (norm != '1', 'O'および'I')
= -2: 入力パラメータ n の誤り (n < 0)
= -3: 入力パラメータ lda の誤り (lda < max(1, n))
= -5: 入力パラメータ anorm の誤り (anorm < 0) |
- 出典
- LAPACK
|