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

added floating-point support

Browse Source
This commit is contained in:
david clarke 2023-10-20 13:54:26 -06:00
parent 7bb8ab7000
commit 21ed6ac0f4
2 changed files with 19 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
Grid/qcd/smearing/HISQSmearing.h
View File

@ -81,7 +81,7 @@ struct SmearingParameters{
/*! @brief create fat links from link variables */ /*! @brief create fat links from link variables */
//template<class LGF, class Gimpl> //template<class GF, class Gimpl>
template<class Gimpl> template<class Gimpl>
class Smear_HISQ_fat : public Gimpl { class Smear_HISQ_fat : public Gimpl {
@ -92,7 +92,8 @@ private:
public: public:
INHERIT_GIMPL_TYPES(Gimpl); INHERIT_GIMPL_TYPES(Gimpl);
typedef typename Gimpl::GaugeField LGF; typedef typename Gimpl::GaugeField GF;
typedef typename Gimpl::GaugeLinkField LF;
// Don't allow default values here. // Don't allow default values here.
Smear_HISQ_fat(GridCartesian* grid, Real c1, Real cnaik, Real c3, Real c5, Real c7, Real clp) Smear_HISQ_fat(GridCartesian* grid, Real c1, Real cnaik, Real c3, Real c5, Real c7, Real clp)
@ -112,20 +113,20 @@ public:
~Smear_HISQ_fat() {} ~Smear_HISQ_fat() {}
void smear(LGF& u_smr, LGF& u_thin) const { void smear(GF& u_smr, GF& u_thin) const {
SmearingParameters lt = this->_linkTreatment; SmearingParameters lt = this->_linkTreatment;
// Create a padded cell of extra padding depth=1 and fill the padding. // Create a padded cell of extra padding depth=1 and fill the padding.
int depth = 1; int depth = 1;
PaddedCell Ghost(depth,this->_grid); PaddedCell Ghost(depth,this->_grid);
LGF Ughost = Ghost.Exchange(u_thin); GF Ughost = Ghost.Exchange(u_thin);
// This is where auxiliary N-link fields and the final smear will be stored. // This is where auxiliary N-link fields and the final smear will be stored.
LGF Ughost_fat(Ughost.Grid()); GF Ughost_fat(Ughost.Grid());
LGF Ughost_3link(Ughost.Grid()); GF Ughost_3link(Ughost.Grid());
LGF Ughost_5linkA(Ughost.Grid()); GF Ughost_5linkA(Ughost.Grid());
LGF Ughost_5linkB(Ughost.Grid()); GF Ughost_5linkB(Ughost.Grid());
// Create 3-link stencil. We allow mu==nu just to make the indexing easier. // Create 3-link stencil. We allow mu==nu just to make the indexing easier.
// Shifts with mu==nu will not be used. // Shifts with mu==nu will not be used.
@ -279,8 +280,8 @@ public:
u_smr = Ghost.Extract(Ughost_fat) + lt.c_1*u_thin; u_smr = Ghost.Extract(Ughost_fat) + lt.c_1*u_thin;
// Load up U and V std::vectors to access thin and smeared links. // Load up U and V std::vectors to access thin and smeared links.
std::vector<LatticeColourMatrix> U(Nd, u_thin.Grid()); std::vector<LF> U(Nd, u_thin.Grid());
std::vector<LatticeColourMatrix> V(Nd, u_smr.Grid()); std::vector<LF> V(Nd, u_smr.Grid());
for (int mu = 0; mu < Nd; mu++) { for (int mu = 0; mu < Nd; mu++) {
U[mu] = PeekIndex<LorentzIndex>(u_thin, mu); U[mu] = PeekIndex<LorentzIndex>(u_thin, mu);
V[mu] = PeekIndex<LorentzIndex>(u_smr, mu); V[mu] = PeekIndex<LorentzIndex>(u_smr, mu);
@ -317,7 +318,7 @@ public:
/*! @brief create long links from link variables. */ /*! @brief create long links from link variables. */
template<class LGF> template<class GF>
class Smear_HISQ_Naik { class Smear_HISQ_Naik {
private: private:
@ -333,7 +334,7 @@ public:
~Smear_HISQ_Naik() {} ~Smear_HISQ_Naik() {}
// void smear(LGF& u_smr, const LGF& U) const { // void smear(GF& u_smr, const GF& U) const {
// }; // };
// void derivative(const GaugeField& Gauge) const { // void derivative(const GaugeField& Gauge) const {

10
tests/smearing/Test_fatLinks.cc
View File

@ -67,6 +67,8 @@ int main (int argc, char** argv) {
std::string conf_out = "nersc.l8t4b3360.357link"; std::string conf_out = "nersc.l8t4b3360.357link";
int threads = GridThread::GetThreads(); int threads = GridThread::GetThreads();
typedef LatticeGaugeFieldD LGF;
// Initialize the Grid // Initialize the Grid
Grid_init(&argc,&argv); Grid_init(&argc,&argv);
Coordinate simd_layout = GridDefaultSimd(Nd,vComplexD::Nsimd()); Coordinate simd_layout = GridDefaultSimd(Nd,vComplexD::Nsimd());
@ -81,9 +83,9 @@ int main (int argc, char** argv) {
ConfParameters param(Reader); ConfParameters param(Reader);
if(param.benchmark) Grid_log(" Nloop = ",param.Nloop); if(param.benchmark) Grid_log(" Nloop = ",param.Nloop);
// Instantiate the LatticeGaugeField objects holding thin (Umu) and fat (U_smr) links // Instantiate the LGF objects holding thin (Umu) and fat (U_smr) links
LatticeGaugeField Umu(&GRID); LGF Umu(&GRID);
LatticeGaugeField U_smr(&GRID); LGF U_smr(&GRID);
// Read the configuration into Umu // Read the configuration into Umu
FieldMetaData header; FieldMetaData header;
@ -101,7 +103,7 @@ int main (int argc, char** argv) {
// Make sure result doesn't change w.r.t. a trusted lattice // Make sure result doesn't change w.r.t. a trusted lattice
NerscIO::readConfiguration(Umu, header, "nersc.l8t4b3360.357link.control"); NerscIO::readConfiguration(Umu, header, "nersc.l8t4b3360.357link.control");
LatticeGaugeField diff(&GRID); LGF diff(&GRID);
diff = Umu-U_smr; diff = Umu-U_smr;
auto absDiff = norm2(diff)/norm2(Umu); auto absDiff = norm2(diff)/norm2(Umu);
Grid_log(" |Umu-U|/|Umu| = ",absDiff); Grid_log(" |Umu-U|/|Umu| = ",absDiff);