Dpocon Dpotrf Dpotri Dpotrs Dppcon Dpptrf Dpptri Dpptrs

Dpotri() Sub Dpotri ( Uplo As  String, N As  Long, A() As  Double, Info As  Long  ) 行列の逆行列 (正定値対称行列) 目的 本ルーチンはDpotrfにより計算されたコレスキー分解 A = U^T*U または A = L*L^T を用いて, 正定値対称行列Aの逆行列を求める. 引数 [in] Uplo = "U": 上三角行列UをA()に格納. = "L": 下三角行列LをA()に格納. [in] N 行列Aの行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in,out] A() 配列 A(LA1 - 1, LA2 - 1) (LA1 >= N, LA2 >= N) [in] Dpotrfにより計算されたコレスキー分解 A = U^T*U または A = L*L^T の三角行列UまたはL. [out] Info = 0の場合, Aの(対称な)逆行列の上あるいは下三角部分. 入力のUまたはLを上書きする. [out] Info = 0: 正常終了. = -1: パラメータ Uplo の誤り. (Uplo <> "U"および"L") = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ A() の誤り. = i > 0: UまたはLのi番目の対角要素が0で, 逆行列を求めることができなかった. 出典 LAPACK 使用例 行列 A の逆行列を求める. ただし, Aは正定値対称行列で ( 2.2 -0.11 -0.32 ) A = ( -0.11 2.93 0.81 ) ( -0.32 0.81 2.37 ) とする. Sub Ex_Dpotri() Const N = 3 Dim A(N - 1, N - 1) As Double, Info As Long A(0, 0) = 2.2 A(1, 0) = -0.11: A(1, 1) = 2.93 A(2, 0) = -0.32: A(2, 1) = 0.81: A(2, 2) = 2.37: Call Dpotrf("L", N, A(), Info) If Info = 0 Then Call Dpotri("L", N, A(), Info) Debug.Print "Inv(A) =" Debug.Print A(0, 0), A(0, 1), A(0, 2) Debug.Print A(1, 0), A(1, 1), A(1, 2) Debug.Print A(2, 0), A(2, 1), A(2, 2) Debug.Print "Info =", Info End Sub 実行結果 Inv(A) = 0.463651335375572 0 0 1.10603849087679E-04 0.376908423331071 0 6.25649106339333E-02 -0.128801869057578 0.474409403096834 Info = 0