Fixed comment on trans=N: we now handle trans=T correctly in multiple

Jacobi sweeps.
stopcriterion
Salvatore Filippone 17 years ago
parent 38ff0f0373
commit 1a00c98ae8

@ -214,7 +214,6 @@ subroutine mld_dsub_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! Apply prec%iprcparm(smooth_sweeps_) sweeps of a block-Jacobi solver ! Apply prec%iprcparm(smooth_sweeps_) sweeps of a block-Jacobi solver
! to compute an approximate solution of a linear system. ! to compute an approximate solution of a linear system.
! !
! Note: trans is always 'N' here.
! !
if (size(prec%av) < mld_ap_nd_) then if (size(prec%av) < mld_ap_nd_) then

@ -215,8 +215,6 @@ subroutine mld_zsub_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
! Apply prec%iprcparm(smooth_sweeps_) sweeps of a block-Jacobi solver ! Apply prec%iprcparm(smooth_sweeps_) sweeps of a block-Jacobi solver
! to compute an approximate solution of a linear system. ! to compute an approximate solution of a linear system.
! !
! Note: trans is always 'N' here.
!
if (size(prec%av) < mld_ap_nd_) then if (size(prec%av) < mld_ap_nd_) then
info = 4011 info = 4011

Loading…
Cancel
Save