psblas3:

base/serial/f03/psb_d_csr_impl.f03 Fixed bug in colsum,aclsum.
14 years ago · 083573809d
parent e17ebeaf89
commit 083573809d
1 changed files with 2 additions and 2 deletions
--- a/base/serial/f03/psb_d_csr_impl.f03
+++ b/base/serial/f03/psb_d_csr_impl.f03
@ -1241,7 +1241,7 @@ subroutine psb_d_csr_colsum(d,a)
  do i=1, m
    do j=a%irp(i),a%irp(i+1)-1
      k = a%ja(j)
-      d(k) = d(k) + (a%val(k))
+      d(k) = d(k) + (a%val(j))
    end do
  end do

@ -1293,7 +1293,7 @@ subroutine psb_d_csr_aclsum(d,a)
  do i=1, m
    do j=a%irp(i),a%irp(i+1)-1
      k = a%ja(j)
-      d(k) = d(k) + abs(a%val(k))
+      d(k) = d(k) + abs(a%val(j))
    end do
  end do