Fix intent in ilu/inv

prec/impl/psb_c_bjacprec_impl.f90

@@ -495,7 +495,7 @@ subroutine psb_c_bjac_precbld(a,desc_a,prec,info,amold,vmold,imold)
   use psb_c_bjacprec, psb_protect_name => psb_c_bjac_precbld
   Implicit None
 
-  type(psb_cspmat_type), intent(in), target :: a
+  type(psb_cspmat_type), intent(inout), target :: a
   type(psb_desc_type), intent(inout), target   :: desc_a
   class(psb_c_bjac_prec_type),intent(inout) :: prec
   integer(psb_ipk_), intent(out)                      :: info

prec/impl/psb_c_invk_fact.f90

@@ -41,7 +41,7 @@ subroutine psb_c_invk_bld(a,fill1, fill2,lmat,d,umat,desc,info,blck)
   implicit none
 
   ! Arguments
-  type(psb_cspmat_type), intent(in), target   :: a
+  type(psb_cspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fill1, fill2
   type(psb_cspmat_type), intent(inout)        :: lmat, umat
   complex(psb_spk_), allocatable                 :: d(:)
@@ -155,7 +155,7 @@ subroutine psb_csparse_invk(n,a,z,fill_in,info,inlevs)
   use psb_c_invk_fact_mod, psb_protect_name => psb_csparse_invk
 
   integer(psb_ipk_), intent(in)           :: n
-  type(psb_cspmat_type), intent(in)       :: a
+  type(psb_cspmat_type), intent(inout)   :: a
   type(psb_cspmat_type), intent(inout)    :: z
   integer(psb_ipk_), intent(in)           :: fill_in
   integer(psb_ipk_), intent(out)          :: info

prec/impl/psb_c_invt_fact.f90

@@ -73,7 +73,7 @@ subroutine psb_c_invt_bld(a,fillin,invfill,thresh,invthresh,&
   implicit none
 
   ! Arguments
-  type(psb_cspmat_type), intent(in), target   :: a
+  type(psb_cspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fillin,invfill
   real(psb_spk_), intent(in)                :: thresh
   real(psb_spk_), intent(in)                  :: invthresh
@@ -198,7 +198,7 @@ subroutine psb_csparse_invt(n,a,z,nzrmax,sp_thresh,info)
 
   implicit none
   integer(psb_ipk_), intent(in)        :: n
-  type(psb_cspmat_type), intent(in)    :: a
+  type(psb_cspmat_type), intent(inout) :: a
   type(psb_cspmat_type), intent(inout) :: z
   integer(psb_ipk_), intent(in)        :: nzrmax
   real(psb_spk_), intent(in)           :: sp_thresh

prec/impl/psb_d_bjacprec_impl.f90

@@ -495,7 +495,7 @@ subroutine psb_d_bjac_precbld(a,desc_a,prec,info,amold,vmold,imold)
   use psb_d_bjacprec, psb_protect_name => psb_d_bjac_precbld
   Implicit None
 
-  type(psb_dspmat_type), intent(in), target :: a
+  type(psb_dspmat_type), intent(inout), target :: a
   type(psb_desc_type), intent(inout), target   :: desc_a
   class(psb_d_bjac_prec_type),intent(inout) :: prec
   integer(psb_ipk_), intent(out)                      :: info

prec/impl/psb_d_iluk_fact.f90

@@ -211,8 +211,7 @@ subroutine psb_diluk_fact(fill_in,ialg,a,l,u,d,info,blck,shft)
   call u%set_triangle()
   call u%set_unit()
   call u%set_upper()
-  write(0,*) 'In iluk__fact:',&
+
-       & l%is_triangle(),l%is_lower(),u%is_triangle(),u%is_upper()
   !
   ! Nullify pointer / deallocate memory
   !

prec/impl/psb_d_invk_fact.f90

@@ -41,7 +41,7 @@ subroutine psb_d_invk_bld(a,fill1, fill2,lmat,d,umat,desc,info,blck)
   implicit none
 
   ! Arguments
-  type(psb_dspmat_type), intent(in), target   :: a
+  type(psb_dspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fill1, fill2
   type(psb_dspmat_type), intent(inout)        :: lmat, umat
   real(psb_dpk_), allocatable                 :: d(:)
@@ -155,7 +155,7 @@ subroutine psb_dsparse_invk(n,a,z,fill_in,info,inlevs)
   use psb_d_invk_fact_mod, psb_protect_name => psb_dsparse_invk
 
   integer(psb_ipk_), intent(in)           :: n
-  type(psb_dspmat_type), intent(in)       :: a
+  type(psb_dspmat_type), intent(inout)   :: a
   type(psb_dspmat_type), intent(inout)    :: z
   integer(psb_ipk_), intent(in)           :: fill_in
   integer(psb_ipk_), intent(out)          :: info

prec/impl/psb_d_invt_fact.f90

@@ -73,7 +73,7 @@ subroutine psb_d_invt_bld(a,fillin,invfill,thresh,invthresh,&
   implicit none
 
   ! Arguments
-  type(psb_dspmat_type), intent(in), target   :: a
+  type(psb_dspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fillin,invfill
   real(psb_dpk_), intent(in)                :: thresh
   real(psb_dpk_), intent(in)                  :: invthresh
@@ -198,7 +198,7 @@ subroutine psb_dsparse_invt(n,a,z,nzrmax,sp_thresh,info)
 
   implicit none
   integer(psb_ipk_), intent(in)        :: n
-  type(psb_dspmat_type), intent(in)    :: a
+  type(psb_dspmat_type), intent(inout) :: a
   type(psb_dspmat_type), intent(inout) :: z
   integer(psb_ipk_), intent(in)        :: nzrmax
   real(psb_dpk_), intent(in)           :: sp_thresh

prec/impl/psb_s_bjacprec_impl.f90

@@ -495,7 +495,7 @@ subroutine psb_s_bjac_precbld(a,desc_a,prec,info,amold,vmold,imold)
   use psb_s_bjacprec, psb_protect_name => psb_s_bjac_precbld
   Implicit None
 
-  type(psb_sspmat_type), intent(in), target :: a
+  type(psb_sspmat_type), intent(inout), target :: a
   type(psb_desc_type), intent(inout), target   :: desc_a
   class(psb_s_bjac_prec_type),intent(inout) :: prec
   integer(psb_ipk_), intent(out)                      :: info

prec/impl/psb_s_invk_fact.f90

@@ -41,7 +41,7 @@ subroutine psb_s_invk_bld(a,fill1, fill2,lmat,d,umat,desc,info,blck)
   implicit none
 
   ! Arguments
-  type(psb_sspmat_type), intent(in), target   :: a
+  type(psb_sspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fill1, fill2
   type(psb_sspmat_type), intent(inout)        :: lmat, umat
   real(psb_spk_), allocatable                 :: d(:)
@@ -155,7 +155,7 @@ subroutine psb_ssparse_invk(n,a,z,fill_in,info,inlevs)
   use psb_s_invk_fact_mod, psb_protect_name => psb_ssparse_invk
 
   integer(psb_ipk_), intent(in)           :: n
-  type(psb_sspmat_type), intent(in)       :: a
+  type(psb_sspmat_type), intent(inout)   :: a
   type(psb_sspmat_type), intent(inout)    :: z
   integer(psb_ipk_), intent(in)           :: fill_in
   integer(psb_ipk_), intent(out)          :: info

prec/impl/psb_s_invt_fact.f90

@@ -73,7 +73,7 @@ subroutine psb_s_invt_bld(a,fillin,invfill,thresh,invthresh,&
   implicit none
 
   ! Arguments
-  type(psb_sspmat_type), intent(in), target   :: a
+  type(psb_sspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fillin,invfill
   real(psb_spk_), intent(in)                :: thresh
   real(psb_spk_), intent(in)                  :: invthresh
@@ -198,7 +198,7 @@ subroutine psb_ssparse_invt(n,a,z,nzrmax,sp_thresh,info)
 
   implicit none
   integer(psb_ipk_), intent(in)        :: n
-  type(psb_sspmat_type), intent(in)    :: a
+  type(psb_sspmat_type), intent(inout) :: a
   type(psb_sspmat_type), intent(inout) :: z
   integer(psb_ipk_), intent(in)        :: nzrmax
   real(psb_spk_), intent(in)           :: sp_thresh

prec/impl/psb_z_bjacprec_impl.f90

@@ -495,7 +495,7 @@ subroutine psb_z_bjac_precbld(a,desc_a,prec,info,amold,vmold,imold)
   use psb_z_bjacprec, psb_protect_name => psb_z_bjac_precbld
   Implicit None
 
-  type(psb_zspmat_type), intent(in), target :: a
+  type(psb_zspmat_type), intent(inout), target :: a
   type(psb_desc_type), intent(inout), target   :: desc_a
   class(psb_z_bjac_prec_type),intent(inout) :: prec
   integer(psb_ipk_), intent(out)                      :: info

prec/impl/psb_z_invk_fact.f90

@@ -41,7 +41,7 @@ subroutine psb_z_invk_bld(a,fill1, fill2,lmat,d,umat,desc,info,blck)
   implicit none
 
   ! Arguments
-  type(psb_zspmat_type), intent(in), target   :: a
+  type(psb_zspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fill1, fill2
   type(psb_zspmat_type), intent(inout)        :: lmat, umat
   complex(psb_dpk_), allocatable                 :: d(:)
@@ -155,7 +155,7 @@ subroutine psb_zsparse_invk(n,a,z,fill_in,info,inlevs)
   use psb_z_invk_fact_mod, psb_protect_name => psb_zsparse_invk
 
   integer(psb_ipk_), intent(in)           :: n
-  type(psb_zspmat_type), intent(in)       :: a
+  type(psb_zspmat_type), intent(inout)   :: a
   type(psb_zspmat_type), intent(inout)    :: z
   integer(psb_ipk_), intent(in)           :: fill_in
   integer(psb_ipk_), intent(out)          :: info

prec/impl/psb_z_invt_fact.f90

@@ -73,7 +73,7 @@ subroutine psb_z_invt_bld(a,fillin,invfill,thresh,invthresh,&
   implicit none
 
   ! Arguments
-  type(psb_zspmat_type), intent(in), target   :: a
+  type(psb_zspmat_type), intent(inout), target   :: a
   integer(psb_ipk_), intent(in)               :: fillin,invfill
   real(psb_dpk_), intent(in)                :: thresh
   real(psb_dpk_), intent(in)                  :: invthresh
@@ -198,7 +198,7 @@ subroutine psb_zsparse_invt(n,a,z,nzrmax,sp_thresh,info)
 
   implicit none
   integer(psb_ipk_), intent(in)        :: n
-  type(psb_zspmat_type), intent(in)    :: a
+  type(psb_zspmat_type), intent(inout) :: a
   type(psb_zspmat_type), intent(inout) :: z
   integer(psb_ipk_), intent(in)        :: nzrmax
   real(psb_dpk_), intent(in)           :: sp_thresh

prec/psb_c_invk_fact_mod.f90

@@ -81,7 +81,7 @@ module psb_c_invk_fact_mod
       ! import
       import psb_cspmat_type, psb_ipk_, psb_spk_, psb_desc_type
       ! Arguments
-      type(psb_cspmat_type), intent(in), target   :: a
+      type(psb_cspmat_type), intent(inout), target   :: a
       integer(psb_ipk_), intent(in)               :: fill1, fill2
       type(psb_cspmat_type), intent(inout)        :: lmat, umat
       complex(psb_spk_), allocatable                 :: d(:)
@@ -98,7 +98,7 @@ module psb_c_invk_fact_mod
       import psb_ipk_, psb_cspmat_type
       ! Arguments
       integer(psb_ipk_), intent(in)           :: n
-      type(psb_cspmat_type), intent(in)       :: a
+      type(psb_cspmat_type), intent(inout)  :: a
       type(psb_cspmat_type), intent(inout)    :: z
       integer(psb_ipk_), intent(in)           :: fill_in
       integer(psb_ipk_), intent(out)          :: info

prec/psb_c_invt_fact_mod.f90

@@ -82,7 +82,7 @@ module psb_c_invt_fact_mod
          ! Import
          import psb_cspmat_type, psb_spk_, psb_ipk_, psb_desc_type
          ! Arguments
-         type(psb_cspmat_type), intent(in), target   :: a
+         type(psb_cspmat_type), intent(inout), target   :: a
          integer(psb_ipk_), intent(in)               :: fillin,invfill
          real(psb_spk_), intent(in)                :: thresh
          real(psb_spk_), intent(in)                  :: invthresh
@@ -102,7 +102,7 @@ module psb_c_invt_fact_mod
          import psb_cspmat_type, psb_spk_, psb_ipk_
          ! Arguments
          integer(psb_ipk_), intent(in)        :: n
-         type(psb_cspmat_type), intent(in)    :: a
+         type(psb_cspmat_type), intent(inout) :: a
          type(psb_cspmat_type), intent(inout) :: z
          integer(psb_ipk_), intent(in)        :: nzrmax
          real(psb_spk_), intent(in)           :: sp_thresh

prec/psb_d_invk_fact_mod.f90

@@ -81,7 +81,7 @@ module psb_d_invk_fact_mod
       ! import
       import psb_dspmat_type, psb_ipk_, psb_dpk_, psb_desc_type
       ! Arguments
-      type(psb_dspmat_type), intent(in), target   :: a
+      type(psb_dspmat_type), intent(inout), target   :: a
       integer(psb_ipk_), intent(in)               :: fill1, fill2
       type(psb_dspmat_type), intent(inout)        :: lmat, umat
       real(psb_dpk_), allocatable                 :: d(:)
@@ -98,7 +98,7 @@ module psb_d_invk_fact_mod
       import psb_ipk_, psb_dspmat_type
       ! Arguments
       integer(psb_ipk_), intent(in)           :: n
-      type(psb_dspmat_type), intent(in)       :: a
+      type(psb_dspmat_type), intent(inout)  :: a
       type(psb_dspmat_type), intent(inout)    :: z
       integer(psb_ipk_), intent(in)           :: fill_in
       integer(psb_ipk_), intent(out)          :: info

prec/psb_d_invt_fact_mod.f90

@@ -82,7 +82,7 @@ module psb_d_invt_fact_mod
          ! Import
          import psb_dspmat_type, psb_dpk_, psb_ipk_, psb_desc_type
          ! Arguments
-         type(psb_dspmat_type), intent(in), target   :: a
+         type(psb_dspmat_type), intent(inout), target   :: a
          integer(psb_ipk_), intent(in)               :: fillin,invfill
          real(psb_dpk_), intent(in)                :: thresh
          real(psb_dpk_), intent(in)                  :: invthresh
@@ -102,7 +102,7 @@ module psb_d_invt_fact_mod
          import psb_dspmat_type, psb_dpk_, psb_ipk_
          ! Arguments
          integer(psb_ipk_), intent(in)        :: n
-         type(psb_dspmat_type), intent(in)    :: a
+         type(psb_dspmat_type), intent(inout) :: a
          type(psb_dspmat_type), intent(inout) :: z
          integer(psb_ipk_), intent(in)        :: nzrmax
          real(psb_dpk_), intent(in)           :: sp_thresh

prec/psb_s_invk_fact_mod.f90

@@ -81,7 +81,7 @@ module psb_s_invk_fact_mod
       ! import
       import psb_sspmat_type, psb_ipk_, psb_spk_, psb_desc_type
       ! Arguments
-      type(psb_sspmat_type), intent(in), target   :: a
+      type(psb_sspmat_type), intent(inout), target   :: a
       integer(psb_ipk_), intent(in)               :: fill1, fill2
       type(psb_sspmat_type), intent(inout)        :: lmat, umat
       real(psb_spk_), allocatable                 :: d(:)
@@ -98,7 +98,7 @@ module psb_s_invk_fact_mod
       import psb_ipk_, psb_sspmat_type
       ! Arguments
       integer(psb_ipk_), intent(in)           :: n
-      type(psb_sspmat_type), intent(in)       :: a
+      type(psb_sspmat_type), intent(inout)  :: a
       type(psb_sspmat_type), intent(inout)    :: z
       integer(psb_ipk_), intent(in)           :: fill_in
       integer(psb_ipk_), intent(out)          :: info

prec/psb_s_invt_fact_mod.f90

@@ -82,7 +82,7 @@ module psb_s_invt_fact_mod
          ! Import
          import psb_sspmat_type, psb_spk_, psb_ipk_, psb_desc_type
          ! Arguments
-         type(psb_sspmat_type), intent(in), target   :: a
+         type(psb_sspmat_type), intent(inout), target   :: a
          integer(psb_ipk_), intent(in)               :: fillin,invfill
          real(psb_spk_), intent(in)                :: thresh
          real(psb_spk_), intent(in)                  :: invthresh
@@ -102,7 +102,7 @@ module psb_s_invt_fact_mod
          import psb_sspmat_type, psb_spk_, psb_ipk_
          ! Arguments
          integer(psb_ipk_), intent(in)        :: n
-         type(psb_sspmat_type), intent(in)    :: a
+         type(psb_sspmat_type), intent(inout) :: a
          type(psb_sspmat_type), intent(inout) :: z
          integer(psb_ipk_), intent(in)        :: nzrmax
          real(psb_spk_), intent(in)           :: sp_thresh

prec/psb_z_invk_fact_mod.f90

@@ -81,7 +81,7 @@ module psb_z_invk_fact_mod
       ! import
       import psb_zspmat_type, psb_ipk_, psb_dpk_, psb_desc_type
       ! Arguments
-      type(psb_zspmat_type), intent(in), target   :: a
+      type(psb_zspmat_type), intent(inout), target   :: a
       integer(psb_ipk_), intent(in)               :: fill1, fill2
       type(psb_zspmat_type), intent(inout)        :: lmat, umat
       complex(psb_dpk_), allocatable                 :: d(:)
@@ -98,7 +98,7 @@ module psb_z_invk_fact_mod
       import psb_ipk_, psb_zspmat_type
       ! Arguments
       integer(psb_ipk_), intent(in)           :: n
-      type(psb_zspmat_type), intent(in)       :: a
+      type(psb_zspmat_type), intent(inout)  :: a
       type(psb_zspmat_type), intent(inout)    :: z
       integer(psb_ipk_), intent(in)           :: fill_in
       integer(psb_ipk_), intent(out)          :: info

prec/psb_z_invt_fact_mod.f90

@@ -82,7 +82,7 @@ module psb_z_invt_fact_mod
          ! Import
          import psb_zspmat_type, psb_dpk_, psb_ipk_, psb_desc_type
          ! Arguments
-         type(psb_zspmat_type), intent(in), target   :: a
+         type(psb_zspmat_type), intent(inout), target   :: a
          integer(psb_ipk_), intent(in)               :: fillin,invfill
          real(psb_dpk_), intent(in)                :: thresh
          real(psb_dpk_), intent(in)                  :: invthresh
@@ -102,7 +102,7 @@ module psb_z_invt_fact_mod
          import psb_zspmat_type, psb_dpk_, psb_ipk_
          ! Arguments
          integer(psb_ipk_), intent(in)        :: n
-         type(psb_zspmat_type), intent(in)    :: a
+         type(psb_zspmat_type), intent(inout) :: a
          type(psb_zspmat_type), intent(inout) :: z
          integer(psb_ipk_), intent(in)        :: nzrmax
          real(psb_dpk_), intent(in)           :: sp_thresh