fdurastante
/
amg4psblas
mirror of https://github.com/sfilippone/amg4psblas.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

mld2p4-2:

Browse Source 
mlprec/mld_c_prec_type.f90
 mlprec/mld_d_prec_type.f90
 mlprec/mld_s_prec_type.f90
 mlprec/mld_z_prec_type.f90
 tests/fileread/runs/dfs.inp

Start work on aggregation control.
stopcriterion
Salvatore Filippone 10 years ago
parent d2e6fcb5b8
commit 46de48ffa9
5 changed files with 24 additions and 16 deletions
Show Stats Download Patch File Download Diff File

4
mlprec/mld_c_prec_type.f90
Unescape Escape View File

@ -81,7 +81,9 @@ module mld_c_prec_type
type, extends(psb_cprec_type) :: mld_cprec_type
integer(psb_ipk_) :: ictxt
integer(psb_ipk_) :: coarse_aggr_size
integer(psb_ipk_) :: coarse_aggr_size = izero
integer(psb_ipk_) :: n_prec_levs = -ione
integer(psb_ipk_) :: max_prec_levs = 20_psb_ipk_
real(psb_spk_) :: op_complexity=szero
type(mld_c_onelev_type), allocatable :: precv(:)
contains

4
mlprec/mld_d_prec_type.f90
Unescape Escape View File

@ -81,7 +81,9 @@ module mld_d_prec_type
type, extends(psb_dprec_type) :: mld_dprec_type
integer(psb_ipk_) :: ictxt
integer(psb_ipk_) :: coarse_aggr_size
integer(psb_ipk_) :: coarse_aggr_size = izero
integer(psb_ipk_) :: n_prec_levs = -ione
integer(psb_ipk_) :: max_prec_levs = 20_psb_ipk_
real(psb_dpk_) :: op_complexity=dzero
type(mld_d_onelev_type), allocatable :: precv(:)
contains

4
mlprec/mld_s_prec_type.f90
Unescape Escape View File

@ -81,7 +81,9 @@ module mld_s_prec_type
type, extends(psb_sprec_type) :: mld_sprec_type
integer(psb_ipk_) :: ictxt
integer(psb_ipk_) :: coarse_aggr_size
integer(psb_ipk_) :: coarse_aggr_size = izero
integer(psb_ipk_) :: n_prec_levs = -ione
integer(psb_ipk_) :: max_prec_levs = 20_psb_ipk_
real(psb_spk_) :: op_complexity=szero
type(mld_s_onelev_type), allocatable :: precv(:)
contains

4
mlprec/mld_z_prec_type.f90
Unescape Escape View File

@ -81,7 +81,9 @@ module mld_z_prec_type
type, extends(psb_zprec_type) :: mld_zprec_type
integer(psb_ipk_) :: ictxt
integer(psb_ipk_) :: coarse_aggr_size
integer(psb_ipk_) :: coarse_aggr_size = izero
integer(psb_ipk_) :: n_prec_levs = -ione
integer(psb_ipk_) :: max_prec_levs = 20_psb_ipk_
real(psb_dpk_) :: op_complexity=dzero
type(mld_z_onelev_type), allocatable :: precv(:)
contains

24
tests/fileread/runs/dfs.inp
Unescape Escape View File

@ -1,12 +1,12 @@
pressmat.mtx ! This matrix (and others) from: http://math.nist.gov/MatrixMarket/ or
pressrhs.mtx ! rhs | http://www.cise.ufl.edu/research/sparse/matrices/index.html
MM !
BICGSTAB ! Iterative method: BiCGSTAB BiCG CGS RGMRES BiCGSTABL CG
RGMRES ! Iterative method: BiCGSTAB BiCG CGS RGMRES BiCGSTABL CG
CSR ! Storage format: CSR COO JAD
0 ! IPART (partition method): 0 (block) 2 (graph, with Metis)
2 ! IPART (partition method): 0 (block) 2 (graph, with Metis)
2 ! ISTOPC
00100 ! ITMAX
01 ! ITRACE
00120 ! ITMAX
30 ! ITRACE
30 ! IRST (restart for RGMRES and BiCGSTABL)
1.d-6 ! EPS
4L-M-RAS-I-UR ! Longer descriptive name for preconditioner (up to 20 chars)
@ -14,21 +14,21 @@ ML ! Preconditioner type: NONE JACOBI BJAC AS ML
0 ! Number of overlap layers for AS preconditioner
HALO ! AS restriction operator: NONE HALO
NONE ! AS prolongation operator: NONE SUM AVG
GS ! AS subdomain solver: DSCALE ILU MILU ILUT UMF SLU MUMPS
ILU ! AS subdomain solver: DSCALE ILU MILU ILUT UMF SLU MUMPS
0 ! Fill level P for ILU(P) and ILU(T,P)
1.d-4 ! Threshold T for ILU(T,P)
4 ! Number of Jacobi sweeps for base smoother
2 ! Number of Jacobi sweeps for base smoother
4 ! Number of levels in a multilevel preconditioner
BJAC ! Smoother type JACOBI BJAC AS ignored for non-ML
NONSMOOTHED ! Type of aggregation: SMOOTHED NONSMOOTHED
SMOOTHED ! Type of aggregation: SMOOTHED NONSMOOTHED
DEC ! Type of aggregation: DEC
NATURAL ! Ordering of aggregation NATURAL DEGREE
DEGREE ! Ordering of aggregation NATURAL DEGREE
MULT ! Type of multilevel correction: ADD MULT
POST ! Side of multiplicative correction PRE POST TWOSIDE (ignored for ADD)
REPL ! Coarsest-level matrix distribution: DIST REPL
TWOSIDE ! Side of multiplicative correction PRE POST TWOSIDE (ignored for ADD)
DIST ! Coarsest-level matrix distribution: DIST REPL
BJAC ! Coarsest-level solver: JACOBI BJAC UMF SLU SLUDIST MUMPS
ILU ! Coarsest-level subsolver: ILU UMF SLU MUMPS SLUDIST
0 ! Coarsest-level fillin P for ILU(P) and ILU(T,P)
UMF ! Coarsest-level subsolver: ILU UMF SLU MUMPS SLUDIST
2 ! Coarsest-level fillin P for ILU(P) and ILU(T,P)
1.d-4 ! Coarsest-level threshold T for ILU(T,P)
4 ! Number of Jacobi sweeps for BJAC/PJAC coarsest-level solver
0.0125d0 ! Smoothed aggregation threshold: >= 0.0

Write Preview
Loading…
Cancel
Save