From 7d4264176afb7614f248e96cbb5615d5bca9b4de Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Fri, 6 Oct 2017 14:11:12 +0100 Subject: [PATCH] Fix zeroing of tx in as_apply_vect. --- mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 | 2 +- mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 | 2 +- mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 | 2 +- mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 b/mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 index 493b6482..3c74d6ae 100644 --- a/mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 @@ -127,7 +127,7 @@ subroutine mld_c_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,trans,& call psb_geasb(tx,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ty,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) - + call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE ! this will save one SPMM when INIT=Z, and will be diff --git a/mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 b/mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 index 3f5d4c5b..3063e398 100644 --- a/mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 @@ -127,7 +127,7 @@ subroutine mld_d_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,trans,& call psb_geasb(tx,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ty,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) - + call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE ! this will save one SPMM when INIT=Z, and will be diff --git a/mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 b/mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 index a0fcd990..ea1d20ba 100644 --- a/mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 @@ -127,7 +127,7 @@ subroutine mld_s_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,trans,& call psb_geasb(tx,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ty,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) - + call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE ! this will save one SPMM when INIT=Z, and will be diff --git a/mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 b/mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 index 45618af0..a0fe467b 100644 --- a/mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 @@ -127,7 +127,7 @@ subroutine mld_z_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,trans,& call psb_geasb(tx,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ty,sm%desc_data,info,mold=x%v,scratch=.true.) call psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) - + call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE ! this will save one SPMM when INIT=Z, and will be