dgbsv dgbsvx dgesv dgesvx dgtsv dgtsvx dsgesv

dgesv() void dgesv ( int  n, int  nrhs, int  lda, double  a[], int  ipiv[], int  ldb, double  b[], int *  info  ) (シンプルドライバ) 連立一次方程式 AX = B の解 (一般行列) 目的 本ルーチンは次の連立一次方程式を解く. A * X = B ここで, Aはn×n行列, また, XおよびBはn×nrhs行列である. まず, 行交換によるピボットの部分選択を行うLU分解を用いて, 次のようにAを分解する. A = P * L * U ここで, Pは置換行列, Lは対角要素が1の下三角行列, そして, Uは上三角行列である. 次に, 分解されたAを用いて連立方程式 A * X = B の解を求める. 引数 [in] n 連立方程式の数, すなわち, 行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] nrhs 右辺の数, すなわち, 行列Bの列数. (nrhs >= 0) (nrhs = 0 の場合, 処理を行わずに戻る) [in] lda 二次元配列a[][]の整合寸法. (lda >= max(1, n)) [in,out] a[][] 配列 a[la][lda] (la >= n) [in] n×n 係数行列 A. [out] 分解 A = P*L*U のLおよびU. Lの対角要素(= 1)は格納されない. [out] ipiv[] 配列 ipiv[lipiv] (lipiv >= n) 置換行列Pを定義するピボットインデックス. 第i行が第ipiv[i-1]行と交換されたことを表す. [in] ldb 二次元配列b[][]の整合寸法. (ldb >= max(1, n)) [in,out] b[][] 配列 b[lb][ldb] (lb >= nrhs) [in] n×nrhs右辺行列 B. [out] info = 0 の場合, n×nrhs解行列 X. [out] info = 0: 正常終了 = -1: 入力パラメータ n の誤り (n < 0) = -2: 入力パラメータ nrhs の誤り (nrhs < 0) = -3: 入力パラメータ lda の誤り (lda < max(1, n)) = -6: 入力パラメータ ldb の誤り (ldb < max(1, n)) = i > 0: Uのi番目の対角要素が0である. 分解を完了したが, Uが特異であるため解を計算できなかった. 出典 LAPACK