diff --git a/prec/impl/psb_c_bjacprec_impl.f90 b/prec/impl/psb_c_bjacprec_impl.f90 index 32baf385..74948530 100644 --- a/prec/impl/psb_c_bjacprec_impl.f90 +++ b/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 diff --git a/prec/impl/psb_c_invk_fact.f90 b/prec/impl/psb_c_invk_fact.f90 index 5c047b18..269070e7 100644 --- a/prec/impl/psb_c_invk_fact.f90 +++ b/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 diff --git a/prec/impl/psb_c_invt_fact.f90 b/prec/impl/psb_c_invt_fact.f90 index 73e242f6..f5b174f9 100644 --- a/prec/impl/psb_c_invt_fact.f90 +++ b/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 diff --git a/prec/impl/psb_d_bjacprec_impl.f90 b/prec/impl/psb_d_bjacprec_impl.f90 index ec866dbe..7142510a 100644 --- a/prec/impl/psb_d_bjacprec_impl.f90 +++ b/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 diff --git a/prec/impl/psb_d_iluk_fact.f90 b/prec/impl/psb_d_iluk_fact.f90 index 6089db7c..dda5b300 100644 --- a/prec/impl/psb_d_iluk_fact.f90 +++ b/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 ! diff --git a/prec/impl/psb_d_invk_fact.f90 b/prec/impl/psb_d_invk_fact.f90 index f1a09985..3ea5868c 100644 --- a/prec/impl/psb_d_invk_fact.f90 +++ b/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 diff --git a/prec/impl/psb_d_invt_fact.f90 b/prec/impl/psb_d_invt_fact.f90 index 83aa73ff..bd25bf54 100644 --- a/prec/impl/psb_d_invt_fact.f90 +++ b/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 diff --git a/prec/impl/psb_s_bjacprec_impl.f90 b/prec/impl/psb_s_bjacprec_impl.f90 index d4a4fd17..4250cc47 100644 --- a/prec/impl/psb_s_bjacprec_impl.f90 +++ b/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 diff --git a/prec/impl/psb_s_invk_fact.f90 b/prec/impl/psb_s_invk_fact.f90 index 4e7bb755..e45dcaea 100644 --- a/prec/impl/psb_s_invk_fact.f90 +++ b/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 diff --git a/prec/impl/psb_s_invt_fact.f90 b/prec/impl/psb_s_invt_fact.f90 index f311860d..eb5dae7d 100644 --- a/prec/impl/psb_s_invt_fact.f90 +++ b/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 diff --git a/prec/impl/psb_z_bjacprec_impl.f90 b/prec/impl/psb_z_bjacprec_impl.f90 index 3533f1e3..fb5d0f90 100644 --- a/prec/impl/psb_z_bjacprec_impl.f90 +++ b/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 diff --git a/prec/impl/psb_z_invk_fact.f90 b/prec/impl/psb_z_invk_fact.f90 index 99489acd..e92f5158 100644 --- a/prec/impl/psb_z_invk_fact.f90 +++ b/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 diff --git a/prec/impl/psb_z_invt_fact.f90 b/prec/impl/psb_z_invt_fact.f90 index bed713ae..bce1d640 100644 --- a/prec/impl/psb_z_invt_fact.f90 +++ b/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 diff --git a/prec/psb_c_ilu_fact_mod.f90 b/prec/psb_c_ilu_fact_mod.f90 index dcaa4b4a..c5065ea6 100644 --- a/prec/psb_c_ilu_fact_mod.f90 +++ b/prec/psb_c_ilu_fact_mod.f90 @@ -79,7 +79,7 @@ module psb_c_ilu_fact_mod use psb_base_mod use psb_prec_const_mod - interface psb_ilu0_fact + interface psb_ilu0_fact subroutine psb_cilu0_fact(ialg,a,l,u,d,info,blck,upd,shft) import psb_cspmat_type, psb_spk_, psb_ipk_ integer(psb_ipk_), intent(in) :: ialg diff --git a/prec/psb_c_invk_fact_mod.f90 b/prec/psb_c_invk_fact_mod.f90 index 620a8adf..62391cdb 100644 --- a/prec/psb_c_invk_fact_mod.f90 +++ b/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 diff --git a/prec/psb_c_invt_fact_mod.f90 b/prec/psb_c_invt_fact_mod.f90 index 841c39b1..b5457e32 100644 --- a/prec/psb_c_invt_fact_mod.f90 +++ b/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 diff --git a/prec/psb_d_ilu_fact_mod.f90 b/prec/psb_d_ilu_fact_mod.f90 index 934cf784..4372bca4 100644 --- a/prec/psb_d_ilu_fact_mod.f90 +++ b/prec/psb_d_ilu_fact_mod.f90 @@ -79,7 +79,7 @@ module psb_d_ilu_fact_mod use psb_base_mod use psb_prec_const_mod - interface psb_ilu0_fact + interface psb_ilu0_fact subroutine psb_dilu0_fact(ialg,a,l,u,d,info,blck,upd,shft) import psb_dspmat_type, psb_dpk_, psb_ipk_ integer(psb_ipk_), intent(in) :: ialg diff --git a/prec/psb_d_invk_fact_mod.f90 b/prec/psb_d_invk_fact_mod.f90 index 2bd97198..c72eba58 100644 --- a/prec/psb_d_invk_fact_mod.f90 +++ b/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 diff --git a/prec/psb_d_invt_fact_mod.f90 b/prec/psb_d_invt_fact_mod.f90 index f38c1c2b..947340ba 100644 --- a/prec/psb_d_invt_fact_mod.f90 +++ b/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 diff --git a/prec/psb_s_ilu_fact_mod.f90 b/prec/psb_s_ilu_fact_mod.f90 index d1b558e6..c35ed5bb 100644 --- a/prec/psb_s_ilu_fact_mod.f90 +++ b/prec/psb_s_ilu_fact_mod.f90 @@ -79,7 +79,7 @@ module psb_s_ilu_fact_mod use psb_base_mod use psb_prec_const_mod - interface psb_ilu0_fact + interface psb_ilu0_fact subroutine psb_silu0_fact(ialg,a,l,u,d,info,blck,upd,shft) import psb_sspmat_type, psb_spk_, psb_ipk_ integer(psb_ipk_), intent(in) :: ialg diff --git a/prec/psb_s_invk_fact_mod.f90 b/prec/psb_s_invk_fact_mod.f90 index 6b0d3553..cfbc59b9 100644 --- a/prec/psb_s_invk_fact_mod.f90 +++ b/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 diff --git a/prec/psb_s_invt_fact_mod.f90 b/prec/psb_s_invt_fact_mod.f90 index 2c9ce38c..b3788f6f 100644 --- a/prec/psb_s_invt_fact_mod.f90 +++ b/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 diff --git a/prec/psb_z_ilu_fact_mod.f90 b/prec/psb_z_ilu_fact_mod.f90 index 394c24d7..829e7f4d 100644 --- a/prec/psb_z_ilu_fact_mod.f90 +++ b/prec/psb_z_ilu_fact_mod.f90 @@ -79,7 +79,7 @@ module psb_z_ilu_fact_mod use psb_base_mod use psb_prec_const_mod - interface psb_ilu0_fact + interface psb_ilu0_fact subroutine psb_zilu0_fact(ialg,a,l,u,d,info,blck,upd,shft) import psb_zspmat_type, psb_dpk_, psb_ipk_ integer(psb_ipk_), intent(in) :: ialg diff --git a/prec/psb_z_invk_fact_mod.f90 b/prec/psb_z_invk_fact_mod.f90 index 0a1e5faf..aafa497e 100644 --- a/prec/psb_z_invk_fact_mod.f90 +++ b/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 diff --git a/prec/psb_z_invt_fact_mod.f90 b/prec/psb_z_invt_fact_mod.f90 index 1cdf32f4..d3326d0d 100644 --- a/prec/psb_z_invt_fact_mod.f90 +++ b/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