Fix silly bug in eigenvalue estimate in aggrmat_smth.

stopcriterion
Salvatore Filippone 6 years ago
parent bfe8c55865
commit f4f8bcc5ff

@ -233,6 +233,7 @@ subroutine mld_caggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
if (parms%aggr_eig == mld_max_norm_) then
anorm = acsr3%spnmi()
call psb_amx(ictxt,anorm)
omega = 4.d0/(3.d0*anorm)
parms%aggr_omega_val = omega

@ -233,6 +233,7 @@ subroutine mld_daggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
if (parms%aggr_eig == mld_max_norm_) then
anorm = acsr3%spnmi()
call psb_amx(ictxt,anorm)
omega = 4.d0/(3.d0*anorm)
parms%aggr_omega_val = omega

@ -233,6 +233,7 @@ subroutine mld_saggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
if (parms%aggr_eig == mld_max_norm_) then
anorm = acsr3%spnmi()
call psb_amx(ictxt,anorm)
omega = 4.d0/(3.d0*anorm)
parms%aggr_omega_val = omega

@ -233,6 +233,7 @@ subroutine mld_zaggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
if (parms%aggr_eig == mld_max_norm_) then
anorm = acsr3%spnmi()
call psb_amx(ictxt,anorm)
omega = 4.d0/(3.d0*anorm)
parms%aggr_omega_val = omega

Loading…
Cancel
Save