XLPack 7.0
XLPack 数値計算ライブラリ (C API) リファレンスマニュアル
読み取り中…
検索中…
一致する文字列を見つけられません

◆ 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]norm1ノルム条件数を求めるか, 無限ノルム条件数を求めるかを指定.
= '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]anormnorm = '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