k_qr_decomp

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

    •  
  • 呼び出し方：call k_qr_decomp(mdim, ndim, adim)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    mdim, ndim 行列 [aij] の寸法．mdim ≥ ndim を満たす必要がある．
    入力
    		 整数型
    adim(mdim, ndim) 行列 [aij] を格納する2次元配列
    入力
    		 倍精度実数型

  • 注　意：
    1. 本ルーチンを実行すると，QR分解情報が module名 mdl_k_qr_decomp 中に格納される．
    2. サブルーチン k_matvec_qtx, k_matvec_qx, k_matvec_ufx を用いる際には，前もって本ルーチンを呼び出しておく必要がある．