Debug check (to be taken out later)

base/tools/psb_c_remap.F90

@@ -203,6 +203,7 @@ subroutine psb_c_remap(np_remap, desc_in, a_in, desc_out, isrc, nrsrc, a_out, in
         call psb_rcv(ictxt,nrsrc(ip),isrc(ip))
         nzl = nzl + nzsrc(ip)
       end do
+      write(0,*) rme,' Check on NR:',newnl(rme+1),sum(nrsrc)
       call acoo_rcv%allocate(newnl(rme+1),newnl(rme+1),nzl)
       nrl = acoo_rcv%get_nrows()
       ncl = acoo_rcv%get_ncols()

base/tools/psb_d_remap.F90

@@ -203,6 +203,7 @@ subroutine psb_d_remap(np_remap, desc_in, a_in, desc_out, isrc, nrsrc, a_out, in
         call psb_rcv(ictxt,nrsrc(ip),isrc(ip))
         nzl = nzl + nzsrc(ip)
       end do
+      write(0,*) rme,' Check on NR:',newnl(rme+1),sum(nrsrc)
       call acoo_rcv%allocate(newnl(rme+1),newnl(rme+1),nzl)
       nrl = acoo_rcv%get_nrows()
       ncl = acoo_rcv%get_ncols()

base/tools/psb_s_remap.F90

@@ -203,6 +203,7 @@ subroutine psb_s_remap(np_remap, desc_in, a_in, desc_out, isrc, nrsrc, a_out, in
         call psb_rcv(ictxt,nrsrc(ip),isrc(ip))
         nzl = nzl + nzsrc(ip)
       end do
+      write(0,*) rme,' Check on NR:',newnl(rme+1),sum(nrsrc)
       call acoo_rcv%allocate(newnl(rme+1),newnl(rme+1),nzl)
       nrl = acoo_rcv%get_nrows()
       ncl = acoo_rcv%get_ncols()

base/tools/psb_z_remap.F90

@@ -203,6 +203,7 @@ subroutine psb_z_remap(np_remap, desc_in, a_in, desc_out, isrc, nrsrc, a_out, in
         call psb_rcv(ictxt,nrsrc(ip),isrc(ip))
         nzl = nzl + nzsrc(ip)
       end do
+      write(0,*) rme,' Check on NR:',newnl(rme+1),sum(nrsrc)
       call acoo_rcv%allocate(newnl(rme+1),newnl(rme+1),nzl)
       nrl = acoo_rcv%get_nrows()
       ncl = acoo_rcv%get_ncols()

test/cdasb/psb_d_pde3d.f90

@@ -647,7 +647,7 @@ program psb_d_pde3d
   integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, ipart
   integer(psb_epk_) :: amatsize, precsize, descsize, d2size
   real(psb_dpk_)   :: err, eps
-
+  integer(psb_ipk_), allocatable :: isrc(:), nrsrc(:)
   ! other variables
   integer(psb_ipk_) :: info, i, rnp
   character(len=20) :: name,ch_err
@@ -701,7 +701,7 @@ program psb_d_pde3d
     if (iam == 0) write(0,*) 'Remapping from ',np,' to ',rnp
     flush(0)
     call psb_barrier(ictxt)
-    call   psb_remap(rnp,desc_blk,a,desc_rmp,aremap,info)
+    call   psb_remap(rnp,desc_blk,a,desc_rmp,isrc,nrsrc,aremap,info)
     flush(0)
     call psb_barrier(ictxt)
     if (iam == 0) write(0,*) '                   Info ',info