◆ _dstev()

(シンプルドライバ) 固有値・固有ベクトル (対称3重対角行列)

目的: 本ルーチンは実対称3重対角行列 A のすべての固有値, および, 必要により固有ベクトルを求める. 固有ベクトルを求める場合には QL法または QR法を使用する.

引数

[in]	jobz	= 'N': 固有値のみ求める. = 'V': 固有値と固有ベクトルを求める.
[in]	n	行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in,out]	d[]	配列 d[ld] (ld >= n) [in] 対称3重対角行列Aのn個の対角要素. [out] info = の場合, 求められた固有値(昇順).
[in,out]	e[]	配列 e[le] (le >= n - 1) [in] 対称3重対角行列Aのn-1個の副対角要素(e[0]〜e[n-2]に格納). [out] e[]の内容は壊される.
[in]	ldz	配列z[][]の整合寸法. (ldz >= 1 (jobz='N'), ldz >= max(1,n) (jobz='V'))
[out]	z[][]	配列 z[lz][ldz] (lz >= n) jobz = 'V': info = の場合, z[][]に正規直交固有ベクトルが入る. w[i]に関連する固有ベクトルがi行に入る. jobz = 'N': z[][]は参照されない.
[out]	work[]	配列 work[lwork] (lwork >= max(1, 2*n-2)) 作業領域. jobz = 'N'の場合, work[]は参照されない.
[out]	info	= 0: 正常終了 = -1: 入力パラメータ jobz の誤り (jobz != 'V'および'N') = -2: 入力パラメータ n の誤り (n < 0) = -5: 入力パラメータ ldz の誤り (ldzが小さすぎる) = i > 0: 収束しなかった. e[]の非対角要素のうちi個が0にならなかった.