Fixed intent(in) on X(:) for preconditioning.

psblas3-type-indexed
Salvatore Filippone 17 years ago
parent f46a8b1131
commit ee61f08f55

@ -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 use psb_prec_mod, psb_protect_name => psb_dbjac_aply
implicit none 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 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(in) :: alpha,beta real(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans real(kind(0.d0)),intent(in) :: alpha,beta
real(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info real(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
! Local variables ! Local variables
integer :: n_row,n_col integer :: n_row,n_col

@ -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 use psb_prec_mod, psb_protect_name => psb_dgprec_aply
implicit none 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 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(in) :: alpha,beta real(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans real(kind(0.d0)),intent(in) :: alpha,beta
real(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info real(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
! Local variables ! Local variables
integer :: n_row,int_err(5) integer :: n_row,int_err(5)

@ -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_desc_type),intent(in) :: desc_data
type(psb_dprec_type), intent(in) :: prec 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 integer, intent(out) :: info
character(len=1), optional :: trans character(len=1), optional :: trans
real(kind(0.d0)), optional, target :: work(:) real(kind(0.d0)), optional, target :: work(:)

@ -118,7 +118,8 @@ module psb_prec_mod
use psb_prec_type use psb_prec_type
type(psb_desc_type),intent(in) :: desc_data type(psb_desc_type),intent(in) :: desc_data
type(psb_dprec_type), intent(in) :: prec 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 integer, intent(out) :: info
character(len=1), optional :: trans character(len=1), optional :: trans
real(kind(0.d0)),intent(inout), optional, target :: work(:) real(kind(0.d0)),intent(inout), optional, target :: work(:)
@ -137,7 +138,8 @@ module psb_prec_mod
use psb_prec_type 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 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 integer, intent(out) :: info
character(len=1), optional :: trans character(len=1), optional :: trans
complex(kind(0.d0)),intent(inout), optional, target :: work(:) 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) subroutine psb_dbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
use psb_base_mod use psb_base_mod
use psb_prec_type use psb_prec_type
type(psb_desc_type), intent(in) :: desc_data type(psb_desc_type), intent(in) :: desc_data
type(psb_dprec_type), intent(in) :: prec 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(in) :: alpha,beta real(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans real(kind(0.d0)),intent(in) :: alpha,beta
real(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info real(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
end subroutine psb_dbjac_aply end subroutine psb_dbjac_aply
subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
use psb_base_mod use psb_base_mod
use psb_prec_type 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 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(in) :: alpha,beta complex(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans complex(kind(0.d0)),intent(in) :: alpha,beta
complex(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info complex(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
end subroutine psb_zbjac_aply end subroutine psb_zbjac_aply
end interface end interface
@ -246,7 +250,8 @@ module psb_prec_mod
use psb_prec_type use psb_prec_type
type(psb_desc_type),intent(in) :: desc_data type(psb_desc_type),intent(in) :: desc_data
type(psb_dprec_type), intent(in) :: prec 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 real(kind(0.d0)),intent(in) :: alpha,beta
character(len=1) :: trans character(len=1) :: trans
real(kind(0.d0)),target :: work(:) 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) subroutine psb_zgprec_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
use psb_base_mod use psb_base_mod
use psb_prec_type 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 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(in) :: alpha,beta complex(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans complex(kind(0.d0)),intent(in) :: alpha,beta
complex(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info complex(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
end subroutine psb_zgprec_aply end subroutine psb_zgprec_aply
end interface end interface

@ -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 use psb_prec_mod, psb_protect_name => psb_zbjac_aply
implicit none 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 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(in) :: alpha,beta complex(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans complex(kind(0.d0)),intent(in) :: alpha,beta
complex(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info complex(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
! Local variables ! Local variables
integer :: n_row,n_col integer :: n_row,n_col

@ -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 use psb_prec_mod, psb_protect_name => psb_zgprec_aply
implicit none 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 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(in) :: alpha,beta complex(kind(0.d0)),intent(inout) :: y(:)
character(len=1) :: trans complex(kind(0.d0)),intent(in) :: alpha,beta
complex(kind(0.d0)),target :: work(:) character(len=1) :: trans
integer, intent(out) :: info complex(kind(0.d0)),target :: work(:)
integer, intent(out) :: info
! Local variables ! Local variables
integer :: n_row,int_err(5) integer :: n_row,int_err(5)

@ -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_desc_type),intent(in) :: desc_data
type(psb_zprec_type), intent(in) :: prec 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 integer, intent(out) :: info
character(len=1), optional :: trans character(len=1), optional :: trans
complex(kind(0.d0)), optional, target :: work(:) complex(kind(0.d0)), optional, target :: work(:)

Loading…
Cancel
Save