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