portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-14 17:55:38 +00:00
Code Releases Activity
eb1759d7ea
Grid/lib/qcd/action/fermion/WilsonKernels.h

60 lines
2.4 KiB
C
Raw Normal View History

Large scale change to support 5d fermion formulations. Have 5d replicated wilson with 4d gauge working and matrix regressing to Ls copies of wilson.
2015-05-31 15:09:02 +01:00
#ifndef GRID_QCD_DHOP_H
#define GRID_QCD_DHOP_H
namespace Grid {
namespace QCD {
////////////////////////////////////////////////////////////////////////////////////////////////////////////////
Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types.
2015-08-10 20:47:44 +01:00
// Helper routines that implement Wilson stencil for a single site.
* Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future.
2015-08-15 23:25:49 +01:00
// Common to both the WilsonFermion and WilsonFermion5D
Large scale change to support 5d fermion formulations. Have 5d replicated wilson with 4d gauge working and matrix regressing to Ls copies of wilson.
2015-05-31 15:09:02 +01:00
////////////////////////////////////////////////////////////////////////////////////////////////////////////////
Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types.
2015-08-10 20:47:44 +01:00
* Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future.
2015-08-15 23:25:49 +01:00
template<class Impl> class WilsonKernels : public FermionOperator<Impl> {
Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types.
2015-08-10 20:47:44 +01:00
public:
* Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future.
2015-08-15 23:25:49 +01:00
INHERIT_IMPL_TYPES(Impl);
typedef FermionOperator<Impl> Base;
Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types.
2015-08-10 20:47:44 +01:00
public:
Stencil class name global change to StencilImpl typedef
2015-11-06 11:30:17 +00:00
void DiracOptDhopSite(StencilImpl &st,DoubledGaugeField &U,
* Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future.
2015-08-15 23:25:49 +01:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in, FermionField &out);
big commit fixing nocompiles in defective C++11 compilers (gcc, icpc). stared getting to near the bleeding edge I guess
2015-06-30 15:03:11 +01:00
Stencil class name global change to StencilImpl typedef
2015-11-06 11:30:17 +00:00
void DiracOptDhopSiteDag(StencilImpl &st,DoubledGaugeField &U,
* Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future.
2015-08-15 23:25:49 +01:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in,FermionField &out);
Stencil class name global change to StencilImpl typedef
2015-11-06 11:30:17 +00:00
void DiracOptDhopDir(StencilImpl &st,DoubledGaugeField &U,
* Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future.
2015-08-15 23:25:49 +01:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in, FermionField &out,int dirdisp,int gamma);
Fix a regression failure on Mobius; chroma regression added
2015-12-10 22:55:00 +00:00
Stencil change
2015-11-07 08:06:31 +00:00
void DiracOptAsmDhopSite(StencilImpl &st,DoubledGaugeField &U,
Changes for AVX512 assembler
2015-11-06 11:45:51 +00:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in, FermionField &out,uint64_t *);
Merge branch 'master' of https://github.com/paboyle/Grid into scidac1_2 Conflicts: lib/qcd/action/fermion/WilsonKernels.h tests/Make.inc
2015-12-15 16:11:59 +00:00
#if 0
Stencil change
2015-11-07 08:06:31 +00:00
void DiracOptAsmDhopSite(StencilImpl &st,DoubledGaugeField &U,
Changes for AVX512 assembler
2015-11-06 11:45:51 +00:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in, FermionField &out,uint64_t *p){
DiracOptDhopSite(st,U,buf,sF,sU,in,out); // will template override for Wilson Nc=3
}
#endif
Fixes for Gparity comparison with CPS (Instantiation, Gamma matrix convention)
2015-12-07 07:04:57 +00:00
// doesn't seem to work with Gparity at the moment
#undef HANDOPT
Merge branch 'master' of https://github.com/paboyle/Grid into scidac1_2 Conflicts: lib/qcd/action/fermion/WilsonKernels.h tests/Make.inc
2015-12-15 16:11:59 +00:00
//#define HANDOPT
Stencil class name global change to StencilImpl typedef
2015-11-06 11:30:17 +00:00
void DiracOptHandDhopSite(StencilImpl &st,DoubledGaugeField &U,
Reintroducing the hand unrolled loops
2015-09-08 17:45:30 +01:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in, FermionField &out);
Stencil class name global change to StencilImpl typedef
2015-11-06 11:30:17 +00:00
void DiracOptHandDhopSiteDag(StencilImpl &st,DoubledGaugeField &U,
Reintroducing the hand unrolled loops
2015-09-08 17:45:30 +01:00
std::vector<SiteHalfSpinor,alignedAllocator<SiteHalfSpinor> > &buf,
int sF,int sU,const FermionField &in, FermionField &out);
Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types.
2015-08-10 20:47:44 +01:00
Fix a regression failure on Mobius; chroma regression added
2015-12-10 22:55:00 +00:00
WilsonKernels(const ImplParams &p= ImplParams());
Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types.
2015-08-10 20:47:44 +01:00
};
Large scale change to support 5d fermion formulations. Have 5d replicated wilson with 4d gauge working and matrix regressing to Ls copies of wilson.
2015-05-31 15:09:02 +01:00
}
}
#endif
Copy Permalink