h_matrix

    実対称密行列を節点情報によって階層行列(\(\mathcal{H}\)-matrix)に変換する倍精度サブルーチン.但し,元の実対称密行列は2次元平面内の節点情報に基づいて構成されている.
     
     
  • サブルーチン名: h_matrix

  •  
  • 呼び出し方:call h_matrix(xnode,nnode,nbmax,nmin,etha,epsln,wdim,ne_aca)

  •  
  • 引 数:
  • 引数
    意 味
    入出力の区分
    備 考
    xnode(nnode,2) 節点の \((x,y)\) 座標を格納する配列.第 \(i\) 番目節点の \(x,y\) 座標は,それぞれxnode(i, 1), xnode(i, 2) に格納される.
    入力
    nnode 全節点数
    入力
    nbmax ブロック行列の生成数の上限
    入力
    nmin クラスター木の葉に含まれる節点の最小数
    入力
    etha 許容条件に現れるパラメタ \(\eta\)
    入力
    \(1.0\) 程度
    epsln ブロック行列にACA近似を行う際の収束判定子 \(\epsilon\)
    入力
    \(10^{-5}\) 程度
    wdim(nnode,nnode) 階層行列(\(\mathcal{H}\)-matrix)に変換する行列を格納する2次元配列
    入力
    ne_aca 低ランク・ブロックに含まれる全要素数
    出力
  • 注 意
    H-行列法を用いて行列・ベクトル積を高速に計算するためには,サブルーチンproduct_hvを呼び出せばよい.その際には,前もって1回だけ本ルーチンを呼び出せばよい.