prec/impl/psb_c_bjacprec_impl.f90
 prec/impl/psb_c_prec_type_impl.f90
 prec/impl/psb_cprecbld.f90
 prec/impl/psb_d_bjacprec_impl.f90
 prec/impl/psb_d_prec_type_impl.f90
 prec/impl/psb_dprecbld.f90
 prec/impl/psb_s_bjacprec_impl.f90
 prec/impl/psb_s_prec_type_impl.f90
 prec/impl/psb_sprecbld.f90
 prec/impl/psb_z_bjacprec_impl.f90
 prec/impl/psb_z_prec_type_impl.f90
 prec/impl/psb_zprecbld.f90
 prec/psb_c_base_prec_mod.f90
 prec/psb_d_base_prec_mod.f90
 prec/psb_s_base_prec_mod.f90
 prec/psb_z_base_prec_mod.f90

Some changes were overdone. Fixed.
psblas3-final
Salvatore Filippone 12 years ago
parent 2bf2843f7a
commit 45c37dbeb0

@ -38,7 +38,7 @@ subroutine psb_c_bjac_dump(prec,info,prefix,head)
integer(psb_ipk_), intent(out) :: info
character(len=*), intent(in), optional :: prefix,head
integer(psb_ipk_) :: i, j, il1, iln, lname, lev
integer(psb_mpik_) :: ictxt,iam, np
integer(psb_ipk_) :: ictxt,iam, np
character(len=80) :: prefix_
character(len=120) :: fname ! len should be at least 20 more than
@ -87,7 +87,7 @@ subroutine psb_c_bjac_apply_vect(alpha,prec,x,beta,y,desc_data,info,trans,work)
integer(psb_ipk_) :: n_row,n_col
complex(psb_spk_), pointer :: ww(:), aux(:)
type(psb_c_vect_type) :: wv, wv1
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -244,7 +244,7 @@ subroutine psb_c_bjac_apply(alpha,prec,x,beta,y,desc_data,info,trans,work)
! Local variables
integer(psb_ipk_) :: n_row,n_col
complex(psb_spk_), pointer :: ww(:), aux(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -444,7 +444,7 @@ subroutine psb_c_bjac_precbld(a,desc_a,prec,info,upd,amold,afmt,vmold)
type(psb_c_csr_sparse_mat), allocatable :: lf, uf
complex(psb_spk_), allocatable :: dd(:)
integer(psb_ipk_) :: nztota, err_act, n_row, nrow_a,n_col, nhalo
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
character(len=20) :: name='c_bjac_precbld'
character(len=20) :: ch_err

@ -76,7 +76,7 @@ subroutine psb_c_apply2_vect(prec,x,y,desc_data,info,trans,work)
character :: trans_
complex(psb_spk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -151,7 +151,7 @@ subroutine psb_c_apply1_vect(prec,x,desc_data,info,trans,work)
type(psb_c_vect_type) :: ww
character :: trans_
complex(psb_spk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -228,7 +228,7 @@ subroutine psb_c_apply2v(prec,x,y,desc_data,info,trans,work)
character :: trans_
complex(psb_spk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -297,7 +297,7 @@ subroutine psb_c_apply1v(prec,x,desc_data,info,trans)
character(len=1), optional :: trans
character :: trans_
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
complex(psb_spk_), pointer :: WW(:), w1(:)
character(len=20) :: name

@ -45,7 +45,7 @@ subroutine psb_cprecbld(a,desc_a,p,info,upd,amold,afmt,vmold)
class(psb_c_base_vect_type), intent(in), optional :: vmold
! Local scalars
integer(psb_mpik_) :: ictxt, me,np
integer(psb_ipk_) :: ictxt, me,np
integer(psb_ipk_) :: err, n_row, n_col,mglob, err_act
integer(psb_ipk_) :: int_err(5)
character :: upd_

@ -38,7 +38,7 @@ subroutine psb_d_bjac_dump(prec,info,prefix,head)
integer(psb_ipk_), intent(out) :: info
character(len=*), intent(in), optional :: prefix,head
integer(psb_ipk_) :: i, j, il1, iln, lname, lev
integer(psb_mpik_) :: ictxt,iam, np
integer(psb_ipk_) :: ictxt,iam, np
character(len=80) :: prefix_
character(len=120) :: fname ! len should be at least 20 more than
@ -87,7 +87,7 @@ subroutine psb_d_bjac_apply_vect(alpha,prec,x,beta,y,desc_data,info,trans,work)
integer(psb_ipk_) :: n_row,n_col
real(psb_dpk_), pointer :: ww(:), aux(:)
type(psb_d_vect_type) :: wv, wv1
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -244,7 +244,7 @@ subroutine psb_d_bjac_apply(alpha,prec,x,beta,y,desc_data,info,trans,work)
! Local variables
integer(psb_ipk_) :: n_row,n_col
real(psb_dpk_), pointer :: ww(:), aux(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -444,7 +444,7 @@ subroutine psb_d_bjac_precbld(a,desc_a,prec,info,upd,amold,afmt,vmold)
type(psb_d_csr_sparse_mat), allocatable :: lf, uf
real(psb_dpk_), allocatable :: dd(:)
integer(psb_ipk_) :: nztota, err_act, n_row, nrow_a,n_col, nhalo
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
character(len=20) :: name='d_bjac_precbld'
character(len=20) :: ch_err

@ -76,7 +76,7 @@ subroutine psb_d_apply2_vect(prec,x,y,desc_data,info,trans,work)
character :: trans_
real(psb_dpk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -151,7 +151,7 @@ subroutine psb_d_apply1_vect(prec,x,desc_data,info,trans,work)
type(psb_d_vect_type) :: ww
character :: trans_
real(psb_dpk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -228,7 +228,7 @@ subroutine psb_d_apply2v(prec,x,y,desc_data,info,trans,work)
character :: trans_
real(psb_dpk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -297,7 +297,7 @@ subroutine psb_d_apply1v(prec,x,desc_data,info,trans)
character(len=1), optional :: trans
character :: trans_
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
real(psb_dpk_), pointer :: WW(:), w1(:)
character(len=20) :: name

@ -45,7 +45,7 @@ subroutine psb_dprecbld(a,desc_a,p,info,upd,amold,afmt,vmold)
class(psb_d_base_vect_type), intent(in), optional :: vmold
! Local scalars
integer(psb_mpik_) :: ictxt, me,np
integer(psb_ipk_) :: ictxt, me,np
integer(psb_ipk_) :: err, n_row, n_col,mglob, err_act
integer(psb_ipk_) :: int_err(5)
character :: upd_

@ -38,7 +38,7 @@ subroutine psb_s_bjac_dump(prec,info,prefix,head)
integer(psb_ipk_), intent(out) :: info
character(len=*), intent(in), optional :: prefix,head
integer(psb_ipk_) :: i, j, il1, iln, lname, lev
integer(psb_mpik_) :: ictxt,iam, np
integer(psb_ipk_) :: ictxt,iam, np
character(len=80) :: prefix_
character(len=120) :: fname ! len should be at least 20 more than
@ -87,7 +87,7 @@ subroutine psb_s_bjac_apply_vect(alpha,prec,x,beta,y,desc_data,info,trans,work)
integer(psb_ipk_) :: n_row,n_col
real(psb_spk_), pointer :: ww(:), aux(:)
type(psb_s_vect_type) :: wv, wv1
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -244,7 +244,7 @@ subroutine psb_s_bjac_apply(alpha,prec,x,beta,y,desc_data,info,trans,work)
! Local variables
integer(psb_ipk_) :: n_row,n_col
real(psb_spk_), pointer :: ww(:), aux(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -444,7 +444,7 @@ subroutine psb_s_bjac_precbld(a,desc_a,prec,info,upd,amold,afmt,vmold)
type(psb_s_csr_sparse_mat), allocatable :: lf, uf
real(psb_spk_), allocatable :: dd(:)
integer(psb_ipk_) :: nztota, err_act, n_row, nrow_a,n_col, nhalo
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
character(len=20) :: name='s_bjac_precbld'
character(len=20) :: ch_err

@ -76,7 +76,7 @@ subroutine psb_s_apply2_vect(prec,x,y,desc_data,info,trans,work)
character :: trans_
real(psb_spk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -151,7 +151,7 @@ subroutine psb_s_apply1_vect(prec,x,desc_data,info,trans,work)
type(psb_s_vect_type) :: ww
character :: trans_
real(psb_spk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -228,7 +228,7 @@ subroutine psb_s_apply2v(prec,x,y,desc_data,info,trans,work)
character :: trans_
real(psb_spk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -297,7 +297,7 @@ subroutine psb_s_apply1v(prec,x,desc_data,info,trans)
character(len=1), optional :: trans
character :: trans_
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
real(psb_spk_), pointer :: WW(:), w1(:)
character(len=20) :: name

@ -45,7 +45,7 @@ subroutine psb_sprecbld(a,desc_a,p,info,upd,amold,afmt,vmold)
class(psb_s_base_vect_type), intent(in), optional :: vmold
! Local scalars
integer(psb_mpik_) :: ictxt, me,np
integer(psb_ipk_) :: ictxt, me,np
integer(psb_ipk_) :: err, n_row, n_col,mglob, err_act
integer(psb_ipk_) :: int_err(5)
character :: upd_

@ -38,7 +38,7 @@ subroutine psb_z_bjac_dump(prec,info,prefix,head)
integer(psb_ipk_), intent(out) :: info
character(len=*), intent(in), optional :: prefix,head
integer(psb_ipk_) :: i, j, il1, iln, lname, lev
integer(psb_mpik_) :: ictxt,iam, np
integer(psb_ipk_) :: ictxt,iam, np
character(len=80) :: prefix_
character(len=120) :: fname ! len should be at least 20 more than
@ -87,7 +87,7 @@ subroutine psb_z_bjac_apply_vect(alpha,prec,x,beta,y,desc_data,info,trans,work)
integer(psb_ipk_) :: n_row,n_col
complex(psb_dpk_), pointer :: ww(:), aux(:)
type(psb_z_vect_type) :: wv, wv1
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -244,7 +244,7 @@ subroutine psb_z_bjac_apply(alpha,prec,x,beta,y,desc_data,info,trans,work)
! Local variables
integer(psb_ipk_) :: n_row,n_col
complex(psb_dpk_), pointer :: ww(:), aux(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act, ierr(5)
integer(psb_ipk_) :: debug_level, debug_unit
character :: trans_
@ -444,7 +444,7 @@ subroutine psb_z_bjac_precbld(a,desc_a,prec,info,upd,amold,afmt,vmold)
type(psb_z_csr_sparse_mat), allocatable :: lf, uf
complex(psb_dpk_), allocatable :: dd(:)
integer(psb_ipk_) :: nztota, err_act, n_row, nrow_a,n_col, nhalo
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
character(len=20) :: name='z_bjac_precbld'
character(len=20) :: ch_err

@ -76,7 +76,7 @@ subroutine psb_z_apply2_vect(prec,x,y,desc_data,info,trans,work)
character :: trans_
complex(psb_dpk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -151,7 +151,7 @@ subroutine psb_z_apply1_vect(prec,x,desc_data,info,trans,work)
type(psb_z_vect_type) :: ww
character :: trans_
complex(psb_dpk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -228,7 +228,7 @@ subroutine psb_z_apply2v(prec,x,y,desc_data,info,trans,work)
character :: trans_
complex(psb_dpk_), pointer :: work_(:)
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
character(len=20) :: name
@ -297,7 +297,7 @@ subroutine psb_z_apply1v(prec,x,desc_data,info,trans)
character(len=1), optional :: trans
character :: trans_
integer(psb_mpik_) :: ictxt,np,me
integer(psb_ipk_) :: ictxt,np,me
integer(psb_ipk_) :: err_act
complex(psb_dpk_), pointer :: WW(:), w1(:)
character(len=20) :: name

@ -45,7 +45,7 @@ subroutine psb_zprecbld(a,desc_a,p,info,upd,amold,afmt,vmold)
class(psb_z_base_vect_type), intent(in), optional :: vmold
! Local scalars
integer(psb_mpik_) :: ictxt, me,np
integer(psb_ipk_) :: ictxt, me,np
integer(psb_ipk_) :: err, n_row, n_col,mglob, err_act
integer(psb_ipk_) :: int_err(5)
character :: upd_

@ -36,7 +36,7 @@
module psb_c_base_prec_mod
! Reduces size of .mod file.
use psb_base_mod, only : psb_spk_, psb_ipk_, psb_long_int_k_, psb_mpik_,&
use psb_base_mod, only : psb_spk_, psb_ipk_, psb_long_int_k_,&
& psb_desc_type, psb_sizeof, psb_free, psb_cdfree, psb_errpush, psb_act_abort_,&
& psb_sizeof_int, psb_sizeof_long_int, psb_sizeof_sp, psb_sizeof_dp, &
& psb_erractionsave, psb_erractionrestore, psb_error, psb_get_errstatus, psb_success_,&
@ -46,7 +46,7 @@ module psb_c_base_prec_mod
use psb_prec_const_mod
type psb_c_base_prec_type
integer(psb_mpik_) :: ictxt
integer(psb_ipk_) :: ictxt
contains
procedure, pass(prec) :: set_ctxt => psb_c_base_set_ctxt
procedure, pass(prec) :: get_ctxt => psb_c_base_get_ctxt
@ -407,7 +407,7 @@ contains
subroutine psb_c_base_set_ctxt(prec,ictxt)
implicit none
class(psb_c_base_prec_type), intent(inout) :: prec
integer(psb_mpik_), intent(in) :: ictxt
integer(psb_ipk_), intent(in) :: ictxt
prec%ictxt = ictxt
@ -423,7 +423,7 @@ contains
function psb_c_base_get_ctxt(prec) result(val)
class(psb_c_base_prec_type), intent(in) :: prec
integer(psb_mpik_) :: val
integer(psb_ipk_) :: val
val = prec%ictxt
return

@ -36,7 +36,7 @@
module psb_d_base_prec_mod
! Reduces size of .mod file.
use psb_base_mod, only : psb_dpk_, psb_ipk_, psb_long_int_k_, psb_mpik_,&
use psb_base_mod, only : psb_dpk_, psb_ipk_, psb_long_int_k_,&
& psb_desc_type, psb_sizeof, psb_free, psb_cdfree, psb_errpush, psb_act_abort_,&
& psb_sizeof_int, psb_sizeof_long_int, psb_sizeof_sp, psb_sizeof_dp, &
& psb_erractionsave, psb_erractionrestore, psb_error, psb_get_errstatus, psb_success_,&
@ -46,7 +46,7 @@ module psb_d_base_prec_mod
use psb_prec_const_mod
type psb_d_base_prec_type
integer(psb_mpik_) :: ictxt
integer(psb_ipk_) :: ictxt
contains
procedure, pass(prec) :: set_ctxt => psb_d_base_set_ctxt
procedure, pass(prec) :: get_ctxt => psb_d_base_get_ctxt
@ -407,7 +407,7 @@ contains
subroutine psb_d_base_set_ctxt(prec,ictxt)
implicit none
class(psb_d_base_prec_type), intent(inout) :: prec
integer(psb_mpik_), intent(in) :: ictxt
integer(psb_ipk_), intent(in) :: ictxt
prec%ictxt = ictxt
@ -423,7 +423,7 @@ contains
function psb_d_base_get_ctxt(prec) result(val)
class(psb_d_base_prec_type), intent(in) :: prec
integer(psb_mpik_) :: val
integer(psb_ipk_) :: val
val = prec%ictxt
return

@ -36,7 +36,7 @@
module psb_s_base_prec_mod
! Reduces size of .mod file.
use psb_base_mod, only : psb_spk_, psb_ipk_, psb_long_int_k_, psb_mpik_,&
use psb_base_mod, only : psb_spk_, psb_ipk_, psb_long_int_k_,&
& psb_desc_type, psb_sizeof, psb_free, psb_cdfree, psb_errpush, psb_act_abort_,&
& psb_sizeof_int, psb_sizeof_long_int, psb_sizeof_sp, psb_sizeof_dp, &
& psb_erractionsave, psb_erractionrestore, psb_error, psb_get_errstatus, psb_success_,&
@ -46,7 +46,7 @@ module psb_s_base_prec_mod
use psb_prec_const_mod
type psb_s_base_prec_type
integer(psb_mpik_) :: ictxt
integer(psb_ipk_) :: ictxt
contains
procedure, pass(prec) :: set_ctxt => psb_s_base_set_ctxt
procedure, pass(prec) :: get_ctxt => psb_s_base_get_ctxt
@ -407,7 +407,7 @@ contains
subroutine psb_s_base_set_ctxt(prec,ictxt)
implicit none
class(psb_s_base_prec_type), intent(inout) :: prec
integer(psb_mpik_), intent(in) :: ictxt
integer(psb_ipk_), intent(in) :: ictxt
prec%ictxt = ictxt
@ -423,7 +423,7 @@ contains
function psb_s_base_get_ctxt(prec) result(val)
class(psb_s_base_prec_type), intent(in) :: prec
integer(psb_mpik_) :: val
integer(psb_ipk_) :: val
val = prec%ictxt
return

@ -36,7 +36,7 @@
module psb_z_base_prec_mod
! Reduces size of .mod file.
use psb_base_mod, only : psb_dpk_, psb_ipk_, psb_long_int_k_, psb_mpik_,&
use psb_base_mod, only : psb_dpk_, psb_ipk_, psb_long_int_k_,&
& psb_desc_type, psb_sizeof, psb_free, psb_cdfree, psb_errpush, psb_act_abort_,&
& psb_sizeof_int, psb_sizeof_long_int, psb_sizeof_sp, psb_sizeof_dp, &
& psb_erractionsave, psb_erractionrestore, psb_error, psb_get_errstatus, psb_success_,&
@ -46,7 +46,7 @@ module psb_z_base_prec_mod
use psb_prec_const_mod
type psb_z_base_prec_type
integer(psb_mpik_) :: ictxt
integer(psb_ipk_) :: ictxt
contains
procedure, pass(prec) :: set_ctxt => psb_z_base_set_ctxt
procedure, pass(prec) :: get_ctxt => psb_z_base_get_ctxt
@ -407,7 +407,7 @@ contains
subroutine psb_z_base_set_ctxt(prec,ictxt)
implicit none
class(psb_z_base_prec_type), intent(inout) :: prec
integer(psb_mpik_), intent(in) :: ictxt
integer(psb_ipk_), intent(in) :: ictxt
prec%ictxt = ictxt
@ -423,7 +423,7 @@ contains
function psb_z_base_get_ctxt(prec) result(val)
class(psb_z_base_prec_type), intent(in) :: prec
integer(psb_mpik_) :: val
integer(psb_ipk_) :: val
val = prec%ictxt
return

Loading…
Cancel
Save