dgees dgees_r dgeesx dgeesx_r dgeev dgeevx

dgeev() void dgeev ( char  jobvl, char  jobvr, int  n, int  lda, double  a[], double  wr[], double  wi[], int  ldvl, double  vl[], int  ldvr, double  vr[], double  work[], int  lwork, int *  info  ) (シンプルドライバ) 固有値・固有ベクトル (一般行列) 目的 本ルーチンは n x n 一般実行列 A のすべての固有値, および, 必要により左および／または右固有ベクトルを求める. はじめに行列 A を上ヘッセンベルグ形に変換し, 次にQR法により固有値, および, 必要によりシュール分解を求める. 固有ベクトルはシュール分解より計算される. A の右固有ベクトル v(j) は次式を満たす. A * v(j) = lambda(j) * v(j) ここで, lambda(j) は固有値である. A の左固有ベクトル u(j) は次式を満たす. u(j)^H * A = lambda(j) * u(j)^H ここで, u(j)^H は u(j) の共役転置を表す。 求められた固有ベクトルはユークリッドノルムが 1 で最大成分が実数になるように正規化される. 引数 [in] jobvl = 'N': A の左固有ベクトルを求めない. = 'V': A の左固有ベクトルを求める. [in] jobvr = 'N': A の右固有ベクトルを求めない. = 'V': A の右固有ベクトルを求める. [in] n 行列 A の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) [in] lda 二次元配列 a[][] の整合寸法. (lda >= max(1, n)) [in,out] a[][] 配列 a[la][lda] (la >= n) [in] n x n 行列 A. [out] a[][] は上書きされる. [out] wr[] 配列 wr[lwr] (lwr >= n) [out] wi[] 配列 wi[lwi] (lwi >= n) wr[] および wi[] にはそれぞれ求められた固有値の実数部および虚数部が入る. 複素共役対の固有値は, 虚数部が正の固有値を先頭に, 続けて現れる. [in] ldvl 二次元配列 vl[][] の整合寸法. (ldvl >= 1 (jobvl = 'N'), ldvl >= n (jobvl = 'V')) [out] vl[][] 配列 vl[lvl][ldvl] (lvl >= n) jobvl = 'V': 左固有ベクトル u(j) が, 対応する固有値と同順に vl[][] の行に順に格納される. j 番目の固有値が実数の場合, u(j) = (vl[][] の j 番目の行). j 番目と j+1 番目の固有値が複素共役対の場合, u(j) = (vl[][] の j 番目の行) + i*(vl[][] の j+1 番目の行), u(j+1) = (vl[][] の j 番目の行) - i*(vl[][] の j+1 番目の行). jobvl = 'N': vl[][] は参照されない. [in] ldvr 二次元配列 vr[][] の整合寸法. (ldvr >= 1 (jobvr = 'N'), ldvr >= n (jobvr = 'V')) [out] vr[][] 配列 vr[lvr][ldvr] (lvr >= n) jobvr = 'V': 右固有ベクトルv(j)が, 対応する固有値と同順にvr[][]の行に順に格納される. j 番目の固有値が実数の場合, v(j) = (vr[][] の j 番目の行). j 番目と j+1 番目の固有値が複素共役対の場合, v(j) = (vr[][] の j 番目の行) + i*(vr[][] の j+1 番目の行), v(j+1) = (vr[][] の j 番目の行) - i*(vr[][] の j+1 番目の行). jobvr = 'N': vr[][]は参照されない. [out] work[] 配列 work[lwork] 作業領域. info = 0 の場合, work[0] に lwork の最適値を返す. [in] lwork 配列work[]のサイズ. (lwork >= max(1, 3*n) (jobvl = 'N' かつ jobvr = 'N'), lwork >= max(1, 4*n) (jobvl = 'V' または jobvr = 'V')) 一般に, パフォーマンスをよくするためには lwork を大きくしなければならない. lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, 最適サイズを求める計算だけを行い, work[0] にその値を返す. [out] info = 0: 正常終了 = -1: 入力パラメータ jobvl の誤り (jobvl != 'V' および 'N') = -2: 入力パラメータ jobvr の誤り (jobvr != 'V' および 'N') = -3: 入力パラメータ n の誤り (n < 0) = -4: 入力パラメータ lda の誤り (lda < max(1, n)) = -8: 入力パラメータ ldvl の誤り (ldvl が小さすぎる) = -10: 入力パラメータ ldvr の誤り (ldvr が小さすぎる) = -13: 入力パラメータ lwork の誤り (lwork が小さすぎる) = i > 0: QRアルゴリズムが失敗しすべての固有値を求めることはできなかった. 固有ベクトルは計算されなかった. wr[] および wi[] の要素 i〜n-1 には収束した固有値が入る. 出典 LAPACK