◆ dopgtr()

三重対角形への変換行列の生成 (対称行列) (圧縮形式)

目的: 本ルーチンは dsptrd が圧縮形式で返す n - 1 個の基本鏡映の積により定義される n 次実直交行列 Q を生成する.
uplo = 'U' の場合, Q = H(n-1) . . . H(2) H(1).

uplo = 'L' の場合, Q = H(1) H(2) . . . H(n-1).

引数

[in]	uplo	= 'U': 圧縮形式の上三角部分が dsptrd の前回呼び出しで使われた. = 'L': 圧縮形式の下三角部分が dsptrd の前回呼び出しで使われた.
[in]	n	行列 Q の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in]	ap[]	配列 ap[lap] (lap >= n(n + 1)/2) dsptrd が返す基本鏡映を定義するベクトル.
[in]	tau[]	配列 tau[ltau] (ltau >= n - 1) tau[i] は dsptrd により返された基本鏡映 H(i) のスカラー係数でなければならない.
[in]	ldq	二次元配列 q[][] の整合寸法. (ldq >= max(1, n))
[out]	q[][]	配列 q[lq][ldq] (lq >= n) n x n 直交行列 Q.
[out]	work[]	配列 work[lwork] (lwork >= n - 1) 作業領域.
[out]	info	= 0: 正常終了. = -1: 入力パラメータ uplo の誤り. (uplo != 'U' および 'L') = -2: 入力パラメータ n の誤り. (n < 0) = -5: 入力パラメータ ldq の誤り. (ldq < max(1, n))