portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Exposed a coulomb/landau enum to the gauge fixing module

Browse Source
This commit is contained in:
fionnoh 2019-05-17 19:01:52 +01:00
parent 48b1c806ed
commit d14512ee03
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Hadrons/Modules/MGauge/GaugeFix.hpp
View File

@ -42,6 +42,8 @@ BEGIN_HADRONS_NAMESPACE
******************************************************************************/ ******************************************************************************/
BEGIN_MODULE_NAMESPACE(MGauge) BEGIN_MODULE_NAMESPACE(MGauge)
GRID_SERIALIZABLE_ENUM(Fix, undef, coulomb, Nd - 1, landau, -1);
class GaugeFixPar: Serializable class GaugeFixPar: Serializable
{ {
public: public:
@ -51,6 +53,7 @@ public:
int, maxiter, int, maxiter,
Real, Omega_tol, Real, Omega_tol,
Real, Phi_tol, Real, Phi_tol,
Fix, gaugeFix,
bool, Fourier); bool, Fourier);
}; };
@ -115,8 +118,8 @@ void TGaugeFix<GImpl>::execute(void)
//Loads the gauge and fixes it //Loads the gauge and fixes it
{ {
std::cout << "executing" << std::endl; std::cout << "executing" << std::endl;
LOG(Message) << "Fixing the Gauge" << std::endl; LOG(Message) << "Fixing the Gauge " << par().gauge << " using "
LOG(Message) << par().gauge << std::endl; << par().gaugeFix << " guage fixing. " << Nd - 1 << std::endl;
auto &U = envGet(GaugeField, par().gauge); auto &U = envGet(GaugeField, par().gauge);
auto &Umu = envGet(GaugeField, getName()); auto &Umu = envGet(GaugeField, getName());
auto &xform = envGet(GaugeMat, getName()+"_xform"); auto &xform = envGet(GaugeMat, getName()+"_xform");
@ -126,9 +129,10 @@ void TGaugeFix<GImpl>::execute(void)
int maxiter = par().maxiter; int maxiter = par().maxiter;
Real Omega_tol = par().Omega_tol; Real Omega_tol = par().Omega_tol;
Real Phi_tol = par().Phi_tol; Real Phi_tol = par().Phi_tol;
int gaugeFix = par().gaugeFix;
bool Fourier = par().Fourier; bool Fourier = par().Fourier;
Umu = U; Umu = U;
FourierAcceleratedGaugeFixer<PeriodicGimplR>::SteepestDescentGaugeFix(Umu,xform,alpha,maxiter,Omega_tol,Phi_tol,Fourier); FourierAcceleratedGaugeFixer<PeriodicGimplR>::SteepestDescentGaugeFix(Umu,xform,alpha,maxiter,Omega_tol,Phi_tol,Fourier,gaugeFix);
LOG(Message) << "Gauge Fixed" << std::endl; LOG(Message) << "Gauge Fixed" << std::endl;
} }