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

Adding Scalar field theory example for Scidac format

Browse Source
This commit is contained in:
Guido Cossu 2018-05-02 14:36:32 +01:00
parent e568c24d1d
commit 4d4ac2517b
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
tests/hmc/Test_hmc_ScalarActionNxN.cc
View File

@ -34,6 +34,8 @@ class ScalarActionParameters : Serializable {
double, lambda, double, lambda,
double, g); double, g);
ScalarActionParameters() = default;
template <class ReaderClass > template <class ReaderClass >
ScalarActionParameters(Reader<ReaderClass>& Reader){ ScalarActionParameters(Reader<ReaderClass>& Reader){
read(Reader, "ScalarAction", *this); read(Reader, "ScalarAction", *this);
@ -124,10 +126,13 @@ int main(int argc, char **argv) {
ScalarGrid.set_rb(new GridRedBlackCartesian(ScalarGrid.get_full())); ScalarGrid.set_rb(new GridRedBlackCartesian(ScalarGrid.get_full()));
TheHMC.Resources.AddGrid("scalar", ScalarGrid); TheHMC.Resources.AddGrid("scalar", ScalarGrid);
std::cout << "Lattice size : " << GridDefaultLatt() << std::endl; std::cout << "Lattice size : " << GridDefaultLatt() << std::endl;
ScalarActionParameters SPar(Reader);
// Checkpointer definition // Checkpointer definition
CheckpointerParameters CPparams(Reader); CheckpointerParameters CPparams(Reader);
TheHMC.Resources.LoadBinaryCheckpointer(CPparams); //TheHMC.Resources.LoadBinaryCheckpointer(CPparams);
TheHMC.Resources.LoadScidacCheckpointer(CPparams, SPar);
RNGModuleParameters RNGpar(Reader); RNGModuleParameters RNGpar(Reader);
TheHMC.Resources.SetRNGSeeds(RNGpar); TheHMC.Resources.SetRNGSeeds(RNGpar);
@ -140,7 +145,6 @@ int main(int argc, char **argv) {
// Collect actions, here use more encapsulation // Collect actions, here use more encapsulation
// Scalar action in adjoint representation // Scalar action in adjoint representation
ScalarActionParameters SPar(Reader);
ScalarAction Saction(SPar.mass_squared, SPar.lambda, SPar.g); ScalarAction Saction(SPar.mass_squared, SPar.lambda, SPar.g);
// Collect actions // Collect actions