From dda6c69d5b32aa7d108b225928b76b8e6748b265 Mon Sep 17 00:00:00 2001 From: Antonin Portelli Date: Mon, 5 Mar 2018 19:58:40 +0000 Subject: [PATCH] Hadrons: scalar SU(N) shift probes --- extras/Hadrons/Modules.hpp | 1 + .../Hadrons/Modules/MScalarSUN/ShiftProbe.hpp | 169 ++++++++++++++++++ extras/Hadrons/modules.inc | 1 + 3 files changed, 171 insertions(+) create mode 100644 extras/Hadrons/Modules/MScalarSUN/ShiftProbe.hpp diff --git a/extras/Hadrons/Modules.hpp b/extras/Hadrons/Modules.hpp index 695dfd02..ee53faa8 100644 --- a/extras/Hadrons/Modules.hpp +++ b/extras/Hadrons/Modules.hpp @@ -61,6 +61,7 @@ See the full license in the file "LICENSE" in the top level distribution directo #include #include #include +#include #include #include #include diff --git a/extras/Hadrons/Modules/MScalarSUN/ShiftProbe.hpp b/extras/Hadrons/Modules/MScalarSUN/ShiftProbe.hpp new file mode 100644 index 00000000..89059180 --- /dev/null +++ b/extras/Hadrons/Modules/MScalarSUN/ShiftProbe.hpp @@ -0,0 +1,169 @@ +/************************************************************************************* + +Grid physics library, www.github.com/paboyle/Grid + +Source file: extras/Hadrons/Modules/MScalarSUN/ShiftProbe.hpp + +Copyright (C) 2015-2018 + +Author: Antonin Portelli + +This program is free software; you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation; either version 2 of the License, or +(at your option) any later version. + +This program is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + +You should have received a copy of the GNU General Public License along +with this program; if not, write to the Free Software Foundation, Inc., +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 */ +#ifndef Hadrons_MScalarSUN_ShiftProbe_hpp_ +#define Hadrons_MScalarSUN_ShiftProbe_hpp_ + +#include +#include +#include +#include + +BEGIN_HADRONS_NAMESPACE + +/****************************************************************************** + * Ward identity phi^n probe with fields at different positions * + ******************************************************************************/ +BEGIN_MODULE_NAMESPACE(MScalarSUN) + +typedef std::pair ShiftPair; + +class ShiftProbePar: Serializable +{ +public: + GRID_SERIALIZABLE_CLASS_MEMBERS(ShiftProbePar, + std::string, field, + std::string, shifts, + std::string, output); +}; + +template +class TShiftProbe: public Module +{ +public: + + typedef typename SImpl::Field Field; + typedef typename SImpl::ComplexField ComplexField; + class Result: Serializable + { + public: + GRID_SERIALIZABLE_CLASS_MEMBERS(Result, + std::string, op, + Complex , value); + }; +public: + // constructor + TShiftProbe(const std::string name); + // destructor + virtual ~TShiftProbe(void) = default; + // dependency relation + virtual std::vector getInput(void); + virtual std::vector getOutput(void); + // setup + virtual void setup(void); + // execution + virtual void execute(void); +}; + +MODULE_REGISTER_NS(ShiftProbeSU2, TShiftProbe>, MScalarSUN); +MODULE_REGISTER_NS(ShiftProbeSU3, TShiftProbe>, MScalarSUN); +MODULE_REGISTER_NS(ShiftProbeSU4, TShiftProbe>, MScalarSUN); +MODULE_REGISTER_NS(ShiftProbeSU5, TShiftProbe>, MScalarSUN); +MODULE_REGISTER_NS(ShiftProbeSU6, TShiftProbe>, MScalarSUN); + +/****************************************************************************** + * TShiftProbe implementation * + ******************************************************************************/ +// constructor ///////////////////////////////////////////////////////////////// +template +TShiftProbe::TShiftProbe(const std::string name) +: Module(name) +{} + +// dependencies/products /////////////////////////////////////////////////////// +template +std::vector TShiftProbe::getInput(void) +{ + std::vector in = {par().field}; + + return in; +} + +template +std::vector TShiftProbe::getOutput(void) +{ + std::vector out = {getName()}; + + return out; +} + +// setup /////////////////////////////////////////////////////////////////////// +template +void TShiftProbe::setup(void) +{ + envTmpLat(Field, "acc"); + envCacheLat(ComplexField, getName()); +} + +// execution /////////////////////////////////////////////////////////////////// +template +void TShiftProbe::execute(void) +{ + LOG(Message) << "Creating shift probe for shifts " << par().shifts + << std::endl; + + std::vector shift; + unsigned int sign; + auto &phi = envGet(Field, par().field); + auto &probe = envGet(ComplexField, getName()); + + shift = strToVec(par().shifts); + if (shift.size() % 2 != 0) + { + HADRON_ERROR(Size, "the number of shifts is odd"); + } + sign = (shift.size() % 4 == 0) ? 1 : -1; + for (auto &s: shift) + { + if (s.first >= env().getNd()) + { + HADRON_ERROR(Size, "dimension to large for shift <" + + std::to_string(s.first) + " " + + std::to_string(s.second) + ">" ); + } + } + envGetTmp(Field, acc); + acc = 1.; + for (unsigned int i = 0; i < shift.size(); ++i) + { + if (shift[i].second == 0) + { + acc *= phi; + } + else + { + acc *= Cshift(phi, shift[i].first, shift[i].second); + } + } + probe = sign*trace(acc); +} + +END_MODULE_NAMESPACE + +END_HADRONS_NAMESPACE + +#endif // Hadrons_MScalarSUN_ShiftProbe_hpp_ diff --git a/extras/Hadrons/modules.inc b/extras/Hadrons/modules.inc index 44df3091..c0be8aef 100644 --- a/extras/Hadrons/modules.inc +++ b/extras/Hadrons/modules.inc @@ -44,6 +44,7 @@ modules_hpp =\ Modules/MAction/Wilson.hpp \ Modules/MAction/WilsonClover.hpp \ Modules/MAction/ZMobiusDWF.hpp \ + Modules/MScalarSUN/ShiftProbe.hpp \ Modules/MScalarSUN/Div.hpp \ Modules/MScalarSUN/TrMag.hpp \ Modules/MScalarSUN/EMT.hpp \