◆ WZhbev2()

固有値・固有ベクトル (エルミート帯行列) (実数部／虚数部分離形式)

目的: WZhbev2はエルミート帯行列 A のすべての固有値, および, 必要により固有ベクトルを求める.

複素数を表現するために実数部と虚数部を隣り合ったセルに格納する(左が実数部, 右が虚数部). 得られた解も実数部と虚数部が隣り合った別々のセルに出力される.

戻り値

N+1×1 (Jobz = "N"の場合), N+1×2N+1 (Jobz = "V" かつ Info = 0の場合)

	列1	列2〜2N+1
行1〜N	固有値 (昇順)	固有ベクトル(Jobz = "V" かつ Info = 0の場合) (実数部と虚数部を隣り合った列に格納(左が実数部, 右が虚数部))
行N+1	リターンコード	0

リターンコード.
= 0: 正常終了.
= i > 0: 中間結果の3重対角形の非対角要素のうちi個が0に収束しなかった.

引数

[in]	Jobz	= "N": 固有値のみ求める. = "V": 固有値と固有ベクトルを求める.
[in]	Uplo	= "U": Aの上三角部分を格納する. = "L": Aの下三角部分を格納する.
[in]	N	行列Aの行および列数. (N >= 1)
[in]	Kd	上または下帯幅. (対角要素を含まない) (Kd >= 0)
[in]	A	(Kd+1 × 2N) N×N対称帯行列 A. (対称帯行列形式. 詳細は下記参照)

詳細

次の例は, 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   *    *

*で示された配列要素は使用されない.

使用例: エルミート帯行列Aの固有値・固有ベクトルを求める. ただし,
( 2.20 -0.32-0.81i 0 )

A = ( -0.32+0.81i 2.11 0.37+0.80i )

( 0 0.37-0.80i 2.93 )

とする.

WZhbev2