From 8d168ded4afeefdd6d3970f791c7c5d98f330b86 Mon Sep 17 00:00:00 2001
From: Guido Cossu <guido.cossu@ed.ac.uk>
Date: Tue, 15 Aug 2017 10:50:44 +0100
Subject: [PATCH] Correction of the dagger version of the Clover

---
 lib/qcd/action/fermion/WilsonCloverFermion.cc | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/lib/qcd/action/fermion/WilsonCloverFermion.cc b/lib/qcd/action/fermion/WilsonCloverFermion.cc
index f8b62ba4..ce120846 100644
--- a/lib/qcd/action/fermion/WilsonCloverFermion.cc
+++ b/lib/qcd/action/fermion/WilsonCloverFermion.cc
@@ -142,19 +142,22 @@ void WilsonCloverFermion<Impl>::Mooee(const FermionField &in, FermionField &out)
 template <class Impl>
 void WilsonCloverFermion<Impl>::MooeeDag(const FermionField &in, FermionField &out)
 {
-  this->MooeeInternal(in, out, DaggerNo, InverseYes);
+  conformable(in,out);
+  this->MooeeInternal(in, out, DaggerYes, InverseYes);
 }
 
 template <class Impl>
 void WilsonCloverFermion<Impl>::MooeeInv(const FermionField &in, FermionField &out)
 {
+  conformable(in,out);
   this->MooeeInternal(in, out, DaggerNo, InverseYes);
 }
 
 template <class Impl>
 void WilsonCloverFermion<Impl>::MooeeInvDag(const FermionField &in, FermionField &out)
 {
-  this->MooeeInternal(in, out, DaggerNo, InverseYes);
+  conformable(in,out);
+  this->MooeeInternal(in, out, DaggerYes, InverseYes);
 }
 
 template <class Impl>
@@ -183,7 +186,10 @@ void WilsonCloverFermion<Impl>::MooeeInternal(const FermionField &in, FermionFie
   }
 
   std::cout << GridLogMessage << "*Clover.checkerboard "  << (*Clover).checkerboard << std::endl;
-  if (dag){ out = adj(*Clover) * in;} else { out = *Clover * in;}
+  std::cout << GridLogMessage << "in.checkerboard "  << in.checkerboard << std::endl;
+  std::cout << GridLogMessage << "out.checkerboard "  << out.checkerboard << std::endl;
+  out = *Clover * in;
+  //if (dag){ out = adj(*Clover) * in;} else { out = *Clover * in;}
 } // MooeeInternal
 
 // Derivative parts