|
|
◆ Drotm()
| Sub Drotm |
( |
N As |
Long, |
|
|
X_I As |
Double, |
|
|
Y_I As |
Double, |
|
|
P() As |
Double, |
|
|
Optional IncX As |
Long = 1, |
|
|
Optional IncY As |
Long = 1 |
|
) |
| |
修正ギブンス変換を適用する (BLAS 1)
- 目的
- 本ルーチンはDrotmgで求めた修正ギブンス変換を適用する.
- 引数
-
| [in] | N | ベクトル x および y の要素数 (N <= 0 の場合, 処理を行わずに戻る) |
| [in,out] | X_I | 配列X()の一要素X(I).
[in] ここを先頭にベクトルxを要素間隔IncXで格納する.
[out] ここを先頭に変換後のベクトル x が要素間隔IncXで格納される. 変換後のx(i) = H*x(i) + H*y(i) (i = 1 〜 N) |
| [in,out] | Y_I | 配列Y()の一要素Y(I).
[in] ここを先頭にベクトルyを要素間隔IncYで格納する.
[out] ここを先頭に変換後のベクトル y が要素間隔IncYで格納される. 変換後のy(i) = H*y(i) - H*x(i) (i = 1 〜 N) |
| [in] | P() | 配列 P(LP - 1) (LP >= 5)
変換行列 H.
P(0): -1, 0, 1 または -2.
P(1), P(2), P(3), P(4): H11, H21, H12, H22 の値. |
| [in] | IncX | (省略可)
配列X()内におけるベクトル x の要素間隔. (省略時 = 1) |
| [in] | IncY | (省略可)
配列Y()内におけるベクトル y の要素間隔. (省略時 = 1) |
- 出典
- BLAS
|