fdurastante
/
psblas3
mirror of https://github.com/sfilippone/psblas3.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Silly bug in coo insert

Browse Source
development
sfilippone 1 year ago
parent 492b28f342
commit 8633e76cb0
4 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File

2
base/serial/impl/psb_c_coo_impl.F90
Unescape Escape View File

@ -2869,6 +2869,7 @@ subroutine psb_c_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
! Hence the call to set_nzeros done here. ! Hence the call to set_nzeros done here.
!$omp critical !$omp critical
nza = a%get_nzeros() nza = a%get_nzeros()
nzaold = nza
isza = a%get_size() isza = a%get_size()
! Build phase. Must handle reallocations in a sensible way. ! Build phase. Must handle reallocations in a sensible way.
if (isza < (nza+nz)) then if (isza < (nza+nz)) then
@ -2879,7 +2880,6 @@ subroutine psb_c_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
info = psb_err_alloc_dealloc_; call psb_errpush(info,name) info = psb_err_alloc_dealloc_; call psb_errpush(info,name)
else else
#if defined(OPENMP) #if defined(OPENMP)
nzaold = nza
nza = nza + nz nza = nza + nz
#endif #endif
call a%set_nzeros(nza) call a%set_nzeros(nza)

2
base/serial/impl/psb_d_coo_impl.F90
Unescape Escape View File

@ -2869,6 +2869,7 @@ subroutine psb_d_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
! Hence the call to set_nzeros done here. ! Hence the call to set_nzeros done here.
!$omp critical !$omp critical
nza = a%get_nzeros() nza = a%get_nzeros()
nzaold = nza
isza = a%get_size() isza = a%get_size()
! Build phase. Must handle reallocations in a sensible way. ! Build phase. Must handle reallocations in a sensible way.
if (isza < (nza+nz)) then if (isza < (nza+nz)) then
@ -2879,7 +2880,6 @@ subroutine psb_d_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
info = psb_err_alloc_dealloc_; call psb_errpush(info,name) info = psb_err_alloc_dealloc_; call psb_errpush(info,name)
else else
#if defined(OPENMP) #if defined(OPENMP)
nzaold = nza
nza = nza + nz nza = nza + nz
#endif #endif
call a%set_nzeros(nza) call a%set_nzeros(nza)

2
base/serial/impl/psb_s_coo_impl.F90
Unescape Escape View File

@ -2869,6 +2869,7 @@ subroutine psb_s_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
! Hence the call to set_nzeros done here. ! Hence the call to set_nzeros done here.
!$omp critical !$omp critical
nza = a%get_nzeros() nza = a%get_nzeros()
nzaold = nza
isza = a%get_size() isza = a%get_size()
! Build phase. Must handle reallocations in a sensible way. ! Build phase. Must handle reallocations in a sensible way.
if (isza < (nza+nz)) then if (isza < (nza+nz)) then
@ -2879,7 +2880,6 @@ subroutine psb_s_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
info = psb_err_alloc_dealloc_; call psb_errpush(info,name) info = psb_err_alloc_dealloc_; call psb_errpush(info,name)
else else
#if defined(OPENMP) #if defined(OPENMP)
nzaold = nza
nza = nza + nz nza = nza + nz
#endif #endif
call a%set_nzeros(nza) call a%set_nzeros(nza)

2
base/serial/impl/psb_z_coo_impl.F90
Unescape Escape View File

@ -2869,6 +2869,7 @@ subroutine psb_z_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
! Hence the call to set_nzeros done here. ! Hence the call to set_nzeros done here.
!$omp critical !$omp critical
nza = a%get_nzeros() nza = a%get_nzeros()
nzaold = nza
isza = a%get_size() isza = a%get_size()
! Build phase. Must handle reallocations in a sensible way. ! Build phase. Must handle reallocations in a sensible way.
if (isza < (nza+nz)) then if (isza < (nza+nz)) then
@ -2879,7 +2880,6 @@ subroutine psb_z_coo_csput_a(nz,ia,ja,val,a,imin,imax,jmin,jmax,info)
info = psb_err_alloc_dealloc_; call psb_errpush(info,name) info = psb_err_alloc_dealloc_; call psb_errpush(info,name)
else else
#if defined(OPENMP) #if defined(OPENMP)
nzaold = nza
nza = nza + nz nza = nza + nz
#endif #endif
call a%set_nzeros(nza) call a%set_nzeros(nza)

Write Preview
Loading…
Cancel
Save