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

Fix silly bug introduced by even more silly mistake.

Browse Source
This commit is contained in:
Salvatore Filippone
2021-09-27 08:55:16 -04:00
parent ba2c6aa721
commit 09928af671
2 changed files with 8 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
amgprec/amg_d_onelev_mod.f90
+1
View File
@@ -57,6 +57,7 @@ module amg_d_onelev_mod
use amg_d_base_smoother_mod
use amg_d_dec_aggregator_mod
use amg_d_parmatch_aggregator_mod
use amg_d_newmatch_aggregator_mod
use psb_base_mod, only : psb_dspmat_type, psb_d_vect_type, &
& psb_d_base_vect_type, psb_ldspmat_type, psb_dlinmap_type, psb_dpk_, &
amgprec/impl/level/amg_d_base_onelev_csetc.F90
+7 -2
View File
@@ -43,6 +43,7 @@ subroutine amg_d_base_onelev_csetc(lv,what,val,info,pos,idx)
use amg_d_dec_aggregator_mod
use amg_d_symdec_aggregator_mod
use amg_d_parmatch_aggregator_mod
use amg_d_newmatch_aggregator_mod
use amg_d_jac_smoother
use amg_d_as_smoother
use amg_d_diag_solver
@@ -252,8 +253,6 @@ subroutine amg_d_base_onelev_csetc(lv,what,val,info,pos,idx)
lv%parms%ml_cycle = amg_stringval(val)
case ('PAR_AGGR_ALG')
ival = amg_stringval(val)
lv%parms%par_aggr_alg = ival
if (allocated(lv%aggr)) then
call lv%aggr%free(info)
if (info == 0) deallocate(lv%aggr,stat=info)
@@ -263,6 +262,9 @@ subroutine amg_d_base_onelev_csetc(lv,what,val,info,pos,idx)
end if
end if
ival = amg_stringval(val)
lv%parms%par_aggr_alg = ival
select case(val)
case('DEC')
allocate(amg_d_dec_aggregator_type :: lv%aggr, stat=info)
@@ -270,9 +272,12 @@ subroutine amg_d_base_onelev_csetc(lv,what,val,info,pos,idx)
allocate(amg_d_symdec_aggregator_type :: lv%aggr, stat=info)
case('COUP','COUPLED')
allocate(amg_d_parmatch_aggregator_type :: lv%aggr, stat=info)
case('NEWMTC')
allocate(amg_d_newmatch_aggregator_type :: lv%aggr, stat=info)
case default
info = psb_err_internal_error_
end select
if (info == psb_success_) call lv%aggr%default()
case ('AGGR_ORD')