fdurastante
/
amg4psblas
mirror of https://github.com/sfilippone/amg4psblas.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Browse Source
stopcriterion
Salvatore Filippone 17 years ago
parent 7204231d15
commit f1c5dbcb3d
2 changed files with 14 additions and 22 deletions
Show Stats Download Patch File Download Diff File

18
mlprec/mld_das_aply.f90
Unescape Escape View File

@ -94,6 +94,7 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
real(kind(1.d0)), pointer :: ww(:), aux(:), tx(:),ty(:) real(kind(1.d0)), pointer :: ww(:), aux(:), tx(:),ty(:)
integer :: ictxt,np,me,isz, err_act integer :: ictxt,np,me,isz, err_act
character(len=20) :: name, ch_err character(len=20) :: name, ch_err
character :: trans_
name='mld_das_aply' name='mld_das_aply'
info = 0 info = 0
@ -103,6 +104,7 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
call psb_info(ictxt, me, np) call psb_info(ictxt, me, np)
trans_ = toupper(trans)
select case(prec%iprcparm(mld_prec_type_)) select case(prec%iprcparm(mld_prec_type_))
@ -115,7 +117,7 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! !
! shortcut: this fixes performance for RAS(0) == BJA ! shortcut: this fixes performance for RAS(0) == BJA
! !
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='psb_bjacaply' ch_err='psb_bjacaply'
@ -124,10 +126,6 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
else else
!
! Note: currently trans is unused
!
n_row = psb_cd_get_local_rows(prec%desc_data) n_row = psb_cd_get_local_rows(prec%desc_data)
n_col = psb_cd_get_local_cols(prec%desc_data) n_col = psb_cd_get_local_cols(prec%desc_data)
nrow_d = psb_cd_get_local_rows(desc_data) nrow_d = psb_cd_get_local_rows(desc_data)
@ -169,9 +167,8 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
tx(1:nrow_d) = x(1:nrow_d) tx(1:nrow_d) = x(1:nrow_d)
tx(nrow_d+1:isz) = dzero tx(nrow_d+1:isz) = dzero
select case(toupper(trans)) select case(trans_)
case('N') case('N')
! !
! Get the overlap entries of tx (tx==x) ! Get the overlap entries of tx (tx==x)
! !
@ -205,7 +202,7 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! block-Jacobi solver can be applied at the coarsest level of a multilevel ! block-Jacobi solver can be applied at the coarsest level of a multilevel
! preconditioner). The resulting vector is ty. ! preconditioner). The resulting vector is ty.
! !
call mld_bjac_aply(done,prec,tx,dzero,ty,prec%desc_data,toupper(trans),aux,info) call mld_bjac_aply(done,prec,tx,dzero,ty,prec%desc_data,trans_,aux,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='mld_bjac_aply' ch_err='mld_bjac_aply'
@ -252,7 +249,6 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
case('T','C') case('T','C')
! !
! With transpose, we have to do it here ! With transpose, we have to do it here
! !
@ -318,7 +314,7 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! block-Jacobi solver can be applied at the coarsest level of a multilevel ! block-Jacobi solver can be applied at the coarsest level of a multilevel
! preconditioner). The resulting vector is ty. ! preconditioner). The resulting vector is ty.
! !
call mld_bjac_aply(done,prec,tx,dzero,ty,prec%desc_data,toupper(trans),aux,info) call mld_bjac_aply(done,prec,tx,dzero,ty,prec%desc_data,trans_,aux,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='mld_bjac_aply' ch_err='mld_bjac_aply'
@ -360,8 +356,6 @@ subroutine mld_das_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
goto 9999 goto 9999
end select end select
! !
! Compute y = beta*y + alpha*ty (ty==K^(-1)*tx) ! Compute y = beta*y + alpha*ty (ty==K^(-1)*tx)
! !

14
mlprec/mld_zas_aply.f90
Unescape Escape View File

@ -94,6 +94,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
complex(kind(1.d0)), pointer :: ww(:), aux(:), tx(:),ty(:) complex(kind(1.d0)), pointer :: ww(:), aux(:), tx(:),ty(:)
integer :: ictxt,np,me,isz, err_act integer :: ictxt,np,me,isz, err_act
character(len=20) :: name, ch_err character(len=20) :: name, ch_err
character :: trans_
name='mld_zas_aply' name='mld_zas_aply'
info = 0 info = 0
@ -103,6 +104,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
call psb_info(ictxt, me, np) call psb_info(ictxt, me, np)
trans_ = toupper(trans)
select case(prec%iprcparm(mld_prec_type_)) select case(prec%iprcparm(mld_prec_type_))
@ -115,7 +117,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! !
! shortcut: this fixes performance for RAS(0) == BJA ! shortcut: this fixes performance for RAS(0) == BJA
! !
call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) call mld_bjac_aply(alpha,prec,x,beta,y,desc_data,trans_,work,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='psb_bjacaply' ch_err='psb_bjacaply'
@ -124,10 +126,6 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
else else
!
! Note: currently trans is unused
!
n_row = psb_cd_get_local_rows(prec%desc_data) n_row = psb_cd_get_local_rows(prec%desc_data)
n_col = psb_cd_get_local_cols(prec%desc_data) n_col = psb_cd_get_local_cols(prec%desc_data)
nrow_d = psb_cd_get_local_rows(desc_data) nrow_d = psb_cd_get_local_rows(desc_data)
@ -169,7 +167,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
tx(1:nrow_d) = x(1:nrow_d) tx(1:nrow_d) = x(1:nrow_d)
tx(nrow_d+1:isz) = dzero tx(nrow_d+1:isz) = dzero
select case(toupper(trans)) select case(trans_)
case('N') case('N')
! !
! Get the overlap entries of tx (tx==x) ! Get the overlap entries of tx (tx==x)
@ -204,7 +202,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! block-Jacobi solver can be applied at the coarsest level of a multilevel ! block-Jacobi solver can be applied at the coarsest level of a multilevel
! preconditioner). The resulting vector is ty. ! preconditioner). The resulting vector is ty.
! !
call mld_bjac_aply(zone,prec,tx,zzero,ty,prec%desc_data,toupper(trans),aux,info) call mld_bjac_aply(zone,prec,tx,zzero,ty,prec%desc_data,trans_,aux,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='mld_bjac_aply' ch_err='mld_bjac_aply'
@ -316,7 +314,7 @@ subroutine mld_zas_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! block-Jacobi solver can be applied at the coarsest level of a multilevel ! block-Jacobi solver can be applied at the coarsest level of a multilevel
! preconditioner). The resulting vector is ty. ! preconditioner). The resulting vector is ty.
! !
call mld_bjac_aply(zone,prec,tx,zzero,ty,prec%desc_data,toupper(trans),aux,info) call mld_bjac_aply(zone,prec,tx,zzero,ty,prec%desc_data,trans_,aux,info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010
ch_err='mld_bjac_aply' ch_err='mld_bjac_aply'

Write Preview
Loading…
Cancel
Save