|
|
◆ zunghr()
| void zunghr |
( |
int |
n, |
|
|
int |
ilo, |
|
|
int |
ihi, |
|
|
int |
lda, |
|
|
doublecomplex |
a[], |
|
|
doublecomplex |
tau[], |
|
|
doublecomplex |
work[], |
|
|
int |
lwork, |
|
|
int * |
info |
|
) |
| |
ヘッセンベルグ形への変換行列の生成 (複素行列)
- 目的
- 本ルーチンは zgehrd により返される ihi-ilo 個の n 次基本鏡映の積として定義されるユニタリ行列 Q を生成する.
Q = H(ilo) H(ilo+1) . . . H(ihi-1).
- 引数
-
| [in] | n | 行列 A の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) |
| [in] | ilo | |
| [in] | ihi | ilo と ihi は zgehrd を呼び出したときと同じ値でなければならない. Q は小行列 Q(ilo+1〜ihi, ilo+1〜ihi) を除いては単位行列に等しい. (1 <= ilo <= ihi <= n (n > 0 の場合). ilo = 1 かつ ihi = 0 (n = 0 の場合)) |
| [in] | lda | 二次元配列 a[][] の整合寸法. (lda >= max(1, n)) |
| [in,out] | a[][] | 配列 a[la][lda] (la >= n)
[in] zgehrd により返された基本鏡映を定義するベクトル.
[out] n x n ユニタリ行列 Q. |
| [in] | tau[] | 配列 tau[ltau] (ltau >= n - 1) tau[i] は zgehrd により返された基本鏡映 H(i) のスカラー係数でなければならない. |
| [out] | work[] | 配列 work[lwork]
作業領域.
info = 0 の場合, work[0] に lwork の最適値を返す. |
| [in] | lwork | 配列work[]のサイズ. (lwork >= ihi-ilo)
最適パフォーマンスのためには lwork >= (ihi-ilo)*nb とせよ. ただし, nb は最適ブロックサイズである.
lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, 最適サイズを求める計算だけを行い, work[0] にその値を返す. |
| [out] | info | = 0: 正常終了
= -1: 入力パラメータ n の誤り (n < 0)
= -2: 入力パラメータ ilo の誤り (ilo < 1 または ilo > n)
= -3: 入力パラメータ ihi の誤り (ihi < min(ilo, n) または ihi > n)
= -4: 入力パラメータ lda の誤り (lda < max(1, n))
= -8: 入力パラメータ lwork の誤り (lwork が小さすぎる)
|
- 出典
- LAPACK
|