WDpbsv WDposv WDptsv WDsysv

WDpbsv() Function WDpbsv ( Uplo As  String, N As  Long, Kd As  Long, Ab As  Variant, B As  Variant, Optional Nrhs As  Long = 1  ) 連立一次方程式 AX = B の解 (正定値対称帯行列) 目的 WDpbsvは次の連立一次方程式を解く. A * X = B ここで, AはN×N正定値対称帯行列, また, XおよびBはN×Nrhs行列である. まず, コレスキー分解によりAを次のように分解する. A = U^T*U (Uplo = "U"の場合) A = L*L^T (Uplo = "L"の場合) ここで, UはAと同じ上帯幅の上三角行列, LはAと同じ下帯幅の下三角行列である. 次に, 分解された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主小行列が正定値でないため分解を完了できなかった. 引数 [in] Uplo = "U": Aの上三角部分を格納. = "L": Aの下三角部分を格納. [in] N 連立方程式の数, すなわち, 行列Aの行および列数 (N >= 1) [in] Kd 上または下帯幅. (対角要素を含まない) (Kd >= 0) [in] Ab (Kd+1 × N) N×N係数行列 A. (対称帯行列形式. 詳細は下記参照) [in] B (N × Nrhs) N×Nrhs右辺行列 B. [in] Nrhs (省略可) 右辺行列Bの列数. (Nrhs >= 1) (省略時 = 1) 詳細 次の例は, n = 6, kd = 2, Uplo = "U" の場合の対称帯行列形式を表す. * * a13 a24 a35 a46 * a12 a23 a34 a45 a56 a11 a22 a33 a44 a55 a66 同様に uplo = "L" の場合, Aの形式は次のようになる. a11 a22 a33 a44 a55 a66 a21 a32 a43 a54 a65 * a31 a42 a53 a64 * * *で示された配列要素は使用されない. 出典 LAPACK 使用例 連立一次方程式 Ax = B を解く. また, Aの条件数の逆数の推定値(RCond)を求める. ただし, Aは正定値対称帯行列で ( 0.61 0.79 0 ) ( 0.3034 ) A = ( 0.79 2.23 0.25 ), B = ( 0.8537 ) ( 0 0.25 2.87 ) ( 0.8000 ) とする.