k_matvec_qtx2

    非正方行列 [aij] の QR 分解で得られたQTとベクトル x の積を計算する倍精度サブルーチン．
     
     
  • サブルーチン名： k_matvec_qtx2

    •  
  • 呼び出し方：call k_matvec_qtx2(mdim, ndim, awork, tau, xdim, qxdim)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    mdim, ndim 行列 [aij] の寸法．mdim ≥ ndim を満たす必要がある．
    入力
    		 整数型
    awork(mdim, ndim) QR 分解情報を格納する2次元配列
    入力
    		 倍精度実数型
    tau(ndim) QR 分解情報を格納する1次元配列
    入力
    		 倍精度実数型
    xdim(mdim) ベクトル x の成分を格納する1次元配列．
    入力
    		 倍精度実数型
    qxdim(ndim) 行列・ベクトル積 QTx の成分を格納する1次元配列．
    出力
    		 倍精度実数型

  • 注　意：本ルーチンを呼び出す前に，サブルーチン k_qr_decomp2 を呼び出しておく必要がある．