k_matvec_ufx2

    非正方行列 [aij] の QR 分解で得られる行列 F ≡ QQT，U ≡ E - Uとベクトル x の積を計算する倍精度サブルーチン．
     
     
  • サブルーチン名： k_matvec_ufx2

    •  
  • 呼び出し方：call k_matvec_ufx2(ndim, kdim, awork, tau, xdim, fxdim, uxdim)

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

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