!SOLVER¶
ソルバーの制御
必須の制御データ．
パラメータ¶
METHOD = 解法 (CG、BiCGSTAB、GMRES、GPBiCG、DIRECT、DIRECTmkl、MUMPS)
DIRECT : 接触解析以外での直接法(逐次処理のみ) (現在使用不可)
DIRECTmkl : Intel MKLによる直接法
MUMPS : MUMPSによる直接法
直接法を選択したとき、データ行は無視される。
3自由度用の反復法はOpenMPによるスレッド並列が利用可能
PRECOND = 反復法の前処理手法 (1､2､3､5､10､11､12)
1, 2 : (Block) SSOR (3自由度用のみマルチカラーオーダリング付き)
3 : (Block) Diagonal Scaling
5 : マルチグリッド前処理パッケージMLによるAMG
10 : Block ILU(0)
11 : Block ILU(1)
12 : Block ILU(2)
10､11､12は3自由度問題でのみ利用可能
OpenMPによるスレッド並列時はSSOR, Diagonal ScalingまたはMLを推奨
ITERLOG = 反復法ソルバー収束履歴出力の有無 (YES/NO)(デフォルト: NO)
TIMELOG = ソルバー計算時間出力の有無 (YES/NO/VERBOSE)(デフォルト: NO）
USEJAD = ベクトル機向けオーダリングの有無 (YES/NO)(デフォルト: NO)
3自由度問題で反復法使用時のみ有効
SCALING = 行列の対角成分を1とするスケーリングの有無 (YES/NO)(デフォルト: NO)
3自由度問題で反復法使用時のみ有効
DUMPTYPE = 行列ダンプ型式(NONE、MM、CSR、BSR) (主にデバッグ用)
NONE : ダンプしない(デフォルト)
MM : マトリックスマーケット型式
CSR : Compressed Sparse Row（CSR）型式
BSR : Blocked CSR型式
DUMPEXIT = 行列ダンプ直後のプログラム終了 (YES/NO)(デフォルト: NO)
MPCMETHOD = 多点拘束条件の処理手法(1､2､3)
1: ペナルティ法（直接法使用時のデフォルト）
2: MPC-CG法 (非推奨)
3: 陽的自由度消去法（反復法使用時のデフォルト）
ESTCOND = 条件数推定の頻度 (試験的)
指定された反復ごと、および、反復終了時に条件数推定を実施
0の場合は推定を行わない
METHOD2 = 第2の解法 (BiCGSTAB、GMRES、GPBiCG) (試験的)
METHODにCGを指定した場合のみ有効
CG法が発散した場合に自動的に切り替えて求解を行う
他のパラメータやデータ行の情報は同じものが利用される
CONTACT_ELIM = 接触解析において自由度消去を行うかどうかの指定 (0,1)
0: 反復法使用時のみ自由度消去を行う(デフォルト)
1: 常に(直接法使用時にも)自由度消去を行う
** 2行目以降 **
| 変数名 | 属性 | 内容 |
|---|---|---|
| NITER | I | 反復回数(デフォルト: 100) |
| iterPREmax | I | Additive Schwarzによる前処理の繰り返し数(デフォルト: 1) (推奨値は1 (並列計算では2が有効な場合もある)) |
| NREST | I | クリロフ部分空間数(デフォルト: 10) (解法としてGMRESを選択したときのみ有効) |
| NCOLOR_IN | I | マルチカラーオーダリングにおける目標色数(デフォルト: 10) (OpenMPのスレッド数が2以上の時のみ有効) |
| RECYCLEPRE | I | 前処理セットアップ情報の再利用回数(デフォルト: 3) (非線形解析でのみ有効) |
| 変数名 | 属性 | 内容 |
|---|---|---|
| RESID | R | 打ち切り誤差(デフォルト値: 1.0e-8) |
| SIGMA_DIAG | R | 前処理行列計算時に対角成分にかける倍率(デフォルト値: 1.0) (ILU前処理でゼロ割や発散が起きる場合に1.0より大きい値を設定すると解ける場合がある) |
| SIGMA | R | 未使用(デフォルト値： 0.0) |
PRECOND=5の場合(省略可)¶
PRECOND にその他の値が指定された場合は、4行目は無視される。
| 変数名 | 属性 | 内容 |
|---|---|---|
| ML_CoarseSolver | I | MLの粗グリッド用ソルバ(1: スムーザーで代用、2: KLU(逐次直接法)、3: MUMPS(並列直接法))(デフォルト値: 1) (推奨値は、解き難い問題では3または2、通常の問題では1) |
| ML_Smoother | I | MLのスムーザー(1: Chebyshev、2: SymBlockGaussSeidel、3: Jacobi)(デフォルト値: 1) (推奨値は1) |
| ML_MGCycle | I | MLのマルチグリッドサイクル(1: V-cycle、2: W-cycle、3: Full-V-cycle)(デフォルト値: 1) (推奨値は、解き難い問題では2、通常の問題では1) |
| ML_MaxLevels | I | MLの最大レベル数(デフォルト値: 10) (推奨値は、特に解き難い問題では、コースソルバを直接法にして、2(メモリが足りない場合は3)、通常の問題では10) |
| ML_CoarseningScheme | I | MLのコースニングスキーム(1: UncoupledMIS、2: METIS、4: Zoltan、5: DD)(デフォルト値: 1) (推奨値は1または5) |
| ML_NumSweep | I | MLのスムーザーのスウィープ数(Chebyshevの場合は多項式の次数)(デフォルト値: 2) (推奨値は、Chebyshevの場合は2、SymBlockGaussSeidelの場合は1) |
使用例¶
SSOR前処理付きCG法を利用し、最大反復階数を10000、打ち切り誤差を1.0e-8に設定する
SSOR前処理付きGMRES法を利用し、クリロフ部分空間数を40、マルチカラーオーダリングの目標色数を100に設定する
ILU(0)前処理付きCG法を利用し、前処理行列計算時に対角成分にかける倍率を1.1に設定する
CG法の前処理をMLによるAMG法とする
CG法の前処理をMLによるAMG法とし、粗グリッド用ソルバをMUMPSとする(解き難い問題向け)
CG法の前処理をMLによるAMG法とし、マルチグリッドサイクルをW-cycleとする(解き難い問題向け)
CG法の前処理をMLによるAMG法とし、粗グリッド用ソルバをMUMPS、最大レベル数を2とする(非常に解き難い問題向け)