diff --git a/base/internals/psi_bld_tmphalo.f90 b/base/internals/psi_bld_tmphalo.f90 index 3fa0f52f..1842f92e 100644 --- a/base/internals/psi_bld_tmphalo.f90 +++ b/base/internals/psi_bld_tmphalo.f90 @@ -100,8 +100,8 @@ subroutine psi_bld_tmphalo(desc,info) helem(i) = n_row+i ! desc%loc_to_glob(n_row+i) end do - call desc%l2gip(helem(1:nh),info) - call desc%fnd_owner(helem(1:nh),hproc,info) + call desc%indxmap%l2gip(helem(1:nh),info) + call desc%indxmap%fnd_owner(helem(1:nh),hproc,info) if (info /= psb_success_) then call psb_errpush(psb_err_from_subroutine_,name,a_err='fnd_owner') diff --git a/base/internals/psi_desc_index.F90 b/base/internals/psi_desc_index.F90 index 1549ccea..d08ea8ed 100644 --- a/base/internals/psi_desc_index.F90 +++ b/base/internals/psi_desc_index.F90 @@ -261,7 +261,7 @@ subroutine psi_desc_index(desc,index_in,dep_list,& end do else - call desc%l2g(index_in(i+1:i+nerv),& + call desc%indxmap%l2g(index_in(i+1:i+nerv),& & sndbuf(bsdindx(proc+1)+1:bsdindx(proc+1)+nerv),& & info) @@ -310,13 +310,13 @@ subroutine psi_desc_index(desc,index_in,dep_list,& i = i + 1 nerv = sdsz(proc+1) desc_index(i) = nerv - call desc%g2l(sndbuf(bsdindx(proc+1)+1:bsdindx(proc+1)+nerv),& + call desc%indxmap%g2l(sndbuf(bsdindx(proc+1)+1:bsdindx(proc+1)+nerv),& & desc_index(i+1:i+nerv),info) i = i + nerv + 1 nesd = rvsz(proc+1) desc_index(i) = nesd - call desc%g2l(rcvbuf(brvindx(proc+1)+1:brvindx(proc+1)+nesd),& + call desc%indxmap%g2l(rcvbuf(brvindx(proc+1)+1:brvindx(proc+1)+nesd),& & desc_index(i+1:i+nesd),info) i = i + nesd + 1 end do diff --git a/base/tools/psb_ccdbldext.F90 b/base/tools/psb_ccdbldext.F90 index 2b01b9be..817fdec4 100644 --- a/base/tools/psb_ccdbldext.F90 +++ b/base/tools/psb_ccdbldext.F90 @@ -259,7 +259,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype) Do j=0,n_elem_recv-1 idx = ovrlap(counter+psb_elem_recv_+j) - call desc_ov%l2g(idx,gidx,info) + call desc_ov%indxmap%l2g(idx,gidx,info) If (gidx < 0) then info=-3 call psb_errpush(info,name) @@ -529,7 +529,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype) & write(debug_unit,*) me,' ',trim(name),& & ': going for first idx_cnv', desc_ov%indxmap%get_state() - call desc_ov%g2l(workr(1:iszr),maskr(1:iszr),info) + call desc_ov%indxmap%g2l(workr(1:iszr),maskr(1:iszr),info) iszs = count(maskr(1:iszr)<=0) if (iszs > size(works)) call psb_realloc(iszs,works,info) j = 0 @@ -559,7 +559,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype) do i=1,iszs idx = works(i) n_col = desc_ov%get_local_cols() - call desc_ov%g2l_ins(idx,lidx,info) + call desc_ov%indxmap%g2l_ins(idx,lidx,info) if (desc_ov%get_local_cols() > n_col ) then ! ! This is a new index. Assigning a local index as diff --git a/base/tools/psb_cins.f90 b/base/tools/psb_cins.f90 index 7fad014f..651dae00 100644 --- a/base/tools/psb_cins.f90 +++ b/base/tools/psb_cins.f90 @@ -134,7 +134,7 @@ subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) case(psb_dupl_ovwrt_) @@ -280,7 +280,7 @@ subroutine psb_cins_vect(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if call x%ins(m,irl,val,dupl_,info) if (info /= 0) then @@ -400,7 +400,7 @@ subroutine psb_cins_vect_r2(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if do i=1,n @@ -569,7 +569,7 @@ subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) diff --git a/base/tools/psb_cspins.f90 b/base/tools/psb_cspins.f90 index 17ad9ff6..f936aa0e 100644 --- a/base/tools/psb_cspins.f90 +++ b/base/tools/psb_cspins.f90 @@ -128,8 +128,8 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_a%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_a%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) if (info /= psb_success_) then ierr(1) = info @@ -174,8 +174,8 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info) - call desc_a%g2l(ja(1:nz),jla(1:nz),info) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info) + call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info) call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -277,8 +277,8 @@ subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) goto 9999 end if - call desc_ar%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_ac%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) + call desc_ar%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_ac%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) if (psb_errstatus_fatal()) then ierr(1) = info diff --git a/base/tools/psb_dcdbldext.F90 b/base/tools/psb_dcdbldext.F90 index 12b33ed8..d1563db8 100644 --- a/base/tools/psb_dcdbldext.F90 +++ b/base/tools/psb_dcdbldext.F90 @@ -259,7 +259,7 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype) Do j=0,n_elem_recv-1 idx = ovrlap(counter+psb_elem_recv_+j) - call desc_ov%l2g(idx,gidx,info) + call desc_ov%indxmap%l2g(idx,gidx,info) If (gidx < 0) then info=-3 call psb_errpush(info,name) @@ -529,7 +529,7 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype) & write(debug_unit,*) me,' ',trim(name),& & ': going for first idx_cnv', desc_ov%indxmap%get_state() - call desc_ov%g2l(workr(1:iszr),maskr(1:iszr),info) + call desc_ov%indxmap%g2l(workr(1:iszr),maskr(1:iszr),info) iszs = count(maskr(1:iszr)<=0) if (iszs > size(works)) call psb_realloc(iszs,works,info) j = 0 @@ -559,7 +559,7 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype) do i=1,iszs idx = works(i) n_col = desc_ov%get_local_cols() - call desc_ov%g2l_ins(idx,lidx,info) + call desc_ov%indxmap%g2l_ins(idx,lidx,info) if (desc_ov%get_local_cols() > n_col ) then ! ! This is a new index. Assigning a local index as diff --git a/base/tools/psb_dins.f90 b/base/tools/psb_dins.f90 index eb2ccb9e..385f0575 100644 --- a/base/tools/psb_dins.f90 +++ b/base/tools/psb_dins.f90 @@ -134,7 +134,7 @@ subroutine psb_dinsvi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) case(psb_dupl_ovwrt_) @@ -280,7 +280,7 @@ subroutine psb_dins_vect(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if call x%ins(m,irl,val,dupl_,info) if (info /= 0) then @@ -400,7 +400,7 @@ subroutine psb_dins_vect_r2(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if do i=1,n @@ -569,7 +569,7 @@ subroutine psb_dinsi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) diff --git a/base/tools/psb_dspins.f90 b/base/tools/psb_dspins.f90 index af05d760..56a3be17 100644 --- a/base/tools/psb_dspins.f90 +++ b/base/tools/psb_dspins.f90 @@ -128,8 +128,8 @@ subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_a%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_a%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) if (info /= psb_success_) then ierr(1) = info @@ -174,8 +174,8 @@ subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info) - call desc_a%g2l(ja(1:nz),jla(1:nz),info) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info) + call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info) call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -277,8 +277,8 @@ subroutine psb_dspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) goto 9999 end if - call desc_ar%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_ac%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) + call desc_ar%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_ac%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) if (psb_errstatus_fatal()) then ierr(1) = info diff --git a/base/tools/psb_loc_to_glob.f90 b/base/tools/psb_loc_to_glob.f90 index 0255ef0d..67e07409 100644 --- a/base/tools/psb_loc_to_glob.f90 +++ b/base/tools/psb_loc_to_glob.f90 @@ -66,6 +66,12 @@ subroutine psb_loc_to_glob2v(x,y,desc_a,info,iact) info=psb_success_ name='psb_loc_to_glob2' call psb_erractionsave(err_act) + if (.not.desc_a%is_valid()) then + info = psb_err_invalid_cd_state_ + call psb_errpush(info,name) + goto 9999 + endif + if (present(iact)) then act=iact @@ -171,6 +177,11 @@ subroutine psb_loc_to_glob1v(x,desc_a,info,iact) info=psb_success_ name='psb_loc_to_glob' call psb_erractionsave(err_act) + if (.not.desc_a%is_valid()) then + info = psb_err_invalid_cd_state_ + call psb_errpush(info,name) + goto 9999 + endif if (present(iact)) then act=iact diff --git a/base/tools/psb_scdbldext.F90 b/base/tools/psb_scdbldext.F90 index cfec406b..fae47496 100644 --- a/base/tools/psb_scdbldext.F90 +++ b/base/tools/psb_scdbldext.F90 @@ -259,7 +259,7 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype) Do j=0,n_elem_recv-1 idx = ovrlap(counter+psb_elem_recv_+j) - call desc_ov%l2g(idx,gidx,info) + call desc_ov%indxmap%l2g(idx,gidx,info) If (gidx < 0) then info=-3 call psb_errpush(info,name) @@ -529,7 +529,7 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype) & write(debug_unit,*) me,' ',trim(name),& & ': going for first idx_cnv', desc_ov%indxmap%get_state() - call desc_ov%g2l(workr(1:iszr),maskr(1:iszr),info) + call desc_ov%indxmap%g2l(workr(1:iszr),maskr(1:iszr),info) iszs = count(maskr(1:iszr)<=0) if (iszs > size(works)) call psb_realloc(iszs,works,info) j = 0 @@ -559,7 +559,7 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype) do i=1,iszs idx = works(i) n_col = desc_ov%get_local_cols() - call desc_ov%g2l_ins(idx,lidx,info) + call desc_ov%indxmap%g2l_ins(idx,lidx,info) if (desc_ov%get_local_cols() > n_col ) then ! ! This is a new index. Assigning a local index as diff --git a/base/tools/psb_sins.f90 b/base/tools/psb_sins.f90 index 9fc964da..6023b626 100644 --- a/base/tools/psb_sins.f90 +++ b/base/tools/psb_sins.f90 @@ -134,7 +134,7 @@ subroutine psb_sinsvi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) case(psb_dupl_ovwrt_) @@ -280,7 +280,7 @@ subroutine psb_sins_vect(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if call x%ins(m,irl,val,dupl_,info) if (info /= 0) then @@ -400,7 +400,7 @@ subroutine psb_sins_vect_r2(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if do i=1,n @@ -569,7 +569,7 @@ subroutine psb_sinsi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) diff --git a/base/tools/psb_sspins.f90 b/base/tools/psb_sspins.f90 index a93eea52..b3110bbe 100644 --- a/base/tools/psb_sspins.f90 +++ b/base/tools/psb_sspins.f90 @@ -128,8 +128,8 @@ subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_a%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_a%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) if (info /= psb_success_) then ierr(1) = info @@ -174,8 +174,8 @@ subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info) - call desc_a%g2l(ja(1:nz),jla(1:nz),info) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info) + call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info) call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -277,8 +277,8 @@ subroutine psb_sspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) goto 9999 end if - call desc_ar%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_ac%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) + call desc_ar%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_ac%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) if (psb_errstatus_fatal()) then ierr(1) = info diff --git a/base/tools/psb_zcdbldext.F90 b/base/tools/psb_zcdbldext.F90 index 0b5a21a2..78211399 100644 --- a/base/tools/psb_zcdbldext.F90 +++ b/base/tools/psb_zcdbldext.F90 @@ -259,7 +259,7 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype) Do j=0,n_elem_recv-1 idx = ovrlap(counter+psb_elem_recv_+j) - call desc_ov%l2g(idx,gidx,info) + call desc_ov%indxmap%l2g(idx,gidx,info) If (gidx < 0) then info=-3 call psb_errpush(info,name) @@ -529,7 +529,7 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype) & write(debug_unit,*) me,' ',trim(name),& & ': going for first idx_cnv', desc_ov%indxmap%get_state() - call desc_ov%g2l(workr(1:iszr),maskr(1:iszr),info) + call desc_ov%indxmap%g2l(workr(1:iszr),maskr(1:iszr),info) iszs = count(maskr(1:iszr)<=0) if (iszs > size(works)) call psb_realloc(iszs,works,info) j = 0 @@ -559,7 +559,7 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype) do i=1,iszs idx = works(i) n_col = desc_ov%get_local_cols() - call desc_ov%g2l_ins(idx,lidx,info) + call desc_ov%indxmap%g2l_ins(idx,lidx,info) if (desc_ov%get_local_cols() > n_col ) then ! ! This is a new index. Assigning a local index as diff --git a/base/tools/psb_zins.f90 b/base/tools/psb_zins.f90 index 4a3426cc..6ded6f52 100644 --- a/base/tools/psb_zins.f90 +++ b/base/tools/psb_zins.f90 @@ -134,7 +134,7 @@ subroutine psb_zinsvi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) case(psb_dupl_ovwrt_) @@ -280,7 +280,7 @@ subroutine psb_zins_vect(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if call x%ins(m,irl,val,dupl_,info) if (info /= 0) then @@ -400,7 +400,7 @@ subroutine psb_zins_vect_r2(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if do i=1,n @@ -569,7 +569,7 @@ subroutine psb_zinsi(m, irw, val, x, desc_a, info, dupl,local) if (local_) then irl(1:m) = irw(1:m) else - call desc_a%g2l(irw(1:m),irl(1:m),info,owned=.true.) + call desc_a%indxmap%g2l(irw(1:m),irl(1:m),info,owned=.true.) end if select case(dupl_) diff --git a/base/tools/psb_zspins.f90 b/base/tools/psb_zspins.f90 index d514059d..5deb87ba 100644 --- a/base/tools/psb_zspins.f90 +++ b/base/tools/psb_zspins.f90 @@ -128,8 +128,8 @@ subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_a%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_a%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info,mask=(ila(1:nz)>0)) if (info /= psb_success_) then ierr(1) = info @@ -174,8 +174,8 @@ subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild,local) goto 9999 end if - call desc_a%g2l(ia(1:nz),ila(1:nz),info) - call desc_a%g2l(ja(1:nz),jla(1:nz),info) + call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info) + call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info) call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info) if (info /= psb_success_) then @@ -277,8 +277,8 @@ subroutine psb_zspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) goto 9999 end if - call desc_ar%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) - call desc_ac%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) + call desc_ar%indxmap%g2l(ia(1:nz),ila(1:nz),info,owned=.true.) + call desc_ac%indxmap%g2l_ins(ja(1:nz),jla(1:nz),info, mask=(ila(1:nz)>0)) if (psb_errstatus_fatal()) then ierr(1) = info