◆ dorgtr()

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

目的: 本ルーチンは dsytrd が返す n - 1 個の基本鏡映 H(i) の積により定義される 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': a[][] の上三角部分に dsytrd から基本鏡映を格納する. = 'L': a[][] の下三角部分に dsytrd から基本鏡映を格納する.
[in]	n	行列 Q の行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in]	lda	二次元配列 a[][] の整合寸法. (lda >= max(1, n))
[in,out]	a[][]	配列 a[la][lda] (la >= n) [in] dsytrd が返す基本鏡映を定義するベクトル. [out] n x n 直交行列 Q.
[in]	tau[]	配列 tau[ltau] (ltau >= n - 1) tau[i] は dsytrd により返された基本鏡映 H(i) のスカラー係数でなければならない.
[out]	work[]	配列 work[lwork] 作業領域. info = 0 の場合, work[0] に lwork の最適値を返す.
[in]	lwork	配列 work[]のサイズ (lwork >= max(1, n-1)) 最適パフォーマンスのためには lwork >= (n-1)*nb とせよ. ただし, nb は最適ブロックサイズである. lwork = -1 の場合, 作業領域サイズの問い合わせとみなし, work[] の最適サイズを求める計算だけを行い, work[0] にその値を返す.
[out]	info	= 0: 正常終了. = -1: 入力パラメータ uplo の誤り. (uplo != 'U' および 'L') = -2: 入力パラメータ n の誤り. (n < 0) = -3: 入力パラメータ lda の誤り. (lda < max(1, n)) = -7: 入力パラメータ lwork の誤り. (lwork が小さすぎる)