From f7373e97a4f32e0b9b07b53f0656decb3b9e7610 Mon Sep 17 00:00:00 2001
From: gfilaci <gianluca.filaci@gmail.com>
Date: Mon, 16 Dec 2019 10:04:44 +0100
Subject: [PATCH] Missing conjugate in MooeeInvDag

---
 Grid/qcd/action/fermion/implementation/CayleyFermion5Dcache.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Grid/qcd/action/fermion/implementation/CayleyFermion5Dcache.h b/Grid/qcd/action/fermion/implementation/CayleyFermion5Dcache.h
index 35402994..dbdf134b 100644
--- a/Grid/qcd/action/fermion/implementation/CayleyFermion5Dcache.h
+++ b/Grid/qcd/action/fermion/implementation/CayleyFermion5Dcache.h
@@ -234,12 +234,12 @@ CayleyFermion5D<Impl>::MooeeInvDag (const FermionField &psi_i, FermionField &chi
     res = psi(ss+Ls-1) - conjugate(puee[Ls-2])*tmp - acc;
     
     // Apply L_m^{-\dagger} D^{-dagger} L^{-dagger}
-    res = (1.0/pdee[Ls-1])*res;
+    res = conjugate(1.0/pdee[Ls-1])*res;
     coalescedWrite(chi[ss+Ls-1],res);
     spProj5m(acc,res);
     spProj5p(tmp,res);
     for (int s=Ls-2;s>=0;s--){
-      res = (1.0/pdee[s])*chi(ss+s) - conjugate(plee[s])*tmp - conjugate(pleem[s])*acc;
+      res = conjugate(1.0/pdee[s])*chi(ss+s) - conjugate(plee[s])*tmp - conjugate(pleem[s])*acc;
       spProj5p(tmp,res);
       coalescedWrite(chi[ss+s],res);
     }