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_cmlprec_aply.f90
 mlprec/impl/mld_dmlprec_aply.f90
 mlprec/impl/mld_smlprec_aply.f90
 mlprec/impl/mld_zmlprec_aply.f90

Fix error and status messages.
stopcriterion
Salvatore Filippone 8 years ago
parent a128383525
commit e8f62cb1af
4 changed files with 28 additions and 12 deletions
Show Stats Download Patch File Download Diff File

10
mlprec/impl/mld_cmlprec_aply.f90
Unescape Escape View File

@ -902,7 +902,7 @@ contains
name = 'inner_inner_add' name = 'inner_kcycle'
info = psb_success_ info = psb_success_
call psb_erractionsave(err_act) call psb_erractionsave(err_act)
debug_unit = psb_get_debug_unit() debug_unit = psb_get_debug_unit()
@ -917,7 +917,7 @@ contains
call psb_info(ictxt, me, np) call psb_info(ictxt, me, np)
if(debug_level > 1) then if(debug_level > 1) then
write(debug_unit,*) me,' inner_add at level ',level write(debug_unit,*) me,name,' at level ',level
end if end if
if ((level<1).or.(level>nlev)) then if ((level<1).or.(level>nlev)) then
@ -992,11 +992,15 @@ contains
!Set the preconditioner !Set the preconditioner
if ((level < nlev - 2)) then if (level < nlev ) then
if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then
call mld_cinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG') call mld_cinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG')
elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then
call mld_cinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR') call mld_cinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR')
else
call psb_errpush(psb_err_internal_error_,name,&
& a_err='Bad value for ml_type')
goto 9999
endif endif
else else
call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info) call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info)

10
mlprec/impl/mld_dmlprec_aply.f90
Unescape Escape View File

@ -902,7 +902,7 @@ contains
name = 'inner_inner_add' name = 'inner_kcycle'
info = psb_success_ info = psb_success_
call psb_erractionsave(err_act) call psb_erractionsave(err_act)
debug_unit = psb_get_debug_unit() debug_unit = psb_get_debug_unit()
@ -917,7 +917,7 @@ contains
call psb_info(ictxt, me, np) call psb_info(ictxt, me, np)
if(debug_level > 1) then if(debug_level > 1) then
write(debug_unit,*) me,' inner_add at level ',level write(debug_unit,*) me,name,' at level ',level
end if end if
if ((level<1).or.(level>nlev)) then if ((level<1).or.(level>nlev)) then
@ -992,11 +992,15 @@ contains
!Set the preconditioner !Set the preconditioner
if ((level < nlev - 2)) then if (level < nlev ) then
if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then
call mld_dinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG') call mld_dinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG')
elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then
call mld_dinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR') call mld_dinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR')
else
call psb_errpush(psb_err_internal_error_,name,&
& a_err='Bad value for ml_type')
goto 9999
endif endif
else else
call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info) call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info)

10
mlprec/impl/mld_smlprec_aply.f90
Unescape Escape View File

@ -902,7 +902,7 @@ contains
name = 'inner_inner_add' name = 'inner_kcycle'
info = psb_success_ info = psb_success_
call psb_erractionsave(err_act) call psb_erractionsave(err_act)
debug_unit = psb_get_debug_unit() debug_unit = psb_get_debug_unit()
@ -917,7 +917,7 @@ contains
call psb_info(ictxt, me, np) call psb_info(ictxt, me, np)
if(debug_level > 1) then if(debug_level > 1) then
write(debug_unit,*) me,' inner_add at level ',level write(debug_unit,*) me,name,' at level ',level
end if end if
if ((level<1).or.(level>nlev)) then if ((level<1).or.(level>nlev)) then
@ -992,11 +992,15 @@ contains
!Set the preconditioner !Set the preconditioner
if ((level < nlev - 2)) then if (level < nlev ) then
if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then
call mld_sinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG') call mld_sinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG')
elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then
call mld_sinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR') call mld_sinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR')
else
call psb_errpush(psb_err_internal_error_,name,&
& a_err='Bad value for ml_type')
goto 9999
endif endif
else else
call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info) call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info)

10
mlprec/impl/mld_zmlprec_aply.f90
Unescape Escape View File

@ -902,7 +902,7 @@ contains
name = 'inner_inner_add' name = 'inner_kcycle'
info = psb_success_ info = psb_success_
call psb_erractionsave(err_act) call psb_erractionsave(err_act)
debug_unit = psb_get_debug_unit() debug_unit = psb_get_debug_unit()
@ -917,7 +917,7 @@ contains
call psb_info(ictxt, me, np) call psb_info(ictxt, me, np)
if(debug_level > 1) then if(debug_level > 1) then
write(debug_unit,*) me,' inner_add at level ',level write(debug_unit,*) me,name,' at level ',level
end if end if
if ((level<1).or.(level>nlev)) then if ((level<1).or.(level>nlev)) then
@ -992,11 +992,15 @@ contains
!Set the preconditioner !Set the preconditioner
if ((level < nlev - 2)) then if (level < nlev ) then
if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then if (p%precv(level)%parms%ml_type == mld_kcyclesym_ml_) then
call mld_zinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG') call mld_zinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'FCG')
elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then elseif (p%precv(level)%parms%ml_type == mld_kcycle_ml_) then
call mld_zinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR') call mld_zinneritkcycle(p, mlprec_wrk, level + 1, trans, work, 'CGR')
else
call psb_errpush(psb_err_internal_error_,name,&
& a_err='Bad value for ml_type')
goto 9999
endif endif
else else
call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info) call inner_ml_aply(level + 1 ,p,mlprec_wrk,trans,work,info)

Write Preview
Loading…
Cancel
Save