portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 15:55:37 +00:00
Code Releases Activity

Correction of the dagger version of the Clover

Browse Source
This commit is contained in:
Guido Cossu 2017-08-15 10:50:44 +01:00
parent 75ee6cfc86
commit 8d168ded4a
1 changed files with 9 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
lib/qcd/action/fermion/WilsonCloverFermion.cc
View File

@ -142,19 +142,22 @@ void WilsonCloverFermion<Impl>::Mooee(const FermionField &in, FermionField &out)
template <class Impl> template <class Impl>
void WilsonCloverFermion<Impl>::MooeeDag(const FermionField &in, FermionField &out) 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> template <class Impl>
void WilsonCloverFermion<Impl>::MooeeInv(const FermionField &in, FermionField &out) void WilsonCloverFermion<Impl>::MooeeInv(const FermionField &in, FermionField &out)
{ {
conformable(in,out);
this->MooeeInternal(in, out, DaggerNo, InverseYes); this->MooeeInternal(in, out, DaggerNo, InverseYes);
} }
template <class Impl> template <class Impl>
void WilsonCloverFermion<Impl>::MooeeInvDag(const FermionField &in, FermionField &out) 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> template <class Impl>
@ -183,7 +186,10 @@ void WilsonCloverFermion<Impl>::MooeeInternal(const FermionField &in, FermionFie
} }
std::cout << GridLogMessage << "*Clover.checkerboard " << (*Clover).checkerboard << std::endl; 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 } // MooeeInternal
// Derivative parts // Derivative parts