fdurastante/psblas3
1
0
Fork 0
mirror of https://github.com/sfilippone/psblas3.git synced 2026-10-17 19:32:08 +00:00
Code Issues Packages Projects Releases Wiki Activity

psblas3:

Browse Source 
base/serial/psb_spdot_srtd.f90

Fixed infinite loop bug in sparse AXPBY.
This commit is contained in:
Salvatore Filippone
2012-11-26 16:12:43 +00:00
parent 084f6d9f1a
commit 7ab5e1e8e0
1 changed files with 36 additions and 44 deletions
Download Patch File Download Diff File Expand all files Collapse all files
base/serial/psb_spdot_srtd.f90
+36 -44
View File
@@ -45,17 +45,15 @@ subroutine psb_s_nspaxpby(nz,iz,z,alpha, nx, ix, x, beta, ny,iy,y,info)
if (ipy > ny) exit
if (ix(ipx) == iy(ipy)) then
acc = beta*y(ipy) + alpha*x(ipx)
if (acc /= szero) then
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
end if
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
else
nz = nz + 1
if (nz > isz) then
@@ -175,17 +173,15 @@ subroutine psb_d_nspaxpby(nz,iz,z,alpha, nx, ix, x, beta, ny,iy,y,info)
if (ipy > ny) exit
if (ix(ipx) == iy(ipy)) then
acc = beta*y(ipy) + alpha*x(ipx)
if (acc /= dzero) then
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
end if
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
else
nz = nz + 1
if (nz > isz) then
@@ -321,17 +317,15 @@ subroutine psb_c_nspaxpby(nz,iz,z,alpha, nx, ix, x, beta, ny,iy,y,info)
if (ipy > ny) exit
if (ix(ipx) == iy(ipy)) then
acc = beta*y(ipy) + alpha*x(ipx)
if (acc /= czero) then
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
end if
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
else
nz = nz + 1
if (nz > isz) then
@@ -451,17 +445,15 @@ subroutine psb_z_nspaxpby(nz,iz,z,alpha, nx, ix, x, beta, ny,iy,y,info)
if (ipy > ny) exit
if (ix(ipx) == iy(ipy)) then
acc = beta*y(ipy) + alpha*x(ipx)
if (acc /= zzero) then
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
end if
nz = nz + 1
if (nz > isz) then
info = -1
return
endif
iz(nz) = ix(ipx)
z(nz) = acc
ipx = ipx + 1
ipy = ipy + 1
else
nz = nz + 1
if (nz > isz) then