k_qr_decomp2

    非正方行列 [aij] を QR 分解する倍精度サブルーチン．
     
     
  • サブルーチン名： k_qr_decomp2

    •  
  • 呼び出し方：call k_qr_decomp2(mdim, ndim, adim, awork, tau, jpvt)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    mdim, ndim 行列 [aij] の寸法．mdim ≥ ndim を満たす必要がある．
    入力
    		 整数型
    adim(mdim, ndim) 行列 [aij] を格納する2次元配列
    入力
    		 倍精度実数型
    awork(mdim, ndim) QR 分解情報を格納する2次元配列
    出力
    		 倍精度実数型
    tau(ndim) QR 分解情報を格納する1次元配列
    出力
    		 倍精度実数型
    jpvt(ndim) QR 分解を行った際の行列 [aij] の列交換に関する情報．
    第 k 列は第jpvt(k)列と交換される．
    出力
    		 整数型

  • 注　意：サブルーチン k_matvec_px2, k_matvec_qtx2, k_matvec_qx2, k_matvec_rmtx2, k_matvec_ufx2 を用いる際には，前もって本ルーチンを呼び出しておく必要がある．