diff --git a/base/tools/psb_cspins.f90 b/base/tools/psb_cspins.f90 index ef0407c2..5cceb8a0 100644 --- a/base/tools/psb_cspins.f90 +++ b/base/tools/psb_cspins.f90 @@ -160,6 +160,8 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) else if (desc_a%is_asb()) then + nrow = desc_a%get_local_rows() + ncol = desc_a%get_local_cols() if (local_) then call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -180,8 +182,6 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) jla(1:nz) = ja(1:nz) call psb_glob_to_loc(ila(1:nz),desc_a,info,iact='I') call psb_glob_to_loc(jla(1:nz),desc_a,info,iact='I') - nrow = desc_a%get_local_rows() - ncol = desc_a%get_local_cols() call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then diff --git a/base/tools/psb_dspins.f90 b/base/tools/psb_dspins.f90 index b7a9bc7d..dc53a6f9 100644 --- a/base/tools/psb_dspins.f90 +++ b/base/tools/psb_dspins.f90 @@ -160,6 +160,8 @@ subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) else if (desc_a%is_asb()) then + nrow = desc_a%get_local_rows() + ncol = desc_a%get_local_cols() if (local_) then call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -180,8 +182,6 @@ subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) jla(1:nz) = ja(1:nz) call psb_glob_to_loc(ila(1:nz),desc_a,info,iact='I') call psb_glob_to_loc(jla(1:nz),desc_a,info,iact='I') - nrow = desc_a%get_local_rows() - ncol = desc_a%get_local_cols() call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then diff --git a/base/tools/psb_sspins.f90 b/base/tools/psb_sspins.f90 index 85ee0a27..841ba52c 100644 --- a/base/tools/psb_sspins.f90 +++ b/base/tools/psb_sspins.f90 @@ -160,6 +160,8 @@ subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) else if (desc_a%is_asb()) then + nrow = desc_a%get_local_rows() + ncol = desc_a%get_local_cols() if (local_) then call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -180,8 +182,6 @@ subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) jla(1:nz) = ja(1:nz) call psb_glob_to_loc(ila(1:nz),desc_a,info,iact='I') call psb_glob_to_loc(jla(1:nz),desc_a,info,iact='I') - nrow = desc_a%get_local_rows() - ncol = desc_a%get_local_cols() call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then diff --git a/base/tools/psb_zspins.f90 b/base/tools/psb_zspins.f90 index ecdcb311..c661216b 100644 --- a/base/tools/psb_zspins.f90 +++ b/base/tools/psb_zspins.f90 @@ -160,6 +160,8 @@ subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) else if (desc_a%is_asb()) then + nrow = desc_a%get_local_rows() + ncol = desc_a%get_local_cols() if (local_) then call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -180,8 +182,6 @@ subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) jla(1:nz) = ja(1:nz) call psb_glob_to_loc(ila(1:nz),desc_a,info,iact='I') call psb_glob_to_loc(jla(1:nz),desc_a,info,iact='I') - nrow = desc_a%get_local_rows() - ncol = desc_a%get_local_cols() call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then diff --git a/test/pargen/runs/ppde.inp b/test/pargen/runs/ppde.inp index 9cc83f55..cd642036 100644 --- a/test/pargen/runs/ppde.inp +++ b/test/pargen/runs/ppde.inp @@ -2,7 +2,7 @@ BICGSTAB Iterative method BICGSTAB CGS BICG BICGSTABL RGMRES BJAC Preconditioner NONE DIAG BJAC CSR Storage format for matrix A: CSR COO JAD -008 Domain size (acutal system is this**3) +040 Domain size (acutal system is this**3) 2 Stopping criterion 1000 MAXIT -2 ITRACE