mirror of
https://github.com/paboyle/Grid.git
synced 2024-09-20 09:15:38 +01:00
92 lines
2.6 KiB
C++
92 lines
2.6 KiB
C++
#include <Grid/Grid.h>
|
|
|
|
int main(int argc, char** argv) {
|
|
using namespace Grid;
|
|
|
|
typedef Representations<SpFundamentalRepresentation,
|
|
SpTwoIndexAntiSymmetricRepresentation>
|
|
TheRepresentations;
|
|
|
|
Grid_init(&argc, &argv);
|
|
|
|
typedef GenericSp2nHMCRunnerHirep<TheRepresentations, MinimumNorm2>
|
|
HMCWrapper;
|
|
typedef SpWilsonTwoIndexAntiSymmetricImplR FermionImplPolicy;
|
|
typedef SpWilsonTwoIndexAntiSymmetricFermionR FermionAction;
|
|
typedef typename FermionAction::FermionField FermionField;
|
|
|
|
//::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
|
|
|
|
HMCWrapper TheHMC;
|
|
|
|
TheHMC.Resources.AddFourDimGrid("gauge");
|
|
|
|
// Checkpointer definition
|
|
CheckpointerParameters CPparams;
|
|
CPparams.config_prefix = "ckpoint_lat";
|
|
CPparams.rng_prefix = "ckpoint_rng";
|
|
CPparams.saveInterval = 100;
|
|
CPparams.format = "IEEE64BIG";
|
|
|
|
TheHMC.Resources.LoadNerscCheckpointer(CPparams);
|
|
|
|
RNGModuleParameters RNGpar;
|
|
RNGpar.serial_seeds = "1 2 3 4 5";
|
|
RNGpar.parallel_seeds = "6 7 8 9 10";
|
|
TheHMC.Resources.SetRNGSeeds(RNGpar);
|
|
|
|
// Construct observables
|
|
typedef PlaquetteMod<HMCWrapper::ImplPolicy> PlaqObs;
|
|
TheHMC.Resources.AddObservable<PlaqObs>();
|
|
|
|
RealD beta = 6.7;
|
|
|
|
SpWilsonGaugeActionR Waction(beta);
|
|
|
|
auto GridPtr = TheHMC.Resources.GetCartesian();
|
|
auto GridRBPtr = TheHMC.Resources.GetRBCartesian();
|
|
|
|
SpTwoIndexAntiSymmetricRepresentation::LatticeField U(GridPtr);
|
|
// LatticeGaugeField U(GridPtr);
|
|
|
|
RealD mass = -0.115;
|
|
|
|
std::vector<Complex> boundary = {-1, -1, -1, -1};
|
|
FermionAction::ImplParams bc(boundary);
|
|
FermionAction FermOp(U, *GridPtr, *GridRBPtr, mass, bc);
|
|
|
|
ConjugateGradient<FermionField> CG(1.0e-8, 2000, false);
|
|
|
|
TwoFlavourPseudoFermionAction<FermionImplPolicy> Nf2(FermOp, CG, CG);
|
|
|
|
Nf2.is_smeared = false;
|
|
std::cout << GridLogMessage << "mass " << mass << std::endl;
|
|
|
|
ActionLevel<HMCWrapper::Field, TheRepresentations> Level1(1);
|
|
Level1.push_back(&Nf2);
|
|
|
|
ActionLevel<HMCWrapper::Field, TheRepresentations> Level2(4);
|
|
Level2.push_back(&Waction);
|
|
|
|
TheHMC.TheAction.push_back(Level1);
|
|
TheHMC.TheAction.push_back(Level2);
|
|
|
|
TheHMC.Parameters.MD.MDsteps = 16;
|
|
TheHMC.Parameters.MD.trajL = 1.0;
|
|
|
|
TheHMC.ReadCommandLine(argc, argv);
|
|
TheHMC.Run();
|
|
|
|
if (std::string(argv[argc - 1]).compare("--assert=true") == 0) {
|
|
// this is the pseudofermion action's value for
|
|
// --grid 4.4.4.4 --Trajectories 1 --Thermalizations 0 --StartingType
|
|
// ColdStart and Nc == 4 if this assertion fails, check that you have run
|
|
// with these exact parameters
|
|
assert(Nc == 4);
|
|
assert(abs(Nf2.S(U) - 6586.6498) < 1.e-4);
|
|
}
|
|
|
|
Grid_finalize();
|
|
}
|
|
|