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

rm unused variables and formatting

Browse Source
This commit is contained in:
Alessandro Lupo 2023-06-23 16:04:18 +01:00
parent f27d2083cd
commit cff1f8d3b8
6 changed files with 51 additions and 69 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
Grid/qcd/QCD.h
View File

@ -47,7 +47,6 @@ static constexpr int Ym = 5;
static constexpr int Zm = 6; static constexpr int Zm = 6;
static constexpr int Tm = 7; static constexpr int Tm = 7;
static constexpr int isSymplectic=Sp2n_config;
static constexpr int Nc=Config_Nc; static constexpr int Nc=Config_Nc;
static constexpr int Ns=4; static constexpr int Ns=4;
static constexpr int Nd=4; static constexpr int Nd=4;

10
configure.ac
View File

@ -146,7 +146,6 @@ case ${ac_TRACING} in
AC_DEFINE([GRID_TRACING_TIMER],[1],[use TIMER]);; AC_DEFINE([GRID_TRACING_TIMER],[1],[use TIMER]);;
*) *)
AC_DEFINE([GRID_TRACING_NONE],[1],[no tracing]);; AC_DEFINE([GRID_TRACING_NONE],[1],[no tracing]);;
esac esac
############### fermions ############### fermions
@ -185,18 +184,15 @@ AC_ARG_ENABLE([Nc],
case ${ac_Nc} in case ${ac_Nc} in
2) 2)
AC_DEFINE([Config_Nc],[2],[Gauge group Nc]) AC_DEFINE([Config_Nc],[2],[Gauge group Nc]);;
AC_DEFINE([Config_nsp],[1],[n of Sp2n]);;
3) 3)
AC_DEFINE([Config_Nc],[3],[Gauge group Nc]);; AC_DEFINE([Config_Nc],[3],[Gauge group Nc]);;
4) 4)
AC_DEFINE([Config_Nc],[4],[Gauge group Nc]) AC_DEFINE([Config_Nc],[4],[Gauge group Nc]);;
AC_DEFINE([Config_nsp],[2],[n of Sp2n]);;
5) 5)
AC_DEFINE([Config_Nc],[5],[Gauge group Nc]);; AC_DEFINE([Config_Nc],[5],[Gauge group Nc]);;
8) 8)
AC_DEFINE([Config_Nc],[8],[Gauge group Nc]) AC_DEFINE([Config_Nc],[8],[Gauge group Nc]);;
AC_DEFINE([Config_nsp],[4],[n of Sp2n]);;
*) *)
AC_MSG_ERROR(["Unsupport gauge group choice Nc = ${ac_Nc}"]);; AC_MSG_ERROR(["Unsupport gauge group choice Nc = ${ac_Nc}"]);;
esac esac

77
tests/sp2n/Test_hmc_Sp_WF_2_Fund_3_2AS.cc
View File

@ -2,34 +2,34 @@
int main(int argc, char **argv) { int main(int argc, char **argv) {
using namespace Grid; using namespace Grid;
typedef Representations< SpFundamentalRepresentation, SpTwoIndexAntiSymmetricRepresentation > TheRepresentations; typedef Representations< SpFundamentalRepresentation, SpTwoIndexAntiSymmetricRepresentation > TheRepresentations;
Grid_init(&argc, &argv); Grid_init(&argc, &argv);
typedef GenericSpHMCRunnerHirep<TheRepresentations, MinimumNorm2> HMCWrapper; typedef GenericSpHMCRunnerHirep<TheRepresentations, MinimumNorm2> HMCWrapper;
typedef SpWilsonTwoIndexAntiSymmetricImplR TwoIndexFermionImplPolicy; typedef SpWilsonTwoIndexAntiSymmetricImplR TwoIndexFermionImplPolicy;
typedef SpWilsonTwoIndexAntiSymmetricFermionD TwoIndexFermionAction; typedef SpWilsonTwoIndexAntiSymmetricFermionD TwoIndexFermionAction;
typedef typename TwoIndexFermionAction::FermionField TwoIndexFermionField; typedef typename TwoIndexFermionAction::FermionField TwoIndexFermionField;
typedef SpWilsonImplR FundFermionImplPolicy; // ok typedef SpWilsonImplR FundFermionImplPolicy; // ok
typedef SpWilsonFermionD FundFermionAction; // ok typedef SpWilsonFermionD FundFermionAction; // ok
typedef typename FundFermionAction::FermionField FundFermionField; typedef typename FundFermionAction::FermionField FundFermionField;
//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: //::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
HMCWrapper TheHMC; HMCWrapper TheHMC;
TheHMC.Resources.AddFourDimGrid("gauge"); TheHMC.Resources.AddFourDimGrid("gauge");
// Checkpointer definition // Checkpointer definition
CheckpointerParameters CPparams; CheckpointerParameters CPparams;
CPparams.config_prefix = "ckpoint_lat"; CPparams.config_prefix = "ckpoint_lat";
CPparams.rng_prefix = "ckpoint_rng"; CPparams.rng_prefix = "ckpoint_rng";
CPparams.saveInterval = 5; CPparams.saveInterval = 5;
CPparams.format = "IEEE64BIG"; CPparams.format = "IEEE64BIG";
TheHMC.Resources.LoadNerscCheckpointer(CPparams); TheHMC.Resources.LoadNerscCheckpointer(CPparams);
RNGModuleParameters RNGpar; RNGModuleParameters RNGpar;
@ -40,63 +40,58 @@ int main(int argc, char **argv) {
// Construct observables // Construct observables
typedef PlaquetteMod<HMCWrapper::ImplPolicy> PlaqObs; typedef PlaquetteMod<HMCWrapper::ImplPolicy> PlaqObs;
TheHMC.Resources.AddObservable<PlaqObs>(); TheHMC.Resources.AddObservable<PlaqObs>();
typedef PolyakovMod<HMCWrapper::ImplPolicy> PolyakovObs;
TheHMC.Resources.AddObservable<PolyakovObs>();
typedef PolyakovMod<HMCWrapper::ImplPolicy> PolyakovObs;
TheHMC.Resources.AddObservable<PolyakovObs>();
RealD beta = 6 ; RealD beta = 6 ;
SpWilsonGaugeActionR Waction(beta); SpWilsonGaugeActionR Waction(beta);
auto GridPtr = TheHMC.Resources.GetCartesian(); auto GridPtr = TheHMC.Resources.GetCartesian();
auto GridRBPtr = TheHMC.Resources.GetRBCartesian(); auto GridRBPtr = TheHMC.Resources.GetRBCartesian();
SpFundamentalRepresentation::LatticeField fundU(GridPtr); SpFundamentalRepresentation::LatticeField fundU(GridPtr);
SpTwoIndexAntiSymmetricRepresentation::LatticeField asU(GridPtr); SpTwoIndexAntiSymmetricRepresentation::LatticeField asU(GridPtr);
//LatticeGaugeField U(GridPtr); //LatticeGaugeField U(GridPtr);
RealD Fundmass = -0.71; RealD Fundmass = -0.71;
RealD ASmass = -0.71; RealD ASmass = -0.71;
std::vector<Complex> boundary = {-1,-1,-1,-1}; std::vector<Complex> boundary = {-1,-1,-1,-1};
FundFermionAction::ImplParams bc(boundary); FundFermionAction::ImplParams bc(boundary);
TwoIndexFermionAction::ImplParams bbc(boundary); TwoIndexFermionAction::ImplParams bbc(boundary);
FundFermionAction FundFermOp(fundU, *GridPtr, *GridRBPtr, Fundmass, bbc); FundFermionAction FundFermOp(fundU, *GridPtr, *GridRBPtr, Fundmass, bbc);
TwoIndexFermionAction TwoIndexFermOp(asU, *GridPtr, *GridRBPtr, ASmass, bbc); TwoIndexFermionAction TwoIndexFermOp(asU, *GridPtr, *GridRBPtr, ASmass, bbc);
ConjugateGradient<FundFermionField> fCG(1.0e-8, 2000, false); ConjugateGradient<FundFermionField> fCG(1.0e-8, 2000, false);
ConjugateGradient<TwoIndexFermionField> asCG(1.0e-8, 2000, false); ConjugateGradient<TwoIndexFermionField> asCG(1.0e-8, 2000, false);
OneFlavourRationalParams Params(1.0e-6, 64.0, 2000, 1.0e-6, 16); OneFlavourRationalParams Params(1.0e-6, 64.0, 2000, 1.0e-6, 16);
TwoFlavourPseudoFermionAction<FundFermionImplPolicy> fundNf2(FundFermOp, fCG, fCG); TwoFlavourPseudoFermionAction<FundFermionImplPolicy> fundNf2(FundFermOp, fCG, fCG);
TwoFlavourPseudoFermionAction<TwoIndexFermionImplPolicy> asNf2(TwoIndexFermOp, asCG, asCG); TwoFlavourPseudoFermionAction<TwoIndexFermionImplPolicy> asNf2(TwoIndexFermOp, asCG, asCG);
OneFlavourRationalPseudoFermionAction<TwoIndexFermionImplPolicy> asNf1(TwoIndexFermOp,Params); OneFlavourRationalPseudoFermionAction<TwoIndexFermionImplPolicy> asNf1(TwoIndexFermOp,Params);
fundNf2.is_smeared = false; fundNf2.is_smeared = false;
asNf2.is_smeared = false; asNf2.is_smeared = false;
asNf1.is_smeared = false; asNf1.is_smeared = false;
ActionLevel<HMCWrapper::Field, TheRepresentations > Level1(1); ActionLevel<HMCWrapper::Field, TheRepresentations > Level1(1);
Level1.push_back(&fundNf2); Level1.push_back(&fundNf2);
Level1.push_back(&asNf2); Level1.push_back(&asNf2);
Level1.push_back(&asNf1); Level1.push_back(&asNf1);
ActionLevel<HMCWrapper::Field, TheRepresentations > Level2(4); ActionLevel<HMCWrapper::Field, TheRepresentations > Level2(4);
Level2.push_back(&Waction); Level2.push_back(&Waction);
TheHMC.TheAction.push_back(Level1); TheHMC.TheAction.push_back(Level1);
TheHMC.TheAction.push_back(Level2); TheHMC.TheAction.push_back(Level2);
TheHMC.Parameters.MD.MDsteps = 28; TheHMC.Parameters.MD.MDsteps = 28;
TheHMC.Parameters.MD.trajL = 1.0; TheHMC.Parameters.MD.trajL = 1.0;
TheHMC.ReadCommandLine(argc, argv); TheHMC.ReadCommandLine(argc, argv);
TheHMC.Run(); TheHMC.Run();
Grid_finalize(); Grid_finalize();
} }

1
tests/sp2n/Test_hmc_Sp_Wilson2ASFermionGauge.cc
View File

@ -79,4 +79,3 @@ int main(int argc, char **argv) {
Grid_finalize(); Grid_finalize();
} }

28
tests/sp2n/Test_hmc_Sp_WilsonFundFermionGauge.cc
View File

@ -43,36 +43,32 @@ int main(int argc, char **argv) {
auto GridPtr = TheHMC.Resources.GetCartesian(); auto GridPtr = TheHMC.Resources.GetCartesian();
auto GridRBPtr = TheHMC.Resources.GetRBCartesian(); auto GridRBPtr = TheHMC.Resources.GetRBCartesian();
SpFundamentalRepresentation::LatticeField U(GridPtr); SpFundamentalRepresentation::LatticeField U(GridPtr);
//LatticeGaugeField U(GridPtr);
RealD mass = -0.76; RealD mass = -0.76;
FermionAction FermOp(U, *GridPtr, *GridRBPtr, mass); FermionAction FermOp(U, *GridPtr, *GridRBPtr, mass);
ConjugateGradient<FermionField> CG(1.0e-8, 2000, false); ConjugateGradient<FermionField> CG(1.0e-8, 2000, false);
TwoFlavourPseudoFermionAction<FermionImplPolicy> Nf2(FermOp, CG, CG); TwoFlavourPseudoFermionAction<FermionImplPolicy> Nf2(FermOp, CG, CG);
Nf2.is_smeared = false; Nf2.is_smeared = false;
ActionLevel<HMCWrapper::Field, TheRepresentations > Level1(1); ActionLevel<HMCWrapper::Field, TheRepresentations > Level1(1);
Level1.push_back(&Nf2); Level1.push_back(&Nf2);
ActionLevel<HMCWrapper::Field, TheRepresentations > Level2(4); ActionLevel<HMCWrapper::Field, TheRepresentations > Level2(4);
Level2.push_back(&Waction); Level2.push_back(&Waction);
TheHMC.TheAction.push_back(Level1); TheHMC.TheAction.push_back(Level1);
TheHMC.TheAction.push_back(Level2); TheHMC.TheAction.push_back(Level2);
TheHMC.Parameters.MD.MDsteps = 36; TheHMC.Parameters.MD.MDsteps = 36;
TheHMC.Parameters.MD.trajL = 1.0; TheHMC.Parameters.MD.trajL = 1.0;
TheHMC.ReadCommandLine(argc, argv);
TheHMC.Run();
TheHMC.ReadCommandLine(argc, argv);
TheHMC.Run();
Grid_finalize(); Grid_finalize();
} }

3
tests/sp2n/Test_hmc_Sp_pureGaugeWilson.cc
View File

@ -40,11 +40,8 @@ int main(int argc, char **argv)
typedef GenericSpHMCRunner<MinimumNorm2> HMCWrapper; typedef GenericSpHMCRunner<MinimumNorm2> HMCWrapper;
HMCWrapper TheHMC; HMCWrapper TheHMC;
TheHMC.Resources.AddFourDimGrid("gauge"); TheHMC.Resources.AddFourDimGrid("gauge");
// Checkpointer definition // Checkpointer definition
CheckpointerParameters CPparams; CheckpointerParameters CPparams;
CPparams.config_prefix = "ckpoint_lat"; CPparams.config_prefix = "ckpoint_lat";