k_product_ax_ccs

    CCS法を用いて格納された非正方行列 \(A=(a_{ij})\) とベクトル \(\boldsymbol{x}=(x_i)\) との積を計算する倍精度サブルーチン．
     
     
  • サブルーチン名： k_product_ax_ccs

    •  
  • 呼び出し方：call k_product_ax_ccs(mdim, ndim, n_nonzero, a_val, row_ind, col_ptr, x, ax)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    mdim, ndim 行列 \(A\) の寸法
    入力
    		 整数型
    n_nonzero 行列 \(A\) の非零要素数
    入力
    		 整数型
    a_val(n_nonzero) 非零要素を格納する1次元配列
    入力
    		 倍精度実数型
    row_ind(n_nonzero) 非零要素の行番号を格納する1次元配列
    入力
    		 整数型
    col_ptr(ndim+1) col_ptr(j)はj-列目上にある最初の非零要素の非零要素番号を示す．
    入力
    		 整数型
    x(ndim) ベクトル \(\boldsymbol{x}\) を格納する1次元配列
    入力
    		 倍精度実数型
    ax(mdim) 行列ベクトル積 \(A\boldsymbol{x}\) を格納する1次元配列
    出力
    		 倍精度実数型

  • 注　意：既に非正方行列 \(A\) が定まっている場合には，サブルーチン k_convm_2_ccs を用いて，n_nonzero, a_val(n_nonzero), row_ind(n_nonzero), col_ptr(ndim+1)を計算することができる．