|
|
◆ _ztrcon()
| void _ztrcon |
( |
char |
norm, |
|
|
char |
uplo, |
|
|
char |
diag, |
|
|
int |
n, |
|
|
int |
lda, |
|
|
doublecomplex |
a[], |
|
|
double * |
rcond, |
|
|
doublecomplex |
work[], |
|
|
double |
rwork[], |
|
|
int * |
info |
|
) |
| |
行列の条件数 (複素三角行列)
- 目的
- 本ルーチンは三角行列Aの(1ノルムまたは無限ノルムによる)条件数の逆数を推定する.
Aのノルムを計算し, norm(inv(A))の推定値を求め, 次のように条件数の逆数を計算する. rcond = 1 / (norm(A) * norm(inv(A)))
- 引数
-
| [in] | norm | 1ノルム条件数を求めるか, 無限ノルム条件数を求めるかを指定.
= '1'または'O': 1ノルム.
= 'I': 無限ノルム. |
| [in] | uplo | = 'U': Aは上三角行列.
= 'L': Aは下三角行列. |
| [in] | diag | = 'N': Aは単位三角行列ではない.
= 'U': Aは単位三角行列である. (a[][]の対角要素を参照せず, 1とみなす) |
| [in] | n | 行列Aの行および列数. (n >= 0) (n = 0 の場合, rcond = 1 を返す) |
| [in] | lda | 二次元配列a[][]の整合寸法. (lda >= max(1, n)) |
| [in] | a[][] | Array a[la][lda] (la >= n)
n×n三角行列 A. uploに従い上三角部分あるいは下三角部分のみ参照される. |
| [out] | rcond | 行列Aの条件数の逆数. (rcond = 1/(norm(A) * norm(inv(A)))) |
| [out] | work[] | 配列 work[lwork] (lwork >= 2*n)
作業領域. |
| [out] | rwork[] | 配列 rwork[lrwork] (lrwork >= n)
作業領域. |
| [out] | info | = 0: 正常終了
= -1: 入力パラメータ norm の誤り (norm != '1', 'O'および'I')
= -2: 入力パラメータ uplo の誤り (uplo != 'U'および'L')
= -3: 入力パラメータ diag の誤り (diag != 'N'および'U')
= -4: 入力パラメータ n の誤り (n < 0)
= -5: 入力パラメータ lda の誤り (lda < max(1, n)) |
- 出典
- LAPACK
|
1.9.6