|
|
◆ csr_ilu_solve()
| void csr_ilu_solve |
( |
char |
trans, |
|
|
int |
n, |
|
|
const double |
val[], |
|
|
const int |
rowptr[], |
|
|
const int |
colind[], |
|
|
int |
base, |
|
|
const double |
d[], |
|
|
const double |
b[], |
|
|
double |
x[], |
|
|
int * |
info |
|
) |
| |
不完全LU分解による前処理 (ILU) (CSR)
- 目的
- 係数行列の不完全LU分解による前処理(ILU)を行う. すなわち, 連立一次方程式 M*x = b または M^T*x = b を解く. ここで, M は前処理行列である.
- 引数
-
| [in] | trans | = 'N': M*x = b を解く.
= 'T' または 'C': M^T*x = b を解く. |
| [in] | n | 前処理行列の次数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る) |
| [in] | val[] | 配列 val[lval] (lval >= nnz) (nnz は非ゼロ要素数)
前処理行列(下三角行列 L と上三角行列 U)の非ゼロ要素の値. |
| [in] | rowptr[] | 配列 rowptr[lrowptr] (lrowptr >= n + 1)
前処理行列の行ポインタ. |
| [in] | colind[] | 配列 colind[lcolind] (lcolind >= nnz)
前処理行列の列インデクス. |
| [in] | base | rowptr[] および colind[] のインデクス形式.
= 0: 0-ベース(C形式): 開始インデクス値が 0.
= 1: 1-ベース(Fortran形式): 開始インデクス値が 1. |
| [in] | d[] | 配列 d[ld] (ld >= n)
上三角行列 U の対角要素. |
| [in] | b[] | 配列 b[lb] (lb >= n)
右辺ベクトル b. |
| [out] | x[] | 配列 x[lx] (lx >= n)
解ベクトル x. |
| [out] | info | = 0: 正常終了.
= i < 0: (-i)番目の入力パラメータの誤り. |
|