From 1a00c98ae8e27394fd4e6ac2b193cb40327de6cd Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Mon, 21 Jan 2008 10:11:28 +0000 Subject: [PATCH] Fixed comment on trans=N: we now handle trans=T correctly in multiple Jacobi sweeps. --- mlprec/mld_dsub_aply.f90 | 1 - mlprec/mld_zsub_aply.f90 | 2 -- 2 files changed, 3 deletions(-) diff --git a/mlprec/mld_dsub_aply.f90 b/mlprec/mld_dsub_aply.f90 index a17880f8..f4c6c8cc 100644 --- a/mlprec/mld_dsub_aply.f90 +++ b/mlprec/mld_dsub_aply.f90 @@ -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 ! to compute an approximate solution of a linear system. ! - ! Note: trans is always 'N' here. ! if (size(prec%av) < mld_ap_nd_) then diff --git a/mlprec/mld_zsub_aply.f90 b/mlprec/mld_zsub_aply.f90 index f4f69c69..9c64eff9 100644 --- a/mlprec/mld_zsub_aply.f90 +++ b/mlprec/mld_zsub_aply.f90 @@ -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 ! to compute an approximate solution of a linear system. ! - ! Note: trans is always 'N' here. - ! if (size(prec%av) < mld_ap_nd_) then info = 4011