Avoid overflow

l1aggregation
sfilippone 8 months ago
parent 12356f65f6
commit 15bb4ac101

@ -275,7 +275,7 @@ subroutine amg_d_soc1_map_bld(iorder,theta,clean_zeros,a,desc_a,nlaggr,ilaggr,in
disjoint = all(ilaggr(icol(1:ip)) == -(nr+1)).or.(ip==0) disjoint = all(ilaggr(icol(1:ip)) == -(nr+1)).or.(ip==0)
if (disjoint) then if (disjoint) then
locnaggr(kk) = locnaggr(kk) + 1 locnaggr(kk) = locnaggr(kk) + 1
itmp = (bnds(kk)-1+locnaggr(kk))*nths+kk itmp = (bnds(kk)-1_psb_lpk_+locnaggr(kk))*nths+kk
if (itmp < (bnds(kk)-1+locnaggr(kk))) then if (itmp < (bnds(kk)-1+locnaggr(kk))) then
!$omp atomic update !$omp atomic update
info = max(12345678,info) info = max(12345678,info)

Loading…
Cancel
Save