k_matvec_rtx2
非正方行列 [
a
ij
] の
QR
分解で得られた
R
T
とベクトル
x
の積を計算する倍精度サブルーチン.
サブルーチン名: k_matvec_rtx2
呼び出し方:call k_matvec_rtx2(mdim, ndim, awork, xdim, ydim)
引 数:
引数
意 味
入出力の区分
備 考
mdim, ndim
行列 [
a
ij
] の寸法.mdim ≥ ndim を満たす必要がある.
入力
整数型
awork(mdim, ndim)
QR
分解情報を格納する2次元配列
入力
倍精度実数型
xdim(ndim)
ベクトル
x
の成分を格納する1次元配列.
入力
倍精度実数型
ydim(ndim)
行列・ベクトル積
R
T
x
の成分を格納する1次元配列.
出力
倍精度実数型
注 意:本ルーチンを呼び出す前に,サブルーチン k_qr_decomp2 を呼び出しておく必要がある.