prod1d

    ICCG法形式で圧縮格納された対称行列[aij] とベクトル x の積 Ax を計算する倍精度サブルーチン.
     
     
  • サブルーチン名: prod1d

  •  
  • 呼び出し方:call prod1d(n, nzu, nzl, iu, il, au, al, ad, x, ap, ksu, ksl)

  •  
  • 引 数:
  • 引数
    意 味
    入入力の区分
    備 考
    n 行列 [aij] の寸法
    入力
    整数型
    nzu(n) 行列 [aij] の上三角部分の第i行目にある非零要素の数
    入力
    整数型
    nzl(n) 行列 [aij] の下三角部分の第i行目にある非零要素の数
    入力
    整数型
    iu(*) iu(k)は上三角部分にある第k番目の非零要素の列番号を表す.(k=1, 2, ..., n_nonzero)
    入力
    整数型
    il(*) il(k)は下三角部分にある第k番目の非零要素の列番号を表す.(k=1, 2, ..., n_nonzero)
    入力
    整数型
    au(*) au(k)は上三角部分にある第k番目の非零要素を表す.(k=1, 2, ..., n_nonzero)
    入力
    倍精度実数型
    al(*) al(k)は下三角部分にある第k番目の非零要素を表す.(k=1, 2, ..., n_nonzero)
    入力
    倍精度実数型
    x(n) ベクトル x の成分を格納する1次元配列.
    入力
    倍精度実数型
    ap(n) 行列・ベクトル積 Ax の成分を格納する1次元配列.
    出力
    倍精度実数型
    ksu(n) 行列 [aij] の上三角部分の第i行目よりも上にある非零要素の数
    入力
    整数型
    ksl(n) 行列 [aij] の下三角部分の第i行目よりも上にある非零要素の数
    入力
    整数型

  • 注 意:正方行列[aij] の情報の情報として,n, n_nonzero, nzu(n), iu(n_nonzero), au(n_nonzero)しか与えられていない場合は,本ルーチンを呼び出す前に,サブルーチンk_lower_from_iccgを用いて,nzl(n), il(n_nonzero), al(n_nonzero), ksu(n), ksl(n)を求めておく必要がある.