Renaming inner precaply.

psblas3-type-indexed
Salvatore Filippone 18 years ago
parent dff1d9c27d
commit 002adffb28

@ -28,7 +28,7 @@
!!$ POSSIBILITY OF SUCH DAMAGE. !!$ POSSIBILITY OF SUCH DAMAGE.
!!$ !!$
!!$ !!$
subroutine psb_dbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) subroutine psb_dgen_precaply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! !
! Compute Y <- beta*Y + alpha*K^-1 X ! Compute Y <- beta*Y + alpha*K^-1 X
! where K is a a basic preconditioner stored in prec ! where K is a a basic preconditioner stored in prec
@ -145,5 +145,5 @@ subroutine psb_dbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
end if end if
return return
end subroutine psb_dbaseprc_aply end subroutine psb_dgen_precaply

@ -48,18 +48,18 @@ subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans, work)
logical,parameter :: debug=.false., debugprt=.false. logical,parameter :: debug=.false., debugprt=.false.
character(len=20) :: name character(len=20) :: name
interface psb_baseprc_aply interface psb_gen_precaply
subroutine psb_dbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) subroutine psb_dgen_precaply(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(inout) :: x(:), 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(:)
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_dbaseprc_aply end subroutine psb_dgen_precaply
end interface end interface
name='psb_prc_aply' name='psb_prc_aply'
@ -86,7 +86,7 @@ subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans, work)
end if end if
call psb_baseprc_aply(done,prec,x,dzero,y,desc_data,trans_,work_,info) call psb_gen_precaply(done,prec,x,dzero,y,desc_data,trans_,work_,info)
if (present(work)) then if (present(work)) then
else else

@ -81,21 +81,6 @@ module psb_prec_mod
end interface end interface
!!$ interface psb_precfree
!!$ subroutine psb_dprecfree(p,info)
!!$ use psb_base_mod
!!$ use psb_prec_type
!!$ type(psb_dprec_type), intent(inout) :: p
!!$ integer, intent(out) :: info
!!$ end subroutine psb_dprecfree
!!$ subroutine psb_zprecfree(p,info)
!!$ use psb_base_mod
!!$ use psb_prec_type
!!$ type(psb_zprec_type), intent(inout) :: p
!!$ integer, intent(out) :: info
!!$ end subroutine psb_zprecfree
!!$ end interface
interface psb_precaply interface psb_precaply
subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans,work) subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans,work)
use psb_base_mod use psb_base_mod

@ -48,8 +48,8 @@ subroutine psb_zprc_aply(prec,x,y,desc_data,info,trans, work)
logical,parameter :: debug=.false., debugprt=.false. logical,parameter :: debug=.false., debugprt=.false.
character(len=20) :: name character(len=20) :: name
interface psb_baseprc_aply interface psb_gen_precaply
subroutine psb_zbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) subroutine psb_zgen_precaply(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
@ -59,7 +59,7 @@ subroutine psb_zprc_aply(prec,x,y,desc_data,info,trans, work)
character(len=1) :: trans character(len=1) :: trans
complex(kind(0.d0)),target :: work(:) complex(kind(0.d0)),target :: work(:)
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_zbaseprc_aply end subroutine psb_zgen_precaply
end interface end interface
name='psb_prc_aply' name='psb_prc_aply'
@ -86,7 +86,7 @@ subroutine psb_zprc_aply(prec,x,y,desc_data,info,trans, work)
end if end if
call psb_baseprc_aply(zone,prec,x,zzero,y,desc_data,trans_, work_,info) call psb_gen_precaply(zone,prec,x,zzero,y,desc_data,trans_, work_,info)
if (present(work)) then if (present(work)) then
else else

Loading…
Cancel
Save