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

mld2p4:

Browse Source 
mlprec/impl/mld_cprecinit.F90
 mlprec/impl/mld_cprecset.F90
 mlprec/impl/mld_dprecinit.F90
 mlprec/impl/mld_dprecset.F90
 mlprec/impl/mld_sprecinit.F90
 mlprec/impl/mld_sprecset.F90
 mlprec/impl/mld_zprecinit.F90
 mlprec/impl/mld_zprecset.F90

Description of defaults.
Single vs. double precision
stopcriterion
Salvatore Filippone 9 years ago
parent ca1c2616e5
commit c9dc1f2d8e
8 changed files with 24 additions and 24 deletions
Show Stats Download Patch File Download Diff File

2
mlprec/impl/mld_cprecinit.F90
Unescape Escape View File

@ -61,7 +61,7 @@
! !
! 'ML' - Multilevel hybrid preconditioner (additive on the ! 'ML' - Multilevel hybrid preconditioner (additive on the
! same level and multiplicative through the levels), ! same level and multiplicative through the levels),
! with 2 levels and post-smoothing only. RAS with ! with 2 levels, pre and post-smoothing, RAS with
! overlap 1 and ILU(0) on the local blocks is ! overlap 1 and ILU(0) on the local blocks is
! applied as post-smoother at each level, but the ! applied as post-smoother at each level, but the
! coarsest one; four sweeps of the block-Jacobi solver, ! coarsest one; four sweeps of the block-Jacobi solver,

10
mlprec/impl/mld_cprecset.F90
Unescape Escape View File

@ -192,6 +192,11 @@ subroutine mld_cprecseti(p,what,val,info,ilev)
select case (val) select case (val)
case(mld_bjac_) case(mld_bjac_)
call onelev_set_smoother(p%precv(nlev_),val,info) call onelev_set_smoother(p%precv(nlev_),val,info)
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)
call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info) call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info)
@ -286,6 +291,11 @@ subroutine mld_cprecseti(p,what,val,info,ilev)
select case (val) select case (val)
case(mld_bjac_) case(mld_bjac_)
call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info) call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info)
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)
call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info) call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info)

2
mlprec/impl/mld_dprecinit.F90
Unescape Escape View File

@ -61,7 +61,7 @@
! !
! 'ML' - Multilevel hybrid preconditioner (additive on the ! 'ML' - Multilevel hybrid preconditioner (additive on the
! same level and multiplicative through the levels), ! same level and multiplicative through the levels),
! with 2 levels and post-smoothing only. RAS with ! with 2 levels, pre and post-smoothing, RAS with
! overlap 1 and ILU(0) on the local blocks is ! overlap 1 and ILU(0) on the local blocks is
! applied as post-smoother at each level, but the ! applied as post-smoother at each level, but the
! coarsest one; four sweeps of the block-Jacobi solver, ! coarsest one; four sweeps of the block-Jacobi solver,

10
mlprec/impl/mld_dprecset.F90
Unescape Escape View File

@ -204,11 +204,6 @@ subroutine mld_dprecseti(p,what,val,info,ilev)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info) call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else #else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info) call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif #endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)
@ -310,11 +305,6 @@ subroutine mld_dprecseti(p,what,val,info,ilev)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info) call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else #else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info) call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif #endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)

2
mlprec/impl/mld_sprecinit.F90
Unescape Escape View File

@ -61,7 +61,7 @@
! !
! 'ML' - Multilevel hybrid preconditioner (additive on the ! 'ML' - Multilevel hybrid preconditioner (additive on the
! same level and multiplicative through the levels), ! same level and multiplicative through the levels),
! with 2 levels and post-smoothing only. RAS with ! with 2 levels, pre and post-smoothing, RAS with
! overlap 1 and ILU(0) on the local blocks is ! overlap 1 and ILU(0) on the local blocks is
! applied as post-smoother at each level, but the ! applied as post-smoother at each level, but the
! coarsest one; four sweeps of the block-Jacobi solver, ! coarsest one; four sweeps of the block-Jacobi solver,

10
mlprec/impl/mld_sprecset.F90
Unescape Escape View File

@ -192,6 +192,11 @@ subroutine mld_sprecseti(p,what,val,info,ilev)
select case (val) select case (val)
case(mld_bjac_) case(mld_bjac_)
call onelev_set_smoother(p%precv(nlev_),val,info) call onelev_set_smoother(p%precv(nlev_),val,info)
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)
call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info) call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info)
@ -286,6 +291,11 @@ subroutine mld_sprecseti(p,what,val,info,ilev)
select case (val) select case (val)
case(mld_bjac_) case(mld_bjac_)
call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info) call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info)
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)
call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info) call onelev_set_smoother(p%precv(nlev_),mld_bjac_,info)

2
mlprec/impl/mld_zprecinit.F90
Unescape Escape View File

@ -61,7 +61,7 @@
! !
! 'ML' - Multilevel hybrid preconditioner (additive on the ! 'ML' - Multilevel hybrid preconditioner (additive on the
! same level and multiplicative through the levels), ! same level and multiplicative through the levels),
! with 2 levels and post-smoothing only. RAS with ! with 2 levels, pre and post-smoothing, RAS with
! overlap 1 and ILU(0) on the local blocks is ! overlap 1 and ILU(0) on the local blocks is
! applied as post-smoother at each level, but the ! applied as post-smoother at each level, but the
! coarsest one; four sweeps of the block-Jacobi solver, ! coarsest one; four sweeps of the block-Jacobi solver,

10
mlprec/impl/mld_zprecset.F90
Unescape Escape View File

@ -204,11 +204,6 @@ subroutine mld_zprecseti(p,what,val,info,ilev)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info) call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else #else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info) call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif #endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)
@ -310,11 +305,6 @@ subroutine mld_zprecseti(p,what,val,info,ilev)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info) call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else #else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info) call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif
#if defined(HAVE_SLU_)
call onelev_set_solver(p%precv(nlev_),mld_slu_,info)
#else
call onelev_set_solver(p%precv(nlev_),mld_ilu_n_,info)
#endif #endif
call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info) call p%precv(nlev_)%set(mld_coarse_mat_,mld_distr_mat_,info)
case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_) case(mld_umf_, mld_slu_,mld_ilu_n_, mld_ilu_t_,mld_milu_n_)

Write Preview
Loading…
Cancel
Save