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