Comment for fixes to apply_vect.

7 years ago · 4c70174c9f
parent 7d4264176a
commit 4c70174c9f
4 changed files with 8 additions and 4 deletions
--- 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
--- 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
--- 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
--- 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