diff --git a/Grid/qcd/action/pseudofermion/ExactOneFlavourRatio.h b/Grid/qcd/action/pseudofermion/ExactOneFlavourRatio.h index 25285565..c6746a88 100644 --- a/Grid/qcd/action/pseudofermion/ExactOneFlavourRatio.h +++ b/Grid/qcd/action/pseudofermion/ExactOneFlavourRatio.h @@ -66,6 +66,7 @@ namespace QCD{ FermionField Phi; // the pseudofermion field for this trajectory public: + ExactOneFlavourRatioPseudoFermionAction(AbstractEOFAFermion& _Lop, AbstractEOFAFermion& _Rop, OperatorFunction& HeatbathCG, diff --git a/HMC/Mobius2p1fEOFA.cc b/HMC/Mobius2p1fEOFA.cc index 61b06829..997e76ab 100644 --- a/HMC/Mobius2p1fEOFA.cc +++ b/HMC/Mobius2p1fEOFA.cc @@ -30,7 +30,9 @@ directory /* END LEGAL */ #include +#ifdef GRID_DEFAULT_PRECISION_DOUBLE #define MIXED_PRECISION +#endif namespace Grid{ namespace QCD{ @@ -346,6 +348,7 @@ int main(int argc, char **argv) { #else ExactOneFlavourRatioPseudoFermionAction EOFA(Strange_Op_L, Strange_Op_R, + ActionCG, ActionCG, ActionCG, DerivativeCG, DerivativeCG, OFRp, true); diff --git a/tests/forces/Test_dwf_force_eofa.cc b/tests/forces/Test_dwf_force_eofa.cc index f17579ae..18e339ab 100644 --- a/tests/forces/Test_dwf_force_eofa.cc +++ b/tests/forces/Test_dwf_force_eofa.cc @@ -84,7 +84,7 @@ int main (int argc, char** argv) DomainWallEOFAFermionR Rop(U, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mb, mf, mb, -1.0, 1, M5); OneFlavourRationalParams Params(0.95, 100.0, 5000, 1.0e-12, 12); ConjugateGradient CG(1.0e-12, 5000); - ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, Params, true); + ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, CG, CG, CG, CG, Params, true); Meofa.refresh(U, RNG5); RealD S = Meofa.S(U); // pdag M p diff --git a/tests/forces/Test_dwf_gpforce_eofa.cc b/tests/forces/Test_dwf_gpforce_eofa.cc index 3afeaa43..32b3ed6b 100644 --- a/tests/forces/Test_dwf_gpforce_eofa.cc +++ b/tests/forces/Test_dwf_gpforce_eofa.cc @@ -89,7 +89,7 @@ int main (int argc, char** argv) FermionAction Rop(U, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mb, mf, mb, -1.0, 1, M5, params); OneFlavourRationalParams Params(0.95, 100.0, 5000, 1.0e-12, 12); ConjugateGradient CG(1.0e-12, 5000); - ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, Params, true); + ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, CG, CG, CG, CG, Params, true); Meofa.refresh(U, RNG5); RealD S = Meofa.S(U); // pdag M p diff --git a/tests/forces/Test_mobius_force_eofa.cc b/tests/forces/Test_mobius_force_eofa.cc index 2a5a7d04..68a2931f 100644 --- a/tests/forces/Test_mobius_force_eofa.cc +++ b/tests/forces/Test_mobius_force_eofa.cc @@ -86,7 +86,7 @@ int main (int argc, char** argv) MobiusEOFAFermionR Rop(U, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mb, mf, mb, -1.0, 1, M5, b, c); OneFlavourRationalParams Params(0.95, 100.0, 5000, 1.0e-12, 12); ConjugateGradient CG(1.0e-12, 5000); - ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, Params, false); + ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, CG, CG, CG, CG, Params, false); Meofa.refresh(U, RNG5); RealD S = Meofa.S(U); // pdag M p diff --git a/tests/forces/Test_mobius_gpforce_eofa.cc b/tests/forces/Test_mobius_gpforce_eofa.cc index 72f1dee2..577b4477 100644 --- a/tests/forces/Test_mobius_gpforce_eofa.cc +++ b/tests/forces/Test_mobius_gpforce_eofa.cc @@ -91,7 +91,7 @@ int main (int argc, char** argv) FermionAction Rop(U, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mb, mf, mb, -1.0, 1, M5, b, c, params); OneFlavourRationalParams Params(0.95, 100.0, 5000, 1.0e-12, 12); ConjugateGradient CG(1.0e-12, 5000); - ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, Params, false); + ExactOneFlavourRatioPseudoFermionAction Meofa(Lop, Rop, CG, CG, CG, CG, CG, Params, false); Meofa.refresh(U, RNG5); RealD S = Meofa.S(U); // pdag M p