diff --git a/mlprec/mld_cprecinit.F90 b/mlprec/mld_cprecinit.F90 index 79164a46..aed40bd5 100644 --- a/mlprec/mld_cprecinit.F90 +++ b/mlprec/mld_cprecinit.F90 @@ -224,6 +224,7 @@ subroutine mld_cprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 1 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.e0/3.e0 + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = szero end do ilev_ = nlev_ if (info == 0) call psb_realloc(mld_ifpsz_,p%baseprecv(ilev_)%iprcparm,info) @@ -251,7 +252,7 @@ subroutine mld_cprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 4 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.e0/3.e0 - + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = szero case default write(0,*) name,': Warning: Unknown preconditioner type request "',ptype,'"' info = 2 diff --git a/mlprec/mld_dprecinit.F90 b/mlprec/mld_dprecinit.F90 index e6f128e4..3559be45 100644 --- a/mlprec/mld_dprecinit.F90 +++ b/mlprec/mld_dprecinit.F90 @@ -224,6 +224,7 @@ subroutine mld_dprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 1 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.d0/3.d0 + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = dzero end do ilev_ = nlev_ if (info == 0) call psb_realloc(mld_ifpsz_,p%baseprecv(ilev_)%iprcparm,info) @@ -253,7 +254,8 @@ subroutine mld_dprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 4 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.d0/3.d0 - + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = dzero + case default write(0,*) name,': Warning: Unknown preconditioner type request "',ptype,'"' info = 2 diff --git a/mlprec/mld_sprecinit.F90 b/mlprec/mld_sprecinit.F90 index 32be91b4..38069712 100644 --- a/mlprec/mld_sprecinit.F90 +++ b/mlprec/mld_sprecinit.F90 @@ -224,6 +224,7 @@ subroutine mld_sprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 1 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.e0/3.e0 + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = szero end do ilev_ = nlev_ if (info == 0) call psb_realloc(mld_ifpsz_,p%baseprecv(ilev_)%iprcparm,info) @@ -251,6 +252,7 @@ subroutine mld_sprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 4 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.e0/3.e0 + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = szero case default write(0,*) name,': Warning: Unknown preconditioner type request "',ptype,'"' diff --git a/mlprec/mld_zprecinit.F90 b/mlprec/mld_zprecinit.F90 index 7eb03596..91c53078 100644 --- a/mlprec/mld_zprecinit.F90 +++ b/mlprec/mld_zprecinit.F90 @@ -224,6 +224,7 @@ subroutine mld_zprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 1 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.d0/3.d0 + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = dzero end do ilev_ = nlev_ if (info == 0) call psb_realloc(mld_ifpsz_,p%baseprecv(ilev_)%iprcparm,info) @@ -253,6 +254,7 @@ subroutine mld_zprecinit(p,ptype,info,nlev) p%baseprecv(ilev_)%iprcparm(mld_sub_fillin_) = 0 p%baseprecv(ilev_)%iprcparm(mld_smoother_sweeps_) = 4 p%baseprecv(ilev_)%rprcparm(mld_aggr_damp_) = 4.d0/3.d0 + p%baseprecv(ilev_)%rprcparm(mld_aggr_thresh_) = dzero case default write(0,*) name,': Warning: Unknown preconditioner type request "',ptype,'"'