k_mod_val

    CRS/CCS法を用いて格納された非正方行列\(A=(a_{ij})\)の情報中に非零要素 \(a_{ij}\)の情報を追加するために，非零要素の値を格納する配列 val に空き要素を追加する．
    推奨される応用例：省メモリ版FEFG，X−EFG用．
     
     
  • サブルーチン名： k_mod_val

    •  
  • 呼び出し方：call k_mod_val(nzdim, val, n_zero, k_zero, icode)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    nzdim 非零要素数の見積値．配列col_ind，valの整合寸法．
    入力
    		 整数型
    val(nzdim) 非零要素の値を格納する1次元配列．
    入出力
    		 倍精度実数型
    n_zero 行列 \(A\) の非零要素数
    入力
    		 整数型
    k_zero 非零要素\(a_{ij}\)の非零要素番号
    入力
    		 整数型
    icode icode=1の場合には，配列valのk_zero番目の要素以降は右側にシフトし，icode=0の場合は，配列valには操作を行わない．
    入力
    		 整数型

  • 注　意：
    1. n_zero, k_zero, icodeを決定するために，本ルーチンを最初に呼び出す前に k_add_nz2_crs または k_add_nz2_ccs を呼び出す必要がある．
    2. 使用例（X-EFG法）
      3.         do n=1,ndom
           i=idom(n)
           do m=1,ndom
              j=idom(m)
              call k_add_nz2_crs(i,j,nnode,n_limit,col_ind,row_ptr,n_zero,k_zero,icode)
              call k_mod_val(n_limit,a_val,n_zero,k_zero,icode)
              a_val(k_zero) = a_val(k_zero)+(psix(i)*phayx(j)+psiy(i)*phayy(j)) &
                                    *weit(ix)*weit(iy)*djaco(nc)
           end do
        end do
    3. X-EFG法やEFG法では，バックグラウンド・セルを用いた数値積分を実行している．それ故，本ルーチンの実行時間は殆ど無視できる． これに対して，行列要素が解析的に計算できるFEMでは，本ルーチンの実行時間は他の計算時間と比べて無視できない．このような場合，本ルーチンは使うべきではない．