Fixes for GLOBTRIU/GLOBTRIL. Still getting spurious warnings.

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.)

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.)

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