dspcon dsptrf dsptri dsptrs dsycon dsytrf dsytri dsytrs

dsytri() void dsytri ( char  uplo, int  n, int  lda, double  a[], int  ipiv[], double  work[], int *  info  ) 行列の逆行列 (対称行列) 目的 本ルーチンはdsytrfにより計算された分解 A = U*D*U^T または A = L*D*L^T を用いて, 対称行列Aの逆行列を求める. 引数 [in] uplo 分解形が上三角行列あるいは下三角行列のどちらで格納されているかを指定. = 'U': 上三角, 分解形は A = U*D*U^T. = 'L': 下三角, 分解形は A = L*D*L^T. [in] n 行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] lda 二次元配列a[][]の整合寸法. (lda >= max(1, n)) [in,out] a[][] 配列 a[la][lda] (la >= n) [in] dsytrfにより求められたUまたはLを得るために使われるブロック対角行列Dおよび乗数. [out] info = 0の場合, 元の行列Aの(対称な)逆行列. uploに従い, 逆行列の上あるいは下三角部分が出力され, a[][]の対角要素より下あるいは上の残りの部分は参照されない. [in] ipiv[] 配列 ipiv[lipiv] (lipiv >= n) dsytrfにより求められた行および列の交換とDのブロック構造の情報. [out] work[] 配列 work[lwork] (lwork >= n) 作業領域. [out] info = 0: 正常終了 = -1: 入力パラメータ uplo の誤り (uplo != 'U'および'L') = -2: 入力パラメータ n の誤り (n < 0) = -3: 入力パラメータ lda の誤り (lda < max(1, n)) = i > 0: Dのi番目の要素が0である. 行列は特異で, その逆行列を求めることはできなかった. 出典 LAPACK