portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 19:25:56 +01:00
Code Releases Activity

Preparing for FTHMC

Browse Source
This commit is contained in:
Peter Boyle 2023-05-19 21:21:55 -04:00
parent d418347d86
commit 4240ad5ca8
1 changed files with 33 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

48
Grid/qcd/smearing/GaugeConfiguration.h
View File

@ -7,26 +7,40 @@
NAMESPACE_BEGIN(Grid); NAMESPACE_BEGIN(Grid);
//trivial class for no smearing
template< class Impl > template< class Impl >
class NoSmearing class ConfigurationBase
{ {
public: public:
INHERIT_FIELD_TYPES(Impl); INHERIT_FIELD_TYPES(Impl);
Field* ThinField; ConfigurationBase() {}
virtual ~ConfigurationBase() {}
virtual void set_Field(Field& U) =0;
virtual void smeared_force(Field&) const = 0;
virtual Field& get_SmearedU() =0;
virtual Field &get_U(bool smeared = false) = 0;
};
NoSmearing(): ThinField(NULL) {} //trivial class for no smearing
template< class Impl >
class NoSmearing : public ConfigurationBase<Impl>
{
public:
INHERIT_FIELD_TYPES(Impl);
void set_Field(Field& U) { ThinField = &U; } Field* ThinLinks;
NoSmearing(): ThinLinks(NULL) {}
void set_Field(Field& U) { ThinLinks = &U; }
void smeared_force(Field&) const {} void smeared_force(Field&) const {}
Field& get_SmearedU() { return *ThinField; } Field& get_SmearedU() { return *ThinLinks; }
Field &get_U(bool smeared = false) Field &get_U(bool smeared = false)
{ {
return *ThinField; return *ThinLinks;
} }
}; };
@ -42,7 +56,7 @@ public:
It stores a list of smeared configurations. It stores a list of smeared configurations.
*/ */
template <class Gimpl> template <class Gimpl>
class SmearedConfiguration class SmearedConfiguration : public ConfigurationBase<Impl>
{ {
public: public:
INHERIT_GIMPL_TYPES(Gimpl); INHERIT_GIMPL_TYPES(Gimpl);
@ -51,10 +65,15 @@ private:
const unsigned int smearingLevels; const unsigned int smearingLevels;
Smear_Stout<Gimpl> *StoutSmearing; Smear_Stout<Gimpl> *StoutSmearing;
std::vector<GaugeField> SmearedSet; std::vector<GaugeField> SmearedSet;
public:
GaugeField* ThinLinks; /* Pointer to the thin links configuration */ // move to base???
private:
// Member functions // Member functions
//==================================================================== //====================================================================
void fill_smearedSet(GaugeField &U)
// Overridden in masked version
virtual void fill_smearedSet(GaugeField &U)
{ {
ThinLinks = &U; // attach the smearing routine to the field U ThinLinks = &U; // attach the smearing routine to the field U
@ -82,9 +101,10 @@ private:
} }
} }
} }
//====================================================================
GaugeField AnalyticSmearedForce(const GaugeField& SigmaKPrime, //overridden in masked verson
const GaugeField& GaugeK) const virtual GaugeField AnalyticSmearedForce(const GaugeField& SigmaKPrime,
const GaugeField& GaugeK) const
{ {
GridBase* grid = GaugeK.Grid(); GridBase* grid = GaugeK.Grid();
GaugeField C(grid), SigmaK(grid), iLambda(grid); GaugeField C(grid), SigmaK(grid), iLambda(grid);
@ -213,8 +233,6 @@ private:
//==================================================================== //====================================================================
public: public:
GaugeField*
ThinLinks; /* Pointer to the thin links configuration */
/* Standard constructor */ /* Standard constructor */
SmearedConfiguration(GridCartesian* UGrid, unsigned int Nsmear, SmearedConfiguration(GridCartesian* UGrid, unsigned int Nsmear,