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

Added KRM options for coarse solver

Browse Source
cmake
Fabio Durastante 1 year ago
parent 31644c972a
commit 5cf7d5b3c1
6 changed files with 232 additions and 26 deletions
Show Stats Download Patch File Download Diff File

58
samples/advanced/pdegen/amg_d_pde2d.F90
Unescape Escape View File

@ -188,6 +188,17 @@ program amg_d_pde2d
integer(psb_ipk_) :: cfill ! fill-in for incomplete LU factorization
real(psb_dpk_) :: cthres ! threshold for ILUT factorization
integer(psb_ipk_) :: cjswp ! sweeps for GS or JAC coarsest-lev subsolver
! settings for the Krylov method
character(len=16) :: krm_method ! Krylov method for coarsest level
character(len=16) :: krm_prec ! Preconditioner for coarsest level
character(len=16) :: krm_subsolve ! Subsolver for coarsest level
character(len=16) :: krm_global ! Is the solver global or local? TRUE or FALSE
real(psb_dpk_) :: krm_eps ! Stopping tolerance
integer(psb_ipk_) :: krm_irst ! Restart for Krylov method
integer(psb_ipk_) :: krm_istop ! Stopping criterion
integer(psb_ipk_) :: krm_itmax ! Maximum number of iterations
integer(psb_ipk_) :: krm_itrace ! Trace of the lower Krylov iterations
integer(psb_ipk_) :: krm_fillin ! Fill-in for incomplete LU factorization
! Dump data
logical :: dump = .false.
@ -424,13 +435,26 @@ program amg_d_pde2d
end if
call prec%set('coarse_solve', p_choice%csolve, info)
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_mat', p_choice%cmat, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
! Set for the case of a KRM solver
if (psb_toupper(p_choice%csolve) == 'KRM') then
call prec%set('krm_method', p_choice%krm_method, info)
call prec%set('krm_kprec', p_choice%krm_prec, info)
call prec%set('krm_sub_solve', p_choice%krm_subsolve, info)
call prec%set('krm_global', p_choice%krm_global, info)
call prec%set('krm_eps', p_choice%krm_eps, info)
call prec%set('krm_irst', p_choice%krm_irst, info)
call prec%set('krm_istopc', p_choice%krm_istop, info)
call prec%set('krm_itmax', p_choice%krm_itmax, info)
call prec%set('krm_itrace', p_choice%krm_itrace, info)
call prec%set('krm_fillin', p_choice%krm_fillin, info)
else
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
end if
end select
! build the preconditioner
@ -681,6 +705,17 @@ contains
call read_data(prec%cfill,inp_unit) ! fill-in for incompl LU
call read_data(prec%cthres,inp_unit) ! Threshold for ILUT
call read_data(prec%cjswp,inp_unit) ! sweeps for GS/JAC subsolver
! Krylov method for coarsest level
call read_data(prec%krm_method,inp_unit) ! Krylov method for coarsest level
call read_data(prec%krm_prec,inp_unit) ! Preconditioner for coarsest level
call read_data(prec%krm_subsolve,inp_unit) ! Subsolver for coarsest level
call read_data(prec%krm_global,inp_unit) ! Is the solver global or local? TRUE or FALSE
call read_data(prec%krm_eps,inp_unit) ! Stopping tolerance
call read_data(prec%krm_irst,inp_unit) ! Restart for Krylov method
call read_data(prec%krm_istop,inp_unit) ! Stopping criterion
call read_data(prec%krm_itmax,inp_unit) ! Maximum number of iterations
call read_data(prec%krm_itrace,inp_unit) ! Trace of the lower Krylov iterations
call read_data(prec%krm_fillin,inp_unit) ! Fill-in for incomplete LU factorization
! dump
call read_data(prec%dump,inp_unit) ! Dump on file?
call read_data(prec%dlmin,inp_unit) ! Minimum level to dump
@ -770,6 +805,17 @@ contains
call psb_bcast(ctxt,prec%cfill)
call psb_bcast(ctxt,prec%cthres)
call psb_bcast(ctxt,prec%cjswp)
! Krylov method for coarsest level (broadcast)
call psb_bcast(ctxt,prec%krm_method)
call psb_bcast(ctxt,prec%krm_prec)
call psb_bcast(ctxt,prec%krm_subsolve)
call psb_bcast(ctxt,prec%krm_global)
call psb_bcast(ctxt,prec%krm_eps)
call psb_bcast(ctxt,prec%krm_irst)
call psb_bcast(ctxt,prec%krm_istop)
call psb_bcast(ctxt,prec%krm_itmax)
call psb_bcast(ctxt,prec%krm_itrace)
call psb_bcast(ctxt,prec%krm_fillin)
! dump
call psb_bcast(ctxt,prec%dump)
call psb_bcast(ctxt,prec%dlmin)

58
samples/advanced/pdegen/amg_d_pde3d.F90
Unescape Escape View File

@ -189,6 +189,17 @@ program amg_d_pde3d
integer(psb_ipk_) :: cfill ! fill-in for incomplete LU factorization
real(psb_dpk_) :: cthres ! threshold for ILUT factorization
integer(psb_ipk_) :: cjswp ! sweeps for GS or JAC coarsest-lev subsolver
! settings for the Krylov method
character(len=16) :: krm_method ! Krylov method for coarsest level
character(len=16) :: krm_prec ! Preconditioner for coarsest level
character(len=16) :: krm_subsolve ! Subsolver for coarsest level
character(len=16) :: krm_global ! Is the solver global or local? TRUE or FALSE
real(psb_dpk_) :: krm_eps ! Stopping tolerance
integer(psb_ipk_) :: krm_irst ! Restart for Krylov method
integer(psb_ipk_) :: krm_istop ! Stopping criterion
integer(psb_ipk_) :: krm_itmax ! Maximum number of iterations
integer(psb_ipk_) :: krm_itrace ! Trace of the lower Krylov iterations
integer(psb_ipk_) :: krm_fillin ! Fill-in for incomplete LU factorization
! Dump data
logical :: dump = .false.
@ -428,13 +439,26 @@ program amg_d_pde3d
end if
call prec%set('coarse_solve', p_choice%csolve, info)
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_mat', p_choice%cmat, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
! Set for the case of a KRM solver
if (psb_toupper(p_choice%csolve) == 'KRM') then
call prec%set('krm_method', p_choice%krm_method, info)
call prec%set('krm_kprec', p_choice%krm_prec, info)
call prec%set('krm_sub_solve', p_choice%krm_subsolve, info)
call prec%set('krm_global', p_choice%krm_global, info)
call prec%set('krm_eps', p_choice%krm_eps, info)
call prec%set('krm_irst', p_choice%krm_irst, info)
call prec%set('krm_istopc', p_choice%krm_istop, info)
call prec%set('krm_itmax', p_choice%krm_itmax, info)
call prec%set('krm_itrace', p_choice%krm_itrace, info)
call prec%set('krm_fillin', p_choice%krm_fillin, info)
else
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
end if
end select
! build the preconditioner
@ -685,6 +709,17 @@ contains
call read_data(prec%cfill,inp_unit) ! fill-in for incompl LU
call read_data(prec%cthres,inp_unit) ! Threshold for ILUT
call read_data(prec%cjswp,inp_unit) ! sweeps for GS/JAC subsolver
! Krylov method for coarsest level
call read_data(prec%krm_method,inp_unit) ! Krylov method for coarsest level
call read_data(prec%krm_prec,inp_unit) ! Preconditioner for coarsest level
call read_data(prec%krm_subsolve,inp_unit) ! Subsolver for coarsest level
call read_data(prec%krm_global,inp_unit) ! Is the solver global or local? TRUE or FALSE
call read_data(prec%krm_eps,inp_unit) ! Stopping tolerance
call read_data(prec%krm_irst,inp_unit) ! Restart for Krylov method
call read_data(prec%krm_istop,inp_unit) ! Stopping criterion
call read_data(prec%krm_itmax,inp_unit) ! Maximum number of iterations
call read_data(prec%krm_itrace,inp_unit) ! Trace of the lower Krylov iterations
call read_data(prec%krm_fillin,inp_unit) ! Fill-in for incomplete LU factorization
! dump
call read_data(prec%dump,inp_unit) ! Dump on file?
call read_data(prec%dlmin,inp_unit) ! Minimum level to dump
@ -774,6 +809,17 @@ contains
call psb_bcast(ctxt,prec%cfill)
call psb_bcast(ctxt,prec%cthres)
call psb_bcast(ctxt,prec%cjswp)
! Krylov method for coarsest level (broadcast)
call psb_bcast(ctxt,prec%krm_method)
call psb_bcast(ctxt,prec%krm_prec)
call psb_bcast(ctxt,prec%krm_subsolve)
call psb_bcast(ctxt,prec%krm_global)
call psb_bcast(ctxt,prec%krm_eps)
call psb_bcast(ctxt,prec%krm_irst)
call psb_bcast(ctxt,prec%krm_istop)
call psb_bcast(ctxt,prec%krm_itmax)
call psb_bcast(ctxt,prec%krm_itrace)
call psb_bcast(ctxt,prec%krm_fillin)
! dump
call psb_bcast(ctxt,prec%dump)
call psb_bcast(ctxt,prec%dlmin)

58
samples/advanced/pdegen/amg_s_pde2d.F90
Unescape Escape View File

@ -188,6 +188,17 @@ program amg_s_pde2d
integer(psb_ipk_) :: cfill ! fill-in for incomplete LU factorization
real(psb_spk_) :: cthres ! threshold for ILUT factorization
integer(psb_ipk_) :: cjswp ! sweeps for GS or JAC coarsest-lev subsolver
! settings for the Krylov method
character(len=16) :: krm_method ! Krylov method for coarsest level
character(len=16) :: krm_prec ! Preconditioner for coarsest level
character(len=16) :: krm_subsolve ! Subsolver for coarsest level
character(len=16) :: krm_global ! Is the solver global or local? TRUE or FALSE
real(psb_spk_) :: krm_eps ! Stopping tolerance
integer(psb_ipk_) :: krm_irst ! Restart for Krylov method
integer(psb_ipk_) :: krm_istop ! Stopping criterion
integer(psb_ipk_) :: krm_itmax ! Maximum number of iterations
integer(psb_ipk_) :: krm_itrace ! Trace of the lower Krylov iterations
integer(psb_ipk_) :: krm_fillin ! Fill-in for incomplete LU factorization
! Dump data
logical :: dump = .false.
@ -424,13 +435,26 @@ program amg_s_pde2d
end if
call prec%set('coarse_solve', p_choice%csolve, info)
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_mat', p_choice%cmat, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
! Set for the case of a KRM solver
if (psb_toupper(p_choice%csolve) == 'KRM') then
call prec%set('krm_method', p_choice%krm_method, info)
call prec%set('krm_kprec', p_choice%krm_prec, info)
call prec%set('krm_sub_solve', p_choice%krm_subsolve, info)
call prec%set('krm_global', p_choice%krm_global, info)
call prec%set('krm_eps', p_choice%krm_eps, info)
call prec%set('krm_irst', p_choice%krm_irst, info)
call prec%set('krm_istopc', p_choice%krm_istop, info)
call prec%set('krm_itmax', p_choice%krm_itmax, info)
call prec%set('krm_itrace', p_choice%krm_itrace, info)
call prec%set('krm_fillin', p_choice%krm_fillin, info)
else
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
end if
end select
! build the preconditioner
@ -681,6 +705,17 @@ contains
call read_data(prec%cfill,inp_unit) ! fill-in for incompl LU
call read_data(prec%cthres,inp_unit) ! Threshold for ILUT
call read_data(prec%cjswp,inp_unit) ! sweeps for GS/JAC subsolver
! Krylov method for coarsest level
call read_data(prec%krm_method,inp_unit) ! Krylov method for coarsest level
call read_data(prec%krm_prec,inp_unit) ! Preconditioner for coarsest level
call read_data(prec%krm_subsolve,inp_unit) ! Subsolver for coarsest level
call read_data(prec%krm_global,inp_unit) ! Is the solver global or local? TRUE or FALSE
call read_data(prec%krm_eps,inp_unit) ! Stopping tolerance
call read_data(prec%krm_irst,inp_unit) ! Restart for Krylov method
call read_data(prec%krm_istop,inp_unit) ! Stopping criterion
call read_data(prec%krm_itmax,inp_unit) ! Maximum number of iterations
call read_data(prec%krm_itrace,inp_unit) ! Trace of the lower Krylov iterations
call read_data(prec%krm_fillin,inp_unit) ! Fill-in for incomplete LU factorization
! dump
call read_data(prec%dump,inp_unit) ! Dump on file?
call read_data(prec%dlmin,inp_unit) ! Minimum level to dump
@ -770,6 +805,17 @@ contains
call psb_bcast(ctxt,prec%cfill)
call psb_bcast(ctxt,prec%cthres)
call psb_bcast(ctxt,prec%cjswp)
! Krylov method for coarsest level (broadcast)
call psb_bcast(ctxt,prec%krm_method)
call psb_bcast(ctxt,prec%krm_prec)
call psb_bcast(ctxt,prec%krm_subsolve)
call psb_bcast(ctxt,prec%krm_global)
call psb_bcast(ctxt,prec%krm_eps)
call psb_bcast(ctxt,prec%krm_irst)
call psb_bcast(ctxt,prec%krm_istop)
call psb_bcast(ctxt,prec%krm_itmax)
call psb_bcast(ctxt,prec%krm_itrace)
call psb_bcast(ctxt,prec%krm_fillin)
! dump
call psb_bcast(ctxt,prec%dump)
call psb_bcast(ctxt,prec%dlmin)

58
samples/advanced/pdegen/amg_s_pde3d.F90
Unescape Escape View File

@ -189,6 +189,17 @@ program amg_s_pde3d
integer(psb_ipk_) :: cfill ! fill-in for incomplete LU factorization
real(psb_spk_) :: cthres ! threshold for ILUT factorization
integer(psb_ipk_) :: cjswp ! sweeps for GS or JAC coarsest-lev subsolver
! settings for the Krylov method
character(len=16) :: krm_method ! Krylov method for coarsest level
character(len=16) :: krm_prec ! Preconditioner for coarsest level
character(len=16) :: krm_subsolve ! Subsolver for coarsest level
character(len=16) :: krm_global ! Is the solver global or local? TRUE or FALSE
real(psb_spk_) :: krm_eps ! Stopping tolerance
integer(psb_ipk_) :: krm_irst ! Restart for Krylov method
integer(psb_ipk_) :: krm_istop ! Stopping criterion
integer(psb_ipk_) :: krm_itmax ! Maximum number of iterations
integer(psb_ipk_) :: krm_itrace ! Trace of the lower Krylov iterations
integer(psb_ipk_) :: krm_fillin ! Fill-in for incomplete LU factorization
! Dump data
logical :: dump = .false.
@ -428,13 +439,26 @@ program amg_s_pde3d
end if
call prec%set('coarse_solve', p_choice%csolve, info)
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_mat', p_choice%cmat, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
! Set for the case of a KRM solver
if (psb_toupper(p_choice%csolve) == 'KRM') then
call prec%set('krm_method', p_choice%krm_method, info)
call prec%set('krm_kprec', p_choice%krm_prec, info)
call prec%set('krm_sub_solve', p_choice%krm_subsolve, info)
call prec%set('krm_global', p_choice%krm_global, info)
call prec%set('krm_eps', p_choice%krm_eps, info)
call prec%set('krm_irst', p_choice%krm_irst, info)
call prec%set('krm_istopc', p_choice%krm_istop, info)
call prec%set('krm_itmax', p_choice%krm_itmax, info)
call prec%set('krm_itrace', p_choice%krm_itrace, info)
call prec%set('krm_fillin', p_choice%krm_fillin, info)
else
if (psb_toupper(p_choice%csolve) == 'BJAC') &
& call prec%set('coarse_subsolve', p_choice%csbsolve, info)
call prec%set('coarse_fillin', p_choice%cfill, info)
call prec%set('coarse_iluthrs', p_choice%cthres, info)
call prec%set('coarse_sweeps', p_choice%cjswp, info)
end if
end select
! build the preconditioner
@ -685,6 +709,17 @@ contains
call read_data(prec%cfill,inp_unit) ! fill-in for incompl LU
call read_data(prec%cthres,inp_unit) ! Threshold for ILUT
call read_data(prec%cjswp,inp_unit) ! sweeps for GS/JAC subsolver
! Krylov method for coarsest level
call read_data(prec%krm_method,inp_unit) ! Krylov method for coarsest level
call read_data(prec%krm_prec,inp_unit) ! Preconditioner for coarsest level
call read_data(prec%krm_subsolve,inp_unit) ! Subsolver for coarsest level
call read_data(prec%krm_global,inp_unit) ! Is the solver global or local? TRUE or FALSE
call read_data(prec%krm_eps,inp_unit) ! Stopping tolerance
call read_data(prec%krm_irst,inp_unit) ! Restart for Krylov method
call read_data(prec%krm_istop,inp_unit) ! Stopping criterion
call read_data(prec%krm_itmax,inp_unit) ! Maximum number of iterations
call read_data(prec%krm_itrace,inp_unit) ! Trace of the lower Krylov iterations
call read_data(prec%krm_fillin,inp_unit) ! Fill-in for incomplete LU factorization
! dump
call read_data(prec%dump,inp_unit) ! Dump on file?
call read_data(prec%dlmin,inp_unit) ! Minimum level to dump
@ -774,6 +809,17 @@ contains
call psb_bcast(ctxt,prec%cfill)
call psb_bcast(ctxt,prec%cthres)
call psb_bcast(ctxt,prec%cjswp)
! Krylov method for coarsest level (broadcast)
call psb_bcast(ctxt,prec%krm_method)
call psb_bcast(ctxt,prec%krm_prec)
call psb_bcast(ctxt,prec%krm_subsolve)
call psb_bcast(ctxt,prec%krm_global)
call psb_bcast(ctxt,prec%krm_eps)
call psb_bcast(ctxt,prec%krm_irst)
call psb_bcast(ctxt,prec%krm_istop)
call psb_bcast(ctxt,prec%krm_itmax)
call psb_bcast(ctxt,prec%krm_itrace)
call psb_bcast(ctxt,prec%krm_fillin)
! dump
call psb_bcast(ctxt,prec%dump)
call psb_bcast(ctxt,prec%dlmin)

13
samples/advanced/pdegen/runs/amg_pde2d.inp
Unescape Escape View File

@ -74,7 +74,18 @@ DIST ! Coarsest-level matrix distribution: DIST REPL
1 ! Coarsest-level fillin P for ILU(P) and ILU(T,P)
1.d-4 ! Coarsest-level threshold T for ILU(T,P)
1 ! Number of sweeps for JACOBI/GS/BJAC coarsest-level solver
%%%%%%%%%%% Dump parms %%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%% Krylov coarse solver options %%%%%%%%%%
FCG ! Krylov method for coarsest level
BJAC ! Preconditioner for coarsest level
ILU ! Subsolver for coarsest level
FALSE ! Is the solver global (TRUE) or local (FALSE)?
1.d-4 ! Tolerance
-1 ! Restart for Krylov method
2 ! Stopping criterion
30 ! Maximum number of iterations
-1 ! Trace of the lower Krylov iterations
0 ! Fill-in for incomplete LU factorization
%%%%%%%%%% Dump parms %%%%%%%%%%%%%%%%%%%%%%%%%%
F ! Dump preconditioner on file
1 ! Min level
20 ! Max level

13
samples/advanced/pdegen/runs/amg_pde3d.inp
Unescape Escape View File

@ -68,12 +68,23 @@ FILTER ! Filtering of matrix: FILTER NOFILTER
-2 ! Number of thresholds in vector, next line ignored if <= 0
0.05 0.025 ! Thresholds
%%%%%%%%%%% Coarse level solver %%%%%%%%%%%%%%%%
BJAC ! Coarsest-level solver: MUMPS UMF SLU SLUDIST JACOBI GS BJAC
BJAC ! Coarsest-level solver: MUMPS UMF SLU SLUDIST JACOBI GS BJAC KRM
ILU ! Coarsest-level subsolver for BJAC: ILU ILUT MILU UMF MUMPS SLU
DIST ! Coarsest-level matrix distribution: DIST REPL
1 ! Coarsest-level fillin P for ILU(P) and ILU(T,P)
1.d-4 ! Coarsest-level threshold T for ILU(T,P)
1 ! Number of sweeps for JACOBI/GS/BJAC coarsest-level solver
%%%%%%%%%% Krylov coarse solver options %%%%%%%%%%
FCG ! Krylov method for coarsest level
BJAC ! Preconditioner for coarsest level
ILU ! Subsolver for coarsest level
FALSE ! Is the solver global (TRUE) or local (FALSE)?
1.d-4 ! Tolerance
-1 ! Restart for Krylov method
2 ! Stopping criterion
30 ! Maximum number of iterations
-1 ! Trace of the lower Krylov iterations
0 ! Fill-in for incomplete LU factorization
%%%%%%%%%%% Dump parms %%%%%%%%%%%%%%%%%%%%%%%%%%
F ! Dump preconditioner on file
1 ! Min level

Write Preview
Loading…
Cancel
Save