diff --git a/base/tools/psb_c_remap.F90 b/base/tools/psb_c_remap.F90 index 710f1719..5b0bfc82 100644 --- a/base/tools/psb_c_remap.F90 +++ b/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() diff --git a/base/tools/psb_d_remap.F90 b/base/tools/psb_d_remap.F90 index 384485a4..4045bf18 100644 --- a/base/tools/psb_d_remap.F90 +++ b/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() diff --git a/base/tools/psb_s_remap.F90 b/base/tools/psb_s_remap.F90 index 199317ff..47af42b8 100644 --- a/base/tools/psb_s_remap.F90 +++ b/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() diff --git a/base/tools/psb_z_remap.F90 b/base/tools/psb_z_remap.F90 index fa3b67c9..b98364dc 100644 --- a/base/tools/psb_z_remap.F90 +++ b/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() diff --git a/test/cdasb/psb_d_pde3d.f90 b/test/cdasb/psb_d_pde3d.f90 index f78dc293..3e947e25 100644 --- a/test/cdasb/psb_d_pde3d.f90 +++ b/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