◆ zgetri()

行列の逆行列 (複素行列)

目的: 本ルーチンはzgetrfにより計算されたLU分解を用いて行列の逆行列を求める.

本ルーチンでは, Uの逆行列を求め, inv(A)についての連立方程式 inv(A)*L = inv(U) を解くことによりinv(A)を求める.

引数

[in]	n	行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in]	lda	二次元配列a[][]の整合寸法. (lda >= max(1, n))
[in,out]	a[][]	配列 a[la][lda] (la >= n) [in] zgetrfにより計算された分解 A = PLU のLおよびU. [out] info = 0の場合, 元の行列Aの逆行列.
[in]	ipiv[]	配列 ipiv[lipiv] (lipiv >= n) zgetrfによるピボットインデックス. 1 <= i <= nに対して, 行列の第i行は第ipiv[i-1]行と交換されたことを表す.
[out]	work[]	配列 work[lwork] 作業領域. info = 0の場合, work[0] にlworkの最適値を返す.
[in]	lwork	配列 work[]のサイズ (lwork >= max(1, n)) 最適効率のためには, lwork >= n*nbとする. ただし, nbは最適ブロックサイズである. lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, 最適サイズを求める計算だけを行い, work[0]にその値を返す.
[out]	info	= 0: 正常終了 = -1: 入力パラメータ n の誤り (n < 0) = -2: 入力パラメータ lda の誤り (lda < max(1, n)) = -6: 入力パラメータ lwork の誤り (lworkが小さすぎる) = i > 0: Uのi番目の対角要素が0である. 行列は特異で, その逆行列を求めることはできなかった.