From 3cb1b545d00fe576a63d8fd437ae9731a2f95578 Mon Sep 17 00:00:00 2001
From: paboyle <paboyle@ph.ed.ac.uk>
Date: Tue, 6 Mar 2018 04:13:26 -0800
Subject: [PATCH] Don't alias the variables with a where statement.

---
 lib/qcd/utils/CovariantCshift.h | 21 ++++++++++++---------
 1 file changed, 12 insertions(+), 9 deletions(-)

diff --git a/lib/qcd/utils/CovariantCshift.h b/lib/qcd/utils/CovariantCshift.h
index ed96f3bf..ed4e2d3f 100644
--- a/lib/qcd/utils/CovariantCshift.h
+++ b/lib/qcd/utils/CovariantCshift.h
@@ -95,11 +95,13 @@ namespace ConjugateBC {
 
     Lattice<iScalar<vInteger> > coor(grid);    LatticeCoordinate(coor,mu);
 
-    Lattice<covariant> field_bc = Cshift(field,mu,1);// moves towards negative mu;
+    Lattice<covariant> tmp1(grid);
+    Lattice<covariant> tmp2(grid);
 
-    field_bc = where(coor==Lmu,conjugate(field_bc),field_bc);
-    //    std::cout<<"Gparity::CovCshiftForward mu="<<mu<<std::endl;
-    return Link*field_bc;
+    tmp1 = Cshift(field,mu,1);// moves towards negative mu;
+    tmp2 = where(coor==Lmu,conjugate(tmp1),tmp1);
+
+    return Link*tmp2;
   }
 
   template<class covariant,class gauge> Lattice<covariant> CovShiftBackward(const Lattice<gauge> &Link, 
@@ -114,12 +116,13 @@ namespace ConjugateBC {
 
     Lattice<iScalar<vInteger> > coor(grid);    LatticeCoordinate(coor,mu);
 
-    Lattice<covariant> tmp(grid);
+    Lattice<covariant> tmp1(grid);
+    Lattice<covariant> tmp2(grid);
 
-    tmp = adj(Link)*field;
-    tmp = where(coor==Lmu,conjugate(tmp),tmp);
-    //    std::cout<<"Gparity::CovCshiftBackward mu="<<mu<<std::endl;
-    return Cshift(tmp,mu,-1);// moves towards positive mu
+    tmp1 = adj(Link)*field;
+    tmp2 = where(coor==Lmu,conjugate(tmp1),tmp1);
+
+    return Cshift(tmp2,mu,-1);// moves towards positive mu
   }
 
 }