fdurastante/amg4psblas
1
0
Fork 0
mirror of https://github.com/sfilippone/amg4psblas.git synced 2026-10-06 14:44:55 +00:00
Code Issues Packages Projects Releases Wiki Activity

Added BJAC_ITRACE and BJAC_RESCHECK options

Browse Source
This commit is contained in:
Cirdans-Home
2021-03-29 16:44:02 +02:00
parent 718b574519
commit baa4e78626
4 changed files with 104 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
amgprec/impl/amg_ccprecset.F90
+26
View File
@@ -193,6 +193,24 @@ subroutine amg_ccprecseti(p,what,val,info,ilev,ilmax,pos,idx)
end if end if
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
case('BJAC_ITRACE')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
case('BJAC_RESCHECK')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
case default case default
do il=ilev_, ilmax_ do il=ilev_, ilmax_
call p%precv(il)%set(what,val,info,pos=pos,idx=idx) call p%precv(il)%set(what,val,info,pos=pos,idx=idx)
@@ -224,6 +242,14 @@ subroutine amg_ccprecseti(p,what,val,info,ilev,ilmax,pos,idx)
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
end if end if
case('BJAC_ITRACE')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
end if
case('BJAC_RESCHECK')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
end if
case default case default
do ilev_=1,nlev_ do ilev_=1,nlev_
call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx) call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx)
amgprec/impl/amg_dcprecset.F90
+26
View File
@@ -199,6 +199,24 @@ subroutine amg_dcprecseti(p,what,val,info,ilev,ilmax,pos,idx)
end if end if
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
case('BJAC_ITRACE')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
case('BJAC_RESCHECK')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
case default case default
do il=ilev_, ilmax_ do il=ilev_, ilmax_
call p%precv(il)%set(what,val,info,pos=pos,idx=idx) call p%precv(il)%set(what,val,info,pos=pos,idx=idx)
@@ -230,6 +248,14 @@ subroutine amg_dcprecseti(p,what,val,info,ilev,ilmax,pos,idx)
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
end if end if
case('BJAC_ITRACE')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
end if
case('BJAC_RESCHECK')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
end if
case default case default
do ilev_=1,nlev_ do ilev_=1,nlev_
call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx) call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx)
amgprec/impl/amg_scprecset.F90
+26
View File
@@ -193,6 +193,24 @@ subroutine amg_scprecseti(p,what,val,info,ilev,ilmax,pos,idx)
end if end if
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
case('BJAC_ITRACE')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
case('BJAC_RESCHECK')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
case default case default
do il=ilev_, ilmax_ do il=ilev_, ilmax_
call p%precv(il)%set(what,val,info,pos=pos,idx=idx) call p%precv(il)%set(what,val,info,pos=pos,idx=idx)
@@ -224,6 +242,14 @@ subroutine amg_scprecseti(p,what,val,info,ilev,ilmax,pos,idx)
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
end if end if
case('BJAC_ITRACE')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
end if
case('BJAC_RESCHECK')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
end if
case default case default
do ilev_=1,nlev_ do ilev_=1,nlev_
call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx) call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx)
amgprec/impl/amg_zcprecset.F90
+26
View File
@@ -199,6 +199,24 @@ subroutine amg_zcprecseti(p,what,val,info,ilev,ilmax,pos,idx)
end if end if
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
case('BJAC_ITRACE')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
case('BJAC_RESCHECK')
if (ilev_ /= nlev_) then
write(psb_err_unit,*) name,&
& ': Error: Inconsistent specification of WHAT vs. ILEV'
info = -2
return
end if
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
case default case default
do il=ilev_, ilmax_ do il=ilev_, ilmax_
call p%precv(il)%set(what,val,info,pos=pos,idx=idx) call p%precv(il)%set(what,val,info,pos=pos,idx=idx)
@@ -230,6 +248,14 @@ subroutine amg_zcprecseti(p,what,val,info,ilev,ilmax,pos,idx)
call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos) call p%precv(nlev_)%set('SUB_FILLIN',val,info,pos=pos)
end if end if
case('BJAC_ITRACE')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_ITRACE',val,info,pos=pos)
end if
case('BJAC_RESCHECK')
if (nlev_ > 1) then
call p%precv(nlev_)%set('SMOOTHER_RESIDUAL',val,info,pos=pos)
end if
case default case default
do ilev_=1,nlev_ do ilev_=1,nlev_
call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx) call p%precv(ilev_)%set(what,val,info,pos=pos,idx=idx)