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

Fix cmp_avg_cr

Browse Source
This commit is contained in:
sfilippone
2026-07-14 10:53:27 +02:00
parent 8ac5b90454
commit 9cf27f6c46
4 changed files with 28 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files
amgprec/amg_c_prec_type.f90
+7 -3
View File
@@ -572,13 +572,17 @@ contains
integer(psb_ipk_) :: il, nl, iam, np
avgcr = szero
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
nl = prec%get_nlevs()
do il=2,nl
avgcr = avgcr + max(szero,prec%precv(il)%szratio)
if (prec%precv(il)%base_desc%is_ok()) then
ctxt = prec%precv(il)%base_desc%get_ctxt()
call psb_info(ctxt,iam,np)
if (iam >=0) avgcr = avgcr + max(szero,prec%precv(il)%szratio)
end if
end do
avgcr = avgcr / (nl-1)
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
call psb_sum(ctxt,avgcr)
prec%ag_data%avg_cr = avgcr/np
end subroutine amg_c_cmp_avg_cr
amgprec/amg_d_prec_type.f90
+7 -3
View File
@@ -572,13 +572,17 @@ contains
integer(psb_ipk_) :: il, nl, iam, np
avgcr = dzero
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
nl = prec%get_nlevs()
do il=2,nl
avgcr = avgcr + max(dzero,prec%precv(il)%szratio)
if (prec%precv(il)%base_desc%is_ok()) then
ctxt = prec%precv(il)%base_desc%get_ctxt()
call psb_info(ctxt,iam,np)
if (iam >=0) avgcr = avgcr + max(dzero,prec%precv(il)%szratio)
end if
end do
avgcr = avgcr / (nl-1)
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
call psb_sum(ctxt,avgcr)
prec%ag_data%avg_cr = avgcr/np
end subroutine amg_d_cmp_avg_cr
amgprec/amg_s_prec_type.f90
+7 -3
View File
@@ -572,13 +572,17 @@ contains
integer(psb_ipk_) :: il, nl, iam, np
avgcr = szero
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
nl = prec%get_nlevs()
do il=2,nl
avgcr = avgcr + max(szero,prec%precv(il)%szratio)
if (prec%precv(il)%base_desc%is_ok()) then
ctxt = prec%precv(il)%base_desc%get_ctxt()
call psb_info(ctxt,iam,np)
if (iam >=0) avgcr = avgcr + max(szero,prec%precv(il)%szratio)
end if
end do
avgcr = avgcr / (nl-1)
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
call psb_sum(ctxt,avgcr)
prec%ag_data%avg_cr = avgcr/np
end subroutine amg_s_cmp_avg_cr
amgprec/amg_z_prec_type.f90
+7 -3
View File
@@ -572,13 +572,17 @@ contains
integer(psb_ipk_) :: il, nl, iam, np
avgcr = dzero
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
nl = prec%get_nlevs()
do il=2,nl
avgcr = avgcr + max(dzero,prec%precv(il)%szratio)
if (prec%precv(il)%base_desc%is_ok()) then
ctxt = prec%precv(il)%base_desc%get_ctxt()
call psb_info(ctxt,iam,np)
if (iam >=0) avgcr = avgcr + max(dzero,prec%precv(il)%szratio)
end if
end do
avgcr = avgcr / (nl-1)
ctxt = prec%ctxt
call psb_info(ctxt,iam,np)
call psb_sum(ctxt,avgcr)
prec%ag_data%avg_cr = avgcr/np
end subroutine amg_z_cmp_avg_cr