MLPREC: merged bjac_aply inside as_aply, and updated baseprec_aply

accordingly.
stopcriterion
Salvatore Filippone 17 years ago
parent 5a0956030d
commit 72f98f8f3a

@ -108,6 +108,15 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
select case(prec%iprcparm(mld_prec_type_)) select case(prec%iprcparm(mld_prec_type_))
case(mld_bjac_)
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if (info /= 0) then
info=4010
ch_err='mld_bjac_aply'
goto 9999
end if
case(mld_as_) case(mld_as_)
! !
! Additive Schwarz preconditioner ! Additive Schwarz preconditioner
@ -120,7 +129,7 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info) call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='psb_bjac_aply' ch_err='mld_bjac_aply'
goto 9999 goto 9999
end if end if

@ -155,19 +155,7 @@ subroutine mld_dbaseprec_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
end if end if
end if end if
case(mld_bjac_) case(mld_bjac_,mld_as_)
!
! Block-Jacobi preconditioner
!
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if(info /= 0) then
info=4010
ch_err='mld_bjac_aply'
goto 9999
end if
case(mld_as_)
! !
! Additive Schwarz preconditioner ! Additive Schwarz preconditioner
! !

@ -108,6 +108,15 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
select case(prec%iprcparm(mld_prec_type_)) select case(prec%iprcparm(mld_prec_type_))
case(mld_bjac_)
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if (info /= 0) then
info=4010
ch_err='mld_bjac_aply'
goto 9999
end if
case(mld_as_) case(mld_as_)
! !
! Additive Schwarz preconditioner ! Additive Schwarz preconditioner
@ -120,7 +129,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info) call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='psb_bjac_aply' ch_err='mld_bjac_aply'
goto 9999 goto 9999
end if end if

@ -159,19 +159,7 @@ subroutine mld_zbaseprec_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
end if end if
end if end if
case(mld_bjac_) case(mld_bjac_,mld_as_)
!
! Block-Jacobi preconditioner
!
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if(info /= 0) then
info=4010
ch_err='mld_bjac_aply'
goto 9999
end if
case(mld_as_)
! !
! Additive Schwarz preconditioner ! Additive Schwarz preconditioner
! !

Loading…
Cancel
Save