From 771235017d5519724b8194848b9c040f747af876 Mon Sep 17 00:00:00 2001 From: neo Date: Fri, 19 Feb 2016 15:30:41 +0900 Subject: [PATCH] Adding smearing routines (development) --- lib/qcd/action/fermion/FermionOperatorImpl.h | 4 +- lib/qcd/action/gauge/GaugeImpl.h | 1 - lib/qcd/hmc/HMC.h | 2 +- lib/qcd/hmc/HmcRunner.h | 6 +- lib/qcd/smearing/GaugeConfiguration.h | 96 ++++++++++++++++++++ 5 files changed, 102 insertions(+), 7 deletions(-) create mode 100644 lib/qcd/smearing/GaugeConfiguration.h diff --git a/lib/qcd/action/fermion/FermionOperatorImpl.h b/lib/qcd/action/fermion/FermionOperatorImpl.h index 05994a39..312a8e41 100644 --- a/lib/qcd/action/fermion/FermionOperatorImpl.h +++ b/lib/qcd/action/fermion/FermionOperatorImpl.h @@ -75,7 +75,7 @@ namespace Grid { // // // template - // class MyOp : pubic { + // class MyOp : public { // public: // // INHERIT_ALL_IMPL_TYPES(Impl); @@ -110,7 +110,7 @@ namespace Grid { // Single flavour four spinors with colour index /////// template - class WilsonImpl : public PeriodicGaugeImpl< GaugeImplTypes< S,Nrepresentation> > { + class WilsonImpl : public PeriodicGaugeImpl< GaugeImplTypes< S, Nrepresentation> > { public: typedef PeriodicGaugeImpl< GaugeImplTypes< S,Nrepresentation> > Gimpl; diff --git a/lib/qcd/action/gauge/GaugeImpl.h b/lib/qcd/action/gauge/GaugeImpl.h index 51a4675e..02a8017b 100644 --- a/lib/qcd/action/gauge/GaugeImpl.h +++ b/lib/qcd/action/gauge/GaugeImpl.h @@ -44,7 +44,6 @@ template class WilsonLoops; typedef typename GImpl::GaugeLinkField GaugeLinkField;\ typedef typename GImpl::GaugeField GaugeField; - // template class GaugeImplTypes { diff --git a/lib/qcd/hmc/HMC.h b/lib/qcd/hmc/HMC.h index 9e762832..0b87b03e 100644 --- a/lib/qcd/hmc/HMC.h +++ b/lib/qcd/hmc/HMC.h @@ -44,7 +44,7 @@ Author: paboyle namespace Grid{ namespace QCD{ - + struct HMCparameters{ diff --git a/lib/qcd/hmc/HmcRunner.h b/lib/qcd/hmc/HmcRunner.h index 4d8bb2dd..ce92b1d2 100644 --- a/lib/qcd/hmc/HmcRunner.h +++ b/lib/qcd/hmc/HmcRunner.h @@ -47,7 +47,7 @@ public: GridRedBlackCartesian * UrbGrid ; GridRedBlackCartesian * FrbGrid ; - virtual void BuildTheAction (int argc, char **argv) = 0; + virtual void BuildTheAction (int argc, char **argv) = 0; // necessary? void Run (int argc, char **argv){ @@ -96,7 +96,7 @@ public: GridSerialRNG sRNG; GridParallelRNG pRNG(UGrid); - LatticeGaugeField U(UGrid); + LatticeGaugeField U(UGrid); // change this to an extended field (smearing class) std::vector SerSeed({1,2,3,4,5}); std::vector ParSeed({6,7,8,9,10}); @@ -129,7 +129,7 @@ public: Checkpoint.CheckpointRestore(StartTraj, U, sRNG, pRNG); } - HybridMonteCarlo HMC(HMCpar, MDynamics,sRNG,pRNG,U); + HybridMonteCarlo HMC(HMCpar, MDynamics,sRNG,pRNG,U); // pass the extended field HMC.AddObservable(&Checkpoint); HMC.AddObservable(&PlaqLog); diff --git a/lib/qcd/smearing/GaugeConfiguration.h b/lib/qcd/smearing/GaugeConfiguration.h new file mode 100644 index 00000000..d040b739 --- /dev/null +++ b/lib/qcd/smearing/GaugeConfiguration.h @@ -0,0 +1,96 @@ +/*! + @file GaugeConfiguration.h + + @brief Declares the GaugeConfiguration class +*/ +#ifndef GAUGE_CONFIG_ +#define GAUGE_CONFIG_ + +namespace Grid { + + namespace QCD { + + /*! + @brief Smeared configuration container + + It will behave like a configuration from the point of view of + the HMC update and integrators. + An "advanced configuration" object that can provide not only the + data to store the gauge configuration but also operations to manipulate + it like smearing. + + It stores a list of smeared configurations. + */ + template + class GaugeConfiguration { + public: + INHERIT_GIMPL_TYPES(Gimpl) + private: + const unsigned int smearingLevels; + Smear_Stout StoutSmearing; + std::vector SmearedSet; + + // Member functions + void fill_smearedSet(); + GaugeField AnalyticSmearedForce(const GaugeField&, + const GaugeField&) const; + const GaugeField& get_smeared_conf(int) const; + + void set_iLambda(GaugeField& iLambda, + GaugeField& e_iQ, + const GaugeField& iQ, + const GaugeField& Sigmap, + const GaugeField& U)const; + + /* Check these types (do I need to pass iQ1,2 ? ) + void set_uw(RealD& u, RealD& w, + const SUNmat& iQ1, const SUNmat& iQ2)const ; + void set_fj(ComplexD& f0, ComplexD& f1, + CompledD& f2, const RealD& u, + const RealD& w)const; + */ + + RealD func_xi0(RealD w)const; + RealD func_xi1(RealD w)const; + + public: + GaugeField* ThinLinks; /*!< @brief Pointer to the thin + links configuration */ + + /*! @brief Standard constructor */ + GaugeConfiguration(GridCartesian * UGrid, + unsigned int Nsmear, + Smear_Stout& Stout): + smearingLevels(Nsmear), + StoutSmearing(Stout), + ThinLinks(new GaugeField){ + for (unsigned int i=0; i< smearingLevels; ++i) + SmearedSet.push_back(*(new GaugeField(UGrid))); + } + + /*! For just thin links */ + GaugeConfiguration(GridCartesian * UGrid): + smearingLevels(0), + StoutSmearing(), + SmearedSet(0), + ThinLinks(new GaugeField(UGrid)){} + + void set_GaugeField(){ fill_smearedSet(); } + void smeared_force(GaugeField&) const; + GaugeField& get_current_conf() const; + GaugeField& select_conf(bool smeared) const { + if (smeared){ + if (smearingLevels) return get_current_conf(); + else return ThinLinks; + } + else return ThinLinks; + } + + }; + + + } + +} + +#endif