|
|
◆ ztbcon()
| void ztbcon |
( |
char |
norm, |
|
|
char |
uplo, |
|
|
char |
diag, |
|
|
int |
n, |
|
|
int |
kd, |
|
|
int |
ldab, |
|
|
doublecomplex |
ab[], |
|
|
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は単位三角行列である. (ab[][]の対角要素を参照せず, 1とみなす) |
| [in] | n | 行列Aの行および列数. (n >= 0) (n = 0 の場合, rcond = 1 を返す) |
| [in] | kd | 三角帯行列Aの上または下帯幅. (kd >= 0) |
| [in] | ldab | 二次元配列ab[][]の整合寸法. (ldab >= kd + 1) |
| [in] | ab[][] | 配列 ab[lab][ldab] (lab >= n)
kd+1×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: 入力パラメータ kd の誤り (kd < 0)
= -6: 入力パラメータ ldab の誤り (ldab < kd+1) |
- 出典
- LAPACK
|
1.9.7