zgges zgges_r zggesx zggesx_r zggev zggevx

zggev() void zggev ( char  jobvl, char  jobvr, int  n, int  lda, doublecomplex  a[], int  ldb, doublecomplex  b[], doublecomplex  alpha[], doublecomplex  beta[], int  ldvl, doublecomplex  vl[], int  ldvr, doublecomplex  vr[], doublecomplex  work[], int  lwork, double  rwork[], int *  info  ) (シンプルドライバ) 一般化固有値問題 (複素行列) 目的 本ルーチンは n x n 複素行列のペア (A, B) の一般化固有値, および, 必要により左および／または右一般化固有ベクトルを求める. 行列のペア (A, B) の一般化固有値は, A - λB が特異となるようなスカラーλあるいは比 α/β = λ である. β = 0 あるいは両方共 0 の場合に適当な解釈ができるため, これは通常ペア (α, β) で表される. (A, B) の一般化固有値 λ(j) に対応する右一般化固有ベクトル v(j) は次式を満たす. A * v(j) = λ(j) * B * v(j) (A, B) の一般化固有値 λ(j) に対応する左一般化固有ベクトル u(j) は次式を満たす. u(j)^H * A = λ(j) * u(j)^H * B ここで, u(j)^H は u(j) の共役転置を表す. 引数 [in] jobvl = 'N': 左一般化固有ベクトルを求めない. = 'V': 左一般化固有ベクトルを求める. [in] jobvr = 'N': 右一般化固有ベクトルを求めない. = 'V': 右一般化固有ベクトルを求める. [in] n 行列 A, B, VL および VR の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] lda 二次元配列 a[][] の整合寸法. (lda >= max(1, n)) [in,out] a[][] 配列 a[la][lda] (la >= n) [in] ペア (A, B) の行列 A. [out] a[][] は上書きされる. [in] ldb 二次元配列 b[][] の整合寸法. (ldb >= max(1, n)) [in,out] b[][] 配列 b[lb][ldb] (lb >= n) [in] ペア (A, B) の行列 B. [out] b[][] は上書きされる. [out] alpha[] 配列 alpha[lalpha] (lalpha >= n) [out] beta[] 配列 beta[lbeta] (lbeta >= n) alpha[j]/beta[j] (j = 0〜n-1) は一般化固有値である. 注: 商 alpha[j]/beta[j] はオーバーフローあるいはアンダーフローしやすく, beta[j] が 0 になることさえある. 従って, 比 α/β を単純に計算することは避けなければならない. しかし, alpha は常に norm(A) より小さく, 通常は同程度の大きさである. また, beta は常に norm(B) より小さく, 通常は同程度の大きさである. [in] ldvl 二次元配列 vl[][] の整合寸法. (ldvl >= 1 (jobvl = 'N'), ldvl >= n (jobvl = 'V')) [out] vl[][] 配列 vl[lvl][ldvl] (lvl >= n) jobvl = 'V': 左一般化固有ベクトル u(j) が固有値と同じ順に vl[][] の列に入る. 固有ベクトルは, 最大要素が |実数部| + |虚数部| = 1 となるように正規化される. jobvl = 'N': 参照されない. [in] ldvr 二次元配列 vr[][] の整合寸法. (ldvr >= 1 (jobvr = 'N'), ldvr >= n (jobvr = 'V')) [out] vr[][] 配列 vr[lvr][ldvr] (lvr >= n) jobvr = 'V': 右一般化固有ベクトル v(j) が固有値と同じ順に vr[][] の列に入る. 固有ベクトルは, 最大要素が |実数部| + |虚数部| = 1 となるように正規化される. jobvr = 'N': 参照されない. [out] work[] 配列 work[lwork] 複素数作業領域. info = 0 の場合, work[0] に lwork の最適値を返す. [in] lwork 配列 work[]のサイズ. (lwork >= max(1, 2*n)) 一般にパフォーマンスをよくするためには lwork を大きくしなければならない. lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, 最適サイズを求める計算だけを行い, work[0] にその値を返す. [out] rwork[] 配列 rwork[lrwork] (lrwork >= 8*n) 作業領域. [out] info = 0: 正常終了 = -1: 入力パラメータ jobvl の誤り (jobvl != 'V' および 'N') = -2: 入力パラメータ jobvr の誤り (jobvr != 'V' および 'N') = -3: 入力パラメータ n の誤り (n < 0) = -4: 入力パラメータ lda の誤り (lda < max(1, n)) = -6: 入力パラメータ ldb の誤り (ldb < max(1, n)) = -10: 入力パラメータ ldvl の誤り (ldvl が小さすぎる) = -12: 入力パラメータ ldvr の誤り (ldvr が小さすぎる) = -15: 入力パラメータ lwork の誤り (lwork が小さすぎる) = i (0 < i <= n): QZ反復が収束しなかった. 固有ベクトルは計算されないが, alpha[j], beta[j] (j = i〜n-1) は正しい値である. = n+1: zhgeqz においてQZ反復の収束以外のエラーが起きた. = n+2: ztgevc においてエラーがあった. 出典 LAPACK