From 4c70174c9fc799d6ff51f67f0529da357c074907 Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Sat, 7 Oct 2017 08:01:28 +0100 Subject: [PATCH] Comment for fixes to apply_vect. --- mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 | 3 ++- mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 | 3 ++- mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 | 3 ++- mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 | 3 ++- 4 files changed, 8 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 3c74d6ae..1f442af4 100644 --- a/mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_c_as_smoother_apply_vect.f90 @@ -126,7 +126,8 @@ 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 psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) + ! Need to zero tx because of the apply_restr call. call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE 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 3063e398..2c8114a7 100644 --- a/mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_d_as_smoother_apply_vect.f90 @@ -126,7 +126,8 @@ 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 psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) + ! Need to zero tx because of the apply_restr call. call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE 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 ea1d20ba..9b89d844 100644 --- a/mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_s_as_smoother_apply_vect.f90 @@ -126,7 +126,8 @@ 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 psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) + ! Need to zero tx because of the apply_restr call. call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE 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 a0fe467b..bb3cce93 100644 --- a/mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 +++ b/mlprec/impl/smoother/mld_z_as_smoother_apply_vect.f90 @@ -126,7 +126,8 @@ 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 psb_geasb(ww,sm%desc_data,info,mold=x%v,scratch=.true.) + ! Need to zero tx because of the apply_restr call. call tx%zero() ! ! Unroll the first iteration and fold it inside SELECT CASE