k_matvec_ptx2
非正方行列 [
a
ij
] の
QR
分解で得られた置換行列
P
の転置行列
P
T
とベクトル
x
の積を計算する倍精度サブルーチン．
サブルーチン名： k_matvec_ptx2
呼び出し方：call k_matvec_ptx2(ndim, jpvt, xdim, pxdim)
引 数：
引数
意 味
入出力の区分
備 考
ndim
行列 [
a
ij
] の列数
入力
整数型
jpvt(ndim)
QR
分解を行った際の行列 [
a
ij
] の列交換に関する情報．
第
k
列は第jpvt(
k
)列と交換される．
入力
整数型
xdim(ndim)
ベクトル
x
の成分を格納する1次元配列
入力
倍精度実数型
pxdim(ndim)
P
T
x
の成分を格納する1次元配列
出力
倍精度実数型
注 意：本ルーチンを呼び出す際には，前もってサブルーチン k_qr_decomp2 を呼び出しておく必要がある．