From ee61f08f55af598f2d83554f2212a5c5d3474d9f Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Wed, 14 Nov 2007 09:55:59 +0000 Subject: [PATCH] Fixed intent(in) on X(:) for preconditioning. --- prec/psb_dbjac_aply.f90 | 13 +++++----- prec/psb_dgprec_aply.f90 | 13 +++++----- prec/psb_dprc_aply.f90 | 3 ++- prec/psb_prec_mod.f90 | 52 ++++++++++++++++++++++------------------ prec/psb_zbjac_aply.f90 | 13 +++++----- prec/psb_zgprec_aply.f90 | 15 ++++++------ prec/psb_zprc_aply.f90 | 3 ++- 7 files changed, 62 insertions(+), 50 deletions(-) diff --git a/prec/psb_dbjac_aply.f90 b/prec/psb_dbjac_aply.f90 index 01f30eea..023c5fa5 100644 --- a/prec/psb_dbjac_aply.f90 +++ b/prec/psb_dbjac_aply.f90 @@ -40,13 +40,14 @@ subroutine psb_dbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_prec_mod, psb_protect_name => psb_dbjac_aply implicit none - type(psb_desc_type), intent(in) :: desc_data + type(psb_desc_type), intent(in) :: desc_data type(psb_dprec_type), intent(in) :: prec - real(kind(0.d0)),intent(inout) :: x(:), y(:) - real(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - real(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + real(kind(0.d0)),intent(in) :: x(:) + real(kind(0.d0)),intent(inout) :: y(:) + real(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + real(kind(0.d0)),target :: work(:) + integer, intent(out) :: info ! Local variables integer :: n_row,n_col diff --git a/prec/psb_dgprec_aply.f90 b/prec/psb_dgprec_aply.f90 index 339d4a3a..1d3985fb 100644 --- a/prec/psb_dgprec_aply.f90 +++ b/prec/psb_dgprec_aply.f90 @@ -38,13 +38,14 @@ subroutine psb_dgprec_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_prec_mod, psb_protect_name => psb_dgprec_aply implicit none - type(psb_desc_type),intent(in) :: desc_data + type(psb_desc_type),intent(in) :: desc_data type(psb_dprec_type), intent(in) :: prec - real(kind(0.d0)),intent(inout) :: x(:), y(:) - real(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - real(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + real(kind(0.d0)),intent(in) :: x(:) + real(kind(0.d0)),intent(inout) :: y(:) + real(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + real(kind(0.d0)),target :: work(:) + integer, intent(out) :: info ! Local variables integer :: n_row,int_err(5) diff --git a/prec/psb_dprc_aply.f90 b/prec/psb_dprc_aply.f90 index 5b98d674..e813ec06 100644 --- a/prec/psb_dprc_aply.f90 +++ b/prec/psb_dprc_aply.f90 @@ -36,7 +36,8 @@ subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans, work) type(psb_desc_type),intent(in) :: desc_data type(psb_dprec_type), intent(in) :: prec - real(kind(0.d0)),intent(inout) :: x(:), y(:) + real(kind(0.d0)),intent(in) :: x(:) + real(kind(0.d0)),intent(inout) :: y(:) integer, intent(out) :: info character(len=1), optional :: trans real(kind(0.d0)), optional, target :: work(:) diff --git a/prec/psb_prec_mod.f90 b/prec/psb_prec_mod.f90 index fedc3797..702f9f32 100644 --- a/prec/psb_prec_mod.f90 +++ b/prec/psb_prec_mod.f90 @@ -118,7 +118,8 @@ module psb_prec_mod use psb_prec_type type(psb_desc_type),intent(in) :: desc_data type(psb_dprec_type), intent(in) :: prec - real(kind(0.d0)),intent(inout) :: x(:), y(:) + real(kind(0.d0)),intent(in) :: x(:) + real(kind(0.d0)),intent(inout) :: y(:) integer, intent(out) :: info character(len=1), optional :: trans real(kind(0.d0)),intent(inout), optional, target :: work(:) @@ -137,7 +138,8 @@ module psb_prec_mod use psb_prec_type type(psb_desc_type),intent(in) :: desc_data type(psb_zprec_type), intent(in) :: prec - complex(kind(0.d0)),intent(inout) :: x(:), y(:) + complex(kind(0.d0)),intent(in) :: x(:) + complex(kind(0.d0)),intent(inout) :: y(:) integer, intent(out) :: info character(len=1), optional :: trans complex(kind(0.d0)),intent(inout), optional, target :: work(:) @@ -158,24 +160,26 @@ module psb_prec_mod subroutine psb_dbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_base_mod use psb_prec_type - type(psb_desc_type), intent(in) :: desc_data - type(psb_dprec_type), intent(in) :: prec - real(kind(0.d0)),intent(inout) :: x(:), y(:) - real(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - real(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + type(psb_desc_type), intent(in) :: desc_data + type(psb_dprec_type), intent(in) :: prec + real(kind(0.d0)),intent(in) :: x(:) + real(kind(0.d0)),intent(inout) :: y(:) + real(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + real(kind(0.d0)),target :: work(:) + integer, intent(out) :: info end subroutine psb_dbjac_aply subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_base_mod use psb_prec_type - type(psb_desc_type), intent(in) :: desc_data + type(psb_desc_type), intent(in) :: desc_data type(psb_zprec_type), intent(in) :: prec - complex(kind(0.d0)),intent(inout) :: x(:), y(:) - complex(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - complex(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + complex(kind(0.d0)),intent(in) :: x(:) + complex(kind(0.d0)),intent(inout) :: y(:) + complex(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + complex(kind(0.d0)),target :: work(:) + integer, intent(out) :: info end subroutine psb_zbjac_aply end interface @@ -246,7 +250,8 @@ module psb_prec_mod use psb_prec_type type(psb_desc_type),intent(in) :: desc_data type(psb_dprec_type), intent(in) :: prec - real(kind(0.d0)),intent(inout) :: x(:), y(:) + real(kind(0.d0)),intent(in) :: x(:) + real(kind(0.d0)),intent(inout) :: y(:) real(kind(0.d0)),intent(in) :: alpha,beta character(len=1) :: trans real(kind(0.d0)),target :: work(:) @@ -256,13 +261,14 @@ module psb_prec_mod subroutine psb_zgprec_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_base_mod use psb_prec_type - type(psb_desc_type),intent(in) :: desc_data - type(psb_zprec_type), intent(in) :: prec - complex(kind(0.d0)),intent(inout) :: x(:), y(:) - complex(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - complex(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + type(psb_desc_type),intent(in) :: desc_data + type(psb_zprec_type), intent(in) :: prec + complex(kind(0.d0)),intent(in) :: x(:) + complex(kind(0.d0)),intent(inout) :: y(:) + complex(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + complex(kind(0.d0)),target :: work(:) + integer, intent(out) :: info end subroutine psb_zgprec_aply end interface diff --git a/prec/psb_zbjac_aply.f90 b/prec/psb_zbjac_aply.f90 index 9f39fca8..917febe5 100644 --- a/prec/psb_zbjac_aply.f90 +++ b/prec/psb_zbjac_aply.f90 @@ -40,13 +40,14 @@ subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_prec_mod, psb_protect_name => psb_zbjac_aply implicit none - type(psb_desc_type), intent(in) :: desc_data + type(psb_desc_type), intent(in) :: desc_data type(psb_zprec_type), intent(in) :: prec - complex(kind(0.d0)),intent(inout) :: x(:), y(:) - complex(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - complex(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + complex(kind(0.d0)),intent(in) :: x(:) + complex(kind(0.d0)),intent(inout) :: y(:) + complex(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + complex(kind(0.d0)),target :: work(:) + integer, intent(out) :: info ! Local variables integer :: n_row,n_col diff --git a/prec/psb_zgprec_aply.f90 b/prec/psb_zgprec_aply.f90 index 06639248..d02df12d 100644 --- a/prec/psb_zgprec_aply.f90 +++ b/prec/psb_zgprec_aply.f90 @@ -38,13 +38,14 @@ subroutine psb_zgprec_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) use psb_prec_mod, psb_protect_name => psb_zgprec_aply implicit none - type(psb_desc_type),intent(in) :: desc_data - type(psb_zprec_type), intent(in) :: prec - complex(kind(0.d0)),intent(inout) :: x(:), y(:) - complex(kind(0.d0)),intent(in) :: alpha,beta - character(len=1) :: trans - complex(kind(0.d0)),target :: work(:) - integer, intent(out) :: info + type(psb_desc_type),intent(in) :: desc_data + type(psb_zprec_type), intent(in) :: prec + complex(kind(0.d0)),intent(in) :: x(:) + complex(kind(0.d0)),intent(inout) :: y(:) + complex(kind(0.d0)),intent(in) :: alpha,beta + character(len=1) :: trans + complex(kind(0.d0)),target :: work(:) + integer, intent(out) :: info ! Local variables integer :: n_row,int_err(5) diff --git a/prec/psb_zprc_aply.f90 b/prec/psb_zprc_aply.f90 index 77ab31e9..3e4734ca 100644 --- a/prec/psb_zprc_aply.f90 +++ b/prec/psb_zprc_aply.f90 @@ -37,7 +37,8 @@ subroutine psb_zprc_aply(prec,x,y,desc_data,info,trans, work) type(psb_desc_type),intent(in) :: desc_data type(psb_zprec_type), intent(in) :: prec - complex(kind(0.d0)),intent(inout) :: x(:), y(:) + complex(kind(0.d0)),intent(in) :: x(:) + complex(kind(0.d0)),intent(inout) :: y(:) integer, intent(out) :: info character(len=1), optional :: trans complex(kind(0.d0)), optional, target :: work(:)