WDgbsv WDgesv WDgtsv WDtrtrs

WDgbsv() Function WDgbsv ( N As  Long, Kl As  Long, Ku As  Long, Ab As  Variant, B As  Variant, Optional Nrhs As  Long = 1  ) 連立一次方程式 AX = B の解 (一般帯行列) 目的 WDgbsvは次の連立一次方程式を解く. A * X = B ここで, Aは下帯幅Kl, 上帯幅KuのN×N帯行列, また, XおよびBはN×Nrhs行列である. 行交換によるピボットの部分選択を行うLU分解を用いて, 次のようにAを分解する. A = L * U ここで, Lは置換行列と下帯幅Klで対角要素が1の下三角行列の積, そして, Uは上帯幅Kl+Kuの上三角行列である. 次に, 分解されたAを用いて連立方程式 A * X = B の解を求める. 戻り値 N+2 × Nrhs 列1 列2 ・・・ 列Nrhs 行1〜N 解行列 X 行N+1 1/条件数 0 ・・・ 0 行N+2 リターンコード 0 ・・・ 0 リターンコード. = 0: 正常終了. = i > 0: 行列のi番目のピボットがゼロになった. (行列 A は特異) 引数 [in] N 連立方程式の数, すなわち, 行列Aの行および列数. (N >= 1) [in] Kl 行列Aの下帯幅(対角要素を除く). (Kl >= 0) [in] Ku 行列Aの上帯幅(対角要素を除く). (Kl >= 0) [in] Ab (Kl+1+Ku×N) N×N係数行列 A. (帯格納形式. 詳細は下記参照) [in] B (N × Nrhs) N×Nrhs右辺行列 B. [in] Nrhs (省略可) 右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) 詳細 次の例は, N = 6, Kl = 2, Ku = 1 の場合の帯行列形式を表す. * a12 a23 a34 a45 a56 a11 a22 a33 a44 a55 a66 a21 a32 a43 a54 a65 * a31 a42 a53 a64 * * *で示されたセルは使用されない. 出典 LAPACK 使用例 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, ( 2.34 0.57 0 ) ( 0.7416 ) A = ( 0.65 1.98 -1.39 ), B = ( 0.7885 ) ( 0 1.50 1.73 ) ( 1.0833 ) とする.