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::Mooee(const FermionField &in, FermionField &out) template void WilsonCloverFermion::MooeeDag(const FermionField &in, FermionField &out) { - this->MooeeInternal(in, out, DaggerNo, InverseYes); + conformable(in,out); + this->MooeeInternal(in, out, DaggerYes, InverseYes); } template void WilsonCloverFermion::MooeeInv(const FermionField &in, FermionField &out) { + conformable(in,out); this->MooeeInternal(in, out, DaggerNo, InverseYes); } template void WilsonCloverFermion::MooeeInvDag(const FermionField &in, FermionField &out) { - this->MooeeInternal(in, out, DaggerNo, InverseYes); + conformable(in,out); + this->MooeeInternal(in, out, DaggerYes, InverseYes); } template @@ -183,7 +186,10 @@ void WilsonCloverFermion::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