You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
69 lines
4.5 KiB
Plaintext
69 lines
4.5 KiB
Plaintext
%%%%%%%%%%% General arguments % Lines starting with % are ignored.
|
|
CSR ! Storage format CSR COO JAD
|
|
0200 ! IDIM; domain size. Linear system size is IDIM**2
|
|
CONST ! PDECOEFF: CONST, EXP, BOX Coefficients of the PDE
|
|
CG ! Iterative method: BiCGSTAB BiCGSTABL BiCG CG CGS FCG GCR RGMRES
|
|
2 ! ISTOPC
|
|
00500 ! ITMAX
|
|
1 ! ITRACE
|
|
30 ! IRST (restart for RGMRES and BiCGSTABL)
|
|
1.d-6 ! EPS
|
|
%%%%%%%%%%% Main preconditioner choices %%%%%%%%%%%%%%%%
|
|
ML-VCYCLE-BJAC-D-BJAC ! Longer descriptive name for preconditioner (up to 20 chars)
|
|
ML ! Preconditioner type: NONE JACOBI GS FBGS BJAC AS ML
|
|
%%%%%%%%%%% First smoother (for all levels but coarsest) %%%%%%%%%%%%%%%%
|
|
BJAC ! Smoother type JACOBI FBGS GS BWGS BJAC AS. For 1-level, repeats previous.
|
|
1 ! Number of sweeps for smoother
|
|
0 ! Number of overlap layers for AS preconditioner
|
|
HALO ! AS restriction operator: NONE HALO
|
|
NONE ! AS prolongation operator: NONE SUM AVG
|
|
ILU ! Subdomain solver for BJAC/AS: JACOBI GS BGS ILU ILUT MILU MUMPS SLU UMF
|
|
LLK ! AINV variant, ignored otherwise
|
|
0 ! Fill level P for ILU(P) and ILU(T,P)
|
|
1 ! Inverse Fill level P for INVK
|
|
1.d-4 ! Threshold T for ILU(T,P)
|
|
%%%%%%%%%%% Second smoother, always ignored for non-ML %%%%%%%%%%%%%%%%
|
|
NONE ! Second (post) smoother, ignored if NONE
|
|
1 ! Number of sweeps for (post) smoother
|
|
0 ! Number of overlap layers for AS preconditioner
|
|
HALO ! AS restriction operator: NONE HALO
|
|
NONE ! AS prolongation operator: NONE SUM AVG
|
|
ILU ! Subdomain solver for BJAC/AS: JACOBI GS BGS ILU ILUT MILU MUMPS SLU UMF
|
|
LLK ! AINV variant, ignored otherwise
|
|
0 ! Fill level P for ILU(P) and ILU(T,P)
|
|
8 ! Inverse Fill level P for INVK
|
|
1.d-4 ! Threshold T for ILU(T,P)
|
|
%%%%%%%%%%% Multilevel parameters %%%%%%%%%%%%%%%%
|
|
VCYCLE ! Type of multilevel CYCLE: VCYCLE WCYCLE KCYCLE MULT ADD
|
|
4 ! Number of outer sweeps for ML
|
|
-3 ! Max Number of levels in a multilevel preconditioner; if <0, lib default
|
|
-3 ! Target coarse matrix size; if <0, lib default
|
|
SMOOTHED ! Type of aggregation: SMOOTHED UNSMOOTHED
|
|
DEC ! Parallel aggregation: DEC, SYMDEC
|
|
NATURAL ! Ordering of aggregation NATURAL DEGREE
|
|
FILTER ! Filtering of matrix: FILTER NOFILTER
|
|
-1.5 ! Coarsening ratio, if < 0 use library default
|
|
-2 ! Number of thresholds in vector, next line ignored if <= 0
|
|
0.05 0.025 ! Thresholds
|
|
-0.0100d0 ! Smoothed aggregation threshold, ignored if < 0
|
|
%%%%%%%%%%% Coarse level solver %%%%%%%%%%%%%%%%
|
|
RKR ! Coarsest-level solver: MUMPS(global) UMF SLU SLUDIST JACOBI GS BJAC RKR(global)
|
|
NONE ! Coarsest-level subsolver for BJAC: ILU ILUT MILU UMF MUMPS(local) SLU RKR(local)
|
|
LLK ! AINV Variant for the coarse solver
|
|
BICGSTAB ! Krylov method for RKR solver/subsolver, ignored otherwise
|
|
DIST ! Coarsest-level matrix distribution: DIST REPL
|
|
0 ! Coarsest-level fillin P for ILU(P) and ILU(T,P)
|
|
1 ! Coarsest-level inverse Fill level P for INVK
|
|
1.d-4 ! Coarsest-level threshold T for ILU(T,P)
|
|
30 ! Number of sweeps for JACOBI/GS/BJAC coarsest-level solver
|
|
30 ! maxit for RKR
|
|
1.d-4 ! eps for RKR
|
|
30 ! itrace for RKR
|
|
T ! Check the BJAC residual
|
|
T ! Print the BJAC residual
|
|
5 ! ITRACE for residual check
|
|
5 ! ITRACE for residual print
|
|
1.d-4 ! Tolerance for exit from BJAC
|
|
% dump
|
|
T ! Dump preconditioner
|