dspcon dsptrf dsptri dsptrs dsycon dsytrf dsytri dsytrs

dsptri() void dsptri ( char  uplo, int  n, double  ap[], int  ipiv[], double  work[], int *  info  ) 行列の逆行列 (対称行列) (圧縮形式) 目的 本ルーチンはdsptrfにより計算された分解 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,out] ap[] 配列 ap[lap] (lap >= n(n + 1)/2) [in] dsptrfにより求められた, UまたはLを得るために使われるブロック対角行列Dおよび乗数. 圧縮形式で三角行列として格納する. [out] info = 0の場合, 圧縮形式で三角行列として格納された元の行列の(対称な)逆行列. uploに従い, 逆行列の上あるいは下三角部分が格納される. [in] ipiv[] 配列 ipiv[lipiv] (lipiv >= n) dsptrfにより求められた行および列の交換とDのブロック構造の情報. [out] work[] 配列 work[lwork] (lwork >= n) 作業領域. [out] info = 0: 正常終了 = -1: 入力パラメータ uplo の誤り (uplo != 'U'および'L') = -2: 入力パラメータ n の誤り (n < 0) = i > 0: Dのi番目の要素が0である. 行列は特異で, その逆行列を求めることはできなかった. 出典 LAPACK