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

Some comments in the hmc files

Browse Source
This commit is contained in:
Guido Cossu 2016-11-10 17:55:58 +00:00
parent 1d666771f9
commit 19b85d8486
1 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
tests/hmc/Test_hmc_WilsonGauge_Binary.cc
View File

@ -58,6 +58,7 @@ class HMCRunnerParameters : Serializable {
HMCRunnerParameters() {} HMCRunnerParameters() {}
}; };
// Derive from the BinaryHmcRunner (templated for gauge fields) // Derive from the BinaryHmcRunner (templated for gauge fields)
class HmcRunner : public BinaryHmcRunner { class HmcRunner : public BinaryHmcRunner {
public: public:
@ -65,15 +66,21 @@ class HmcRunner : public BinaryHmcRunner {
void BuildTheAction(int argc, char **argv) void BuildTheAction(int argc, char **argv)
{ {
// Typedefs to simplify notation
typedef WilsonGaugeActionR GaugeAction;
typedef WilsonImplR ImplPolicy; typedef WilsonImplR ImplPolicy;
typedef WilsonFermionR FermionAction; typedef WilsonFermionR FermionAction;
typedef typename FermionAction::FermionField FermionField; typedef typename FermionAction::FermionField FermionField;
UGrid = SpaceTimeGrid::makeFourDimGrid(GridDefaultLatt(), GridDefaultSimd(Nd, vComplex::Nsimd()), GridDefaultMpi()); // this can be simplified too. MakeDefaultGrid(Nd)
UGrid = SpaceTimeGrid::makeFourDimGrid(
GridDefaultLatt(),
GridDefaultSimd(Nd, vComplex::Nsimd()),
GridDefaultMpi());
// Gauge action // Gauge action
WilsonGaugeActionR Waction(HMCPar.beta); GaugeAction Waction(HMCPar.beta);
// Collect actions // Collect actions
ActionLevel<Field> Level1(1); ActionLevel<Field> Level1(1);
@ -113,12 +120,14 @@ int main(int argc, char **argv) {
<< " threads" << std::endl; << " threads" << std::endl;
HmcRunner TheHMC; HmcRunner TheHMC;
// make input file name general
InputFileReader Reader("input.wilson_gauge.params"); InputFileReader Reader("input.wilson_gauge.params");
read(Reader, "HMC", TheHMC.HMCPar); read(Reader, "HMC", TheHMC.HMCPar);
std::cout << GridLogMessage << TheHMC.HMCPar << std::endl; std::cout << GridLogMessage << TheHMC.HMCPar << std::endl;
// Seeds for the random number generators // Seeds for the random number generators
// generalise
std::vector<int> SerSeed = strToVec<int>(TheHMC.HMCPar.serial_seeds); std::vector<int> SerSeed = strToVec<int>(TheHMC.HMCPar.serial_seeds);
std::vector<int> ParSeed = strToVec<int>(TheHMC.HMCPar.parallel_seeds); std::vector<int> ParSeed = strToVec<int>(TheHMC.HMCPar.parallel_seeds);