◆ zgtcon()

行列の条件数 (複素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] (ldl >= n - 1) zgttrfにより計算されたAのLU分解結果の行列Lを定義するn-1個の乗数.
[in]	d[]	配列 d[ld] (ld >= n) AのLU分解結果の上三角行列Uの対角要素.
[in]	du[]	配列 du[ldu] (ldu >= n - 1) Uのn-1個の上副対角要素.
[in]	du2[]	配列 du2[ldu2] (ldu2 >= n - 2) Uのn-2個の第2上副対角要素.
[out]	ipiv[]	配列 ipiv[lipiv] (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]	work[]	配列 work[lwork] (lwork >= 2*n) 作業領域.
[out]	info	= 0: 正常終了 = -1: 入力パラメータ norm の誤り (norm != '1', 'O'および'I') = -2: 入力パラメータ n の誤り (n < 0) = -8: 入力パラメータ anorm の誤り (anorm < 0)