Do not reinit CD when coming from UPD state.

3 years ago · 42190f17b7
parent 71293a2b2a
commit 42190f17b7
5 changed files with 5103 additions and 6442 deletions
--- a/base/tools/psb_cspasb.f90
+++ b/base/tools/psb_cspasb.f90
@ -146,12 +146,6 @@ subroutine psb_cspasb(a,desc_a, info, afmt, upd, mold)
        call psb_remote_mat(a%rmta,desc_a,a_add,info)
        nz = a_add%get_nzeros()
 !!$        write(0,*) me,name,' Nz to be added',nz
        nzt = nz
        call psb_sum(ctxt,nzt)
        if (nzt>0) then
          allocate(ivm, mold=desc_a%v_halo_index%v)
          call psb_cd_reinit(desc_a, info)
        end if
        if (nz > 0) then
          !
          ! Should we check for new indices here?
--- a/base/tools/psb_dspasb.f90
+++ b/base/tools/psb_dspasb.f90
@ -146,12 +146,6 @@ subroutine psb_dspasb(a,desc_a, info, afmt, upd, mold)
        call psb_remote_mat(a%rmta,desc_a,a_add,info)
        nz = a_add%get_nzeros()
 !!$        write(0,*) me,name,' Nz to be added',nz
        nzt = nz
        call psb_sum(ctxt,nzt)
        if (nzt>0) then
          allocate(ivm, mold=desc_a%v_halo_index%v)
          call psb_cd_reinit(desc_a, info)
        end if
        if (nz > 0) then
          !
          ! Should we check for new indices here?
--- a/base/tools/psb_sspasb.f90
+++ b/base/tools/psb_sspasb.f90
@ -146,12 +146,6 @@ subroutine psb_sspasb(a,desc_a, info, afmt, upd, mold)
        call psb_remote_mat(a%rmta,desc_a,a_add,info)
        nz = a_add%get_nzeros()
 !!$        write(0,*) me,name,' Nz to be added',nz
        nzt = nz
        call psb_sum(ctxt,nzt)
        if (nzt>0) then
          allocate(ivm, mold=desc_a%v_halo_index%v)
          call psb_cd_reinit(desc_a, info)
        end if
        if (nz > 0) then
          !
          ! Should we check for new indices here?
--- a/base/tools/psb_zspasb.f90
+++ b/base/tools/psb_zspasb.f90
@ -146,12 +146,6 @@ subroutine psb_zspasb(a,desc_a, info, afmt, upd, mold)
        call psb_remote_mat(a%rmta,desc_a,a_add,info)
        nz = a_add%get_nzeros()
 !!$        write(0,*) me,name,' Nz to be added',nz
        nzt = nz
        call psb_sum(ctxt,nzt)
        if (nzt>0) then
          allocate(ivm, mold=desc_a%v_halo_index%v)
          call psb_cd_reinit(desc_a, info)
        end if
        if (nz > 0) then
          !
          ! Should we check for new indices here?
--- a/docs/psblas-3.8.pdf
+++ b/docs/psblas-3.8.pdf