_dsbgv _dsbgvd _dsbgvx _dspgv _dspgvd _dspgvx _dsygv _dsygvd _dsygvx

_dsygvd() void _dsygvd ( int  itype, char  jobz, char  uplo, int  n, int  lda, double  a[], int  ldb, double  b[], double  w[], double  work[], int  lwork, int  iwork[], int  liwork, int *  info  ) (分割統治法ドライバ) 一般化固有値問題 (対称行列) 目的 本ルーチンは実対称行列の一般化固有値問題 Ax = λBx, ABx = λx または BAx = λx のすべての固有値, および, 必要により固有ベクトルを求める. ここで, A と B は対称行列で, さらに B は正定値である. 固有ベクトルも求める場合, 分割統治法を使用する. 引数 [in] itype 解くべき問題のタイプを指定. = 1: Ax = λBx. = 2: ABx = λx. = 3: BAx = λx. [in] jobz = 'N': 固有値のみ求める. = 'V': 固有値と固有ベクトルを求める. [in] uplo = 'U': A および B の上三角部分を格納する. = 'L': A および B の下三角部分を格納する. [in] n 行列 A および B の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] lda 二次元配列 a[][] の整合寸法. (lda >= max(1, n)) [in,out] a[][] 配列 a[la][lda] (la >= n) [in] 対称行列 A. uplo = 'U' の場合, a[][] の n x n 上三角部分に行列 A の上三角部分を格納する. uplo = 'L' の場合, a[][] の n x n 下三角部分に行列 A の下三角部分を格納する. [out] jobz = 'V': info = 0 の場合, 次のように正規化された固有ベクトルからなる行列 Z を a[][] に返す.   itype = 1 または 2: Z^T*B*Z = I   itype = 3: Z^T*inv(B)*Z = I jobz = 'N': a[][] の上三角部分(uplo = 'U' の場合)あるいは下三角部分(uplo = 'L' の場合)は, 対角部分を含め壊される. [in] ldb 二次元配列 b[][] の整合寸法. (ldb >= max(1, n)) [in,out] b[][] 配列 b[lb][ldb] (lb >= n) [in] 正定値対称行列 B. uplo = 'U' の場合, b[][] の n x n 上三角部分に行列 B の上三角部分を格納する. uplo = 'L' の場合, b[][] の n x n 下三角部分に行列 B の下三角部分を格納する. [out] info <= n の場合, b[][] の行列を格納している部分はコレスキー分解 B = U^T*U あるいは B = L*L^T の三角行列 U あるいは L により上書きされる. [out] w[] 配列 w[lw] (lw >= n) info = 0 の場合, 求められた固有値(昇順). [out] work[] 配列 work[lwork] 作業領域. info = 0 の場合, work[0] に lwork の最適値を返す. [in] lwork 配列 work[] のサイズ (lwork >= 1 (n <= 1 の場合), 2*n + 1 (jobz = 'N' の場合), 2*n^2 + 6*n + 1 (jobz = 'V' の場合)) lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, work[] と iwork[] の最適サイズを求める計算だけを行い, work[0] と iwork[0] にそれらの値を返す. [out] iwork[] 配列 iwork[liwork] 整数作業領域. info = 0の場合, iwork[0] に liwork の最適値を返す. [in] liwork 配列 iwork[] のサイズ (liwork >= 1 (n <= 1 の場合), 1 (jobz = 'N' の場合), 5*n + 3 (jobz = 'V' の場合)) liwork = -1 の場合, 作業領域サイズの問い合わせとみなし, work[] と iwork[] の最適サイズを求める計算だけを行い, work[0] と iwork[0] にそれらの値を返す. [out] info = 0: 正常終了 = -1: 入力パラメータ itype の誤り (itype < 1 または itype > 3) = -2: 入力パラメータ jobz の誤り (jobz != 'V' および 'N') = -3: 入力パラメータ uplo の誤り (uplo != 'U' および 'L') = -4: 入力パラメータ n の誤り (n < 0) = -5: 入力パラメータ lda の誤り (lda < max(1, n)) = -7: 入力パラメータ ldb の誤り (ldb < max(1, n)) = -11: 入力パラメータ lwork の誤り (lworkが小さすぎる) = -13: 入力パラメータ liwork の誤り (liworkが小さすぎる) = i (0 < i <= n): jobz = 'N' の場合, アルゴリズムが収束しなかった. 中間の3重対角形の非対角要素のうち i 個が 0 に収束しなかった. jobz = 'V' の場合, 行および列 i/(n + 1) 〜 mod(i, n + 1) にある小行列において固有値を求める事ができなかった. = i (i > n): B の i-n 次小行列が正定値でない. B の分解が完了できず, 固有値・固有ベクトルは計算されなかった. 出典 LAPACK