diff --git a/base/tools/psb_dglobtril.F90 b/base/tools/psb_dglobtril.F90 index 3fb7617a4..c4dcd819c 100644 --- a/base/tools/psb_dglobtril.F90 +++ b/base/tools/psb_dglobtril.F90 @@ -90,7 +90,7 @@ Subroutine psb_dglobtril(a,desc_a,b,info,diag,imin,imax,jmin,jmax) call ldglob%tril(ldtril,info,& & diag=diag,imin=imin,imax=imax,jmin=jmin,jmax=jmax) call ldglob%free() - call ldglob%mv_to(ldtcoo) + call ldtril%mv_to(ldtcoo) lnz = ldtcoo%get_nzeros() call desc_a%g2lip(ldtcoo%ia(1:lnz),info,owned=.false.) call desc_a%g2lip(ldtcoo%ja(1:lnz),info,owned=.false.) diff --git a/base/tools/psb_dglobtriu.F90 b/base/tools/psb_dglobtriu.F90 index d784a51f4..34a52da72 100644 --- a/base/tools/psb_dglobtriu.F90 +++ b/base/tools/psb_dglobtriu.F90 @@ -90,7 +90,7 @@ Subroutine psb_dglobtriu(a,desc_a,b,info,diag,imin,imax,jmin,jmax) call ldglob%triu(ldtriu,info,& & diag=diag,imin=imin,imax=imax,jmin=jmin,jmax=jmax) call ldglob%free() - call ldglob%mv_to(ldtcoo) + call ldtriu%mv_to(ldtcoo) lnz = ldtcoo%get_nzeros() call desc_a%g2lip(ldtcoo%ia(1:lnz),info,owned=.false.) call desc_a%g2lip(ldtcoo%ja(1:lnz),info,owned=.false.) diff --git a/test/pdegen/psb_d_pde3d.F90 b/test/pdegen/psb_d_pde3d.F90 index 982e582cf..3a1ff94c5 100644 --- a/test/pdegen/psb_d_pde3d.F90 +++ b/test/pdegen/psb_d_pde3d.F90 @@ -686,7 +686,7 @@ program psb_d_pde3d real(psb_dpk_) :: t1, t2, tprec ! sparse matrix and preconditioner - type(psb_dspmat_type) :: a + type(psb_dspmat_type) :: a, atril type(psb_dprec_type) :: prec ! descriptor type(psb_desc_type) :: desc_a @@ -881,6 +881,17 @@ program psb_d_pde3d write(psb_out_unit,'("Storage format for DESC_A: ",a)') desc_a%get_fmt() end if + block + integer(psb_lpk_), allocatable :: iv(:) + call psb_globtril(a,desc_a,atril,info) + write(fname,'(a,i0,a)') "amat-loc-",iam,".mtx" + call a%print(fname) + iv = desc_a%get_global_indices(owned=.false.) + write(fname,'(a,i0,a)') "amat-glb-",iam,".mtx" + call a%print(fname,iv=iv) + end block + + ! ! cleanup storage and exit