From 97019d2997c72165448553186b86802a28469be6 Mon Sep 17 00:00:00 2001 From: paboyle Date: Sun, 14 Jan 2018 23:39:57 +0000 Subject: [PATCH] Namespace, format --- .../fermion/ContinuedFractionFermion5D.h | 105 +++++++++--------- 1 file changed, 51 insertions(+), 54 deletions(-) diff --git a/lib/qcd/action/fermion/ContinuedFractionFermion5D.h b/lib/qcd/action/fermion/ContinuedFractionFermion5D.h index e1e50aa5..43ed4840 100644 --- a/lib/qcd/action/fermion/ContinuedFractionFermion5D.h +++ b/lib/qcd/action/fermion/ContinuedFractionFermion5D.h @@ -1,4 +1,4 @@ - /************************************************************************************* +/************************************************************************************* Grid physics library, www.github.com/paboyle/Grid @@ -24,76 +24,73 @@ Author: Peter Boyle 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. See the full license in the file "LICENSE" in the top level distribution directory - *************************************************************************************/ - /* END LEGAL */ +*************************************************************************************/ +/* END LEGAL */ #ifndef GRID_QCD_CONTINUED_FRACTION_H #define GRID_QCD_CONTINUED_FRACTION_H #include -namespace Grid { +NAMESPACE_BEGIN(Grid); - namespace QCD { +template +class ContinuedFractionFermion5D : public WilsonFermion5D +{ +public: + INHERIT_IMPL_TYPES(Impl); +public: - template - class ContinuedFractionFermion5D : public WilsonFermion5D - { - public: - INHERIT_IMPL_TYPES(Impl); - public: + // override multiply + virtual RealD M (const FermionField &in, FermionField &out); + virtual RealD Mdag (const FermionField &in, FermionField &out); - // override multiply - virtual RealD M (const FermionField &in, FermionField &out); - virtual RealD Mdag (const FermionField &in, FermionField &out); + // half checkerboard operaions + virtual void Meooe (const FermionField &in, FermionField &out); + virtual void MeooeDag (const FermionField &in, FermionField &out); + virtual void Mooee (const FermionField &in, FermionField &out); + virtual void MooeeDag (const FermionField &in, FermionField &out); + virtual void MooeeInv (const FermionField &in, FermionField &out); + virtual void MooeeInvDag (const FermionField &in, FermionField &out); - // half checkerboard operaions - virtual void Meooe (const FermionField &in, FermionField &out); - virtual void MeooeDag (const FermionField &in, FermionField &out); - virtual void Mooee (const FermionField &in, FermionField &out); - virtual void MooeeDag (const FermionField &in, FermionField &out); - virtual void MooeeInv (const FermionField &in, FermionField &out); - virtual void MooeeInvDag (const FermionField &in, FermionField &out); + // force terms; five routines; default to Dhop on diagonal + virtual void MDeriv (GaugeField &mat,const FermionField &U,const FermionField &V,int dag); + virtual void MoeDeriv(GaugeField &mat,const FermionField &U,const FermionField &V,int dag); + virtual void MeoDeriv(GaugeField &mat,const FermionField &U,const FermionField &V,int dag); - // force terms; five routines; default to Dhop on diagonal - virtual void MDeriv (GaugeField &mat,const FermionField &U,const FermionField &V,int dag); - virtual void MoeDeriv(GaugeField &mat,const FermionField &U,const FermionField &V,int dag); - virtual void MeoDeriv(GaugeField &mat,const FermionField &U,const FermionField &V,int dag); + // virtual void Instantiatable(void)=0; + virtual void Instantiatable(void) =0; - // virtual void Instantiatable(void)=0; - virtual void Instantiatable(void) =0; + // Efficient support for multigrid coarsening + virtual void Mdir (const FermionField &in, FermionField &out,int dir,int disp); - // Efficient support for multigrid coarsening - virtual void Mdir (const FermionField &in, FermionField &out,int dir,int disp); + // Constructors + ContinuedFractionFermion5D(GaugeField &_Umu, + GridCartesian &FiveDimGrid, + GridRedBlackCartesian &FiveDimRedBlackGrid, + GridCartesian &FourDimGrid, + GridRedBlackCartesian &FourDimRedBlackGrid, + RealD _mass,RealD M5,const ImplParams &p= ImplParams()); - // Constructors - ContinuedFractionFermion5D(GaugeField &_Umu, - GridCartesian &FiveDimGrid, - GridRedBlackCartesian &FiveDimRedBlackGrid, - GridCartesian &FourDimGrid, - GridRedBlackCartesian &FourDimRedBlackGrid, - RealD _mass,RealD M5,const ImplParams &p= ImplParams()); +protected: - protected: + void SetCoefficientsTanh(Approx::zolotarev_data *zdata,RealD scale); + void SetCoefficientsZolotarev(RealD zolo_hi,Approx::zolotarev_data *zdata);; - void SetCoefficientsTanh(Approx::zolotarev_data *zdata,RealD scale); - void SetCoefficientsZolotarev(RealD zolo_hi,Approx::zolotarev_data *zdata);; + // Cont frac + RealD dw_diag; + RealD mass; + RealD R; + RealD ZoloHiInv; + std::vector Beta; + std::vector cc;; + std::vector cc_d;; + std::vector sqrt_cc; + std::vector See; + std::vector Aee; - // Cont frac - RealD dw_diag; - RealD mass; - RealD R; - RealD ZoloHiInv; - std::vector Beta; - std::vector cc;; - std::vector cc_d;; - std::vector sqrt_cc; - std::vector See; - std::vector Aee; - - }; +}; - } -} +NAMESPACE_END(Grid); #endif