diff --git a/lib/qcd/action/Actions.h b/lib/qcd/action/Actions.h index 819e4f10..d4a605aa 100644 --- a/lib/qcd/action/Actions.h +++ b/lib/qcd/action/Actions.h @@ -206,15 +206,17 @@ typedef OverlapWilsonPartialFractionZolotarevFermion OverlapWilsonP typedef WilsonFermion GparityWilsonFermionR; typedef WilsonFermion GparityWilsonFermionF; typedef WilsonFermion GparityWilsonFermionD; -typedef WilsonTMFermion GparityWilsonTMFermionR; -typedef WilsonTMFermion GparityWilsonTMFermionF; -typedef WilsonTMFermion GparityWilsonTMFermionD; typedef DomainWallFermion GparityDomainWallFermionR; typedef DomainWallFermion GparityDomainWallFermionF; typedef DomainWallFermion GparityDomainWallFermionD; +#if 0 +typedef WilsonTMFermion GparityWilsonTMFermionR; +typedef WilsonTMFermion GparityWilsonTMFermionF; +typedef WilsonTMFermion GparityWilsonTMFermionD; typedef MobiusFermion GparityMobiusFermionR; typedef MobiusFermion GparityMobiusFermionF; typedef MobiusFermion GparityMobiusFermionD; +#endif }} /////////////////////////////////////////////////////////////////////////////// diff --git a/lib/qcd/action/fermion/ScaledShamirFermion.h b/lib/qcd/action/fermion/ScaledShamirFermion.h index 53acb694..c0deef3e 100644 --- a/lib/qcd/action/fermion/ScaledShamirFermion.h +++ b/lib/qcd/action/fermion/ScaledShamirFermion.h @@ -48,6 +48,7 @@ namespace Grid { GridCartesian &FourDimGrid, GridRedBlackCartesian &FourDimRedBlackGrid, RealD _mass,RealD _M5, +// RealD scale): RealD scale,const ImplParams &p= ImplParams()) : // b+c=scale, b-c = 1 <=> 2b = scale+1; 2c = scale-1 @@ -55,7 +56,8 @@ namespace Grid { FiveDimGrid, FiveDimRedBlackGrid, FourDimGrid, - FourDimRedBlackGrid,_mass,_M5,0.5*(scale+1.0),0.5*(scale-1.0),p) + FourDimRedBlackGrid,_mass,_M5,0.5*(scale+1.0),0.5*(scale-1.0),p) +// FourDimRedBlackGrid,_mass,_M5,0.5*(scale+1.0),0.5*(scale-1.0)) { } diff --git a/lib/qcd/action/fermion/WilsonKernels.cc b/lib/qcd/action/fermion/WilsonKernels.cc index 37b8577f..65013a17 100644 --- a/lib/qcd/action/fermion/WilsonKernels.cc +++ b/lib/qcd/action/fermion/WilsonKernels.cc @@ -596,7 +596,8 @@ void WilsonKernels::DiracOptDhopDir(StencilImpl &st,DoubledGaugeField &U, vstream(out._odata[sF],result*(-0.5)); } -#if ( ! defined(AVX512) ) && ( ! defined(IMCI) ) +//#if ( ! defined(AVX512) ) && ( ! defined(IMCI) ) +#if 1 template void WilsonKernels::DiracOptAsmDhopSite(StencilImpl &st,DoubledGaugeField &U, std::vector > &buf, diff --git a/lib/qcd/action/fermion/WilsonKernels.h b/lib/qcd/action/fermion/WilsonKernels.h index 1308aec1..d48d5f88 100644 --- a/lib/qcd/action/fermion/WilsonKernels.h +++ b/lib/qcd/action/fermion/WilsonKernels.h @@ -60,7 +60,7 @@ namespace Grid { void DiracOptAsmDhopSite(StencilImpl &st,DoubledGaugeField &U, std::vector > &buf, -#if 1 +#if 0 //<<<<<<< HEAD int sF,int sU,const FermionField &in, FermionField &out,bool local= true, bool nonlocal=true); // int sF,int sU,const FermionField &in, FermionField &out,uint64_t *); diff --git a/lib/qcd/action/fermion/WilsonKernelsHand.cc b/lib/qcd/action/fermion/WilsonKernelsHand.cc index d3d0a6bc..a2e0622c 100644 --- a/lib/qcd/action/fermion/WilsonKernelsHand.cc +++ b/lib/qcd/action/fermion/WilsonKernelsHand.cc @@ -1091,12 +1091,17 @@ void WilsonKernels::DiracOptHandDhopSite(StencilImpl &st,DoubledGaugeFiel //////////////////////////////////////////////// // Specialise Gparity to simple implementation //////////////////////////////////////////////// -template<> +//template +//int WilsonKernels::DiracOptHandDhopSiteDag(StencilImpl &st,DoubledGaugeField &U, +// std::vector > &buf, +// int ss,int sU,const FermionField &in, FermionField &out, bool Local, bool Nonlocal) +#if 0 +template void WilsonKernels::DiracOptHandDhopSite(StencilImpl &st,DoubledGaugeField &U, std::vector > &buf, int sF,int sU,const FermionField &in, FermionField &out, bool Local, bool Nonlocal) { - DiracOptDhopSite(st,U,buf,sF,sU,in,out); // will template override for Wilson Nc=3 + DiracOptDhopSite(st,U,buf,sF,sU,in,out,Local,Nonlocal); // will template override for Wilson Nc=3 } template<> @@ -1104,7 +1109,7 @@ void WilsonKernels::DiracOptHandDhopSiteDag(StencilImpl &st, std::vector > &buf, int sF,int sU,const FermionField &in, FermionField &out, bool Local, bool Nonlocal) { - DiracOptDhopSiteDag(st,U,buf,sF,sU,in,out); // will template override for Wilson Nc=3 + DiracOptDhopSiteDag(st,U,buf,sF,sU,in,out,Local,Nonlocal); // will template override for Wilson Nc=3 } template<> @@ -1112,7 +1117,7 @@ void WilsonKernels::DiracOptHandDhopSite(StencilImpl &st,Dou std::vector > &buf, int sF,int sU,const FermionField &in, FermionField &out, bool Local, bool Nonlocal) { - DiracOptDhopSite(st,U,buf,sF,sU,in,out); // will template override for Wilson Nc=3 + DiracOptDhopSite(st,U,buf,sF,sU,in,out,Local,Nonlocal); // will template override for Wilson Nc=3 } template<> @@ -1120,8 +1125,9 @@ void WilsonKernels::DiracOptHandDhopSiteDag(StencilImpl &st, std::vector > &buf, int sF,int sU,const FermionField &in, FermionField &out, bool Local, bool Nonlocal) { - DiracOptDhopSiteDag(st,U,buf,sF,sU,in,out); // will template override for Wilson Nc=3 + DiracOptDhopSiteDag(st,U,buf,sF,sU,in,out,Local,Nonlocal); // will template override for Wilson Nc=3 } +#endif