Bridge++  Version 1.5.4
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
test_HMC_Wilson_Nf2.cpp
Go to the documentation of this file.
1 
14 #include "test.h"
15 
17 
18 #include "Fopr/fopr_Smeared.h"
19 
22 
23 #include "HMC/hmc_General.h"
24 #include "HMC/builder_Integrator.h"
25 
26 #include "IO/gaugeConfig.h"
27 
29 
30 #include "Tools/file_utils.h"
33 
34 //====================================================================
36 
50 namespace Test_HMC_Wilson {
51  const std::string test_name = "HMC.Wilson.Nf2";
52 
53  //- test-private parameters
54  namespace {
55  const std::string filename_input = "test_HMC_Wilson_Nf2.yaml";
56  }
57 
58  //- prototype declaration
59  int update_Nf2(void);
60 
61 #ifdef USE_TESTMANAGER_AUTOREGISTER
62  namespace {
63 #if defined(USE_GROUP_SU2)
64  // Nc=2 is not available.
65 #else
66  static const bool is_registered = TestManager::RegisterTest(
67  test_name,
69  );
70 #endif
71  }
72 #endif
73 
74  //====================================================================
75  int update_Nf2(void)
76  {
77  // ##### parameter setup #####
78  const int Nc = CommonParameters::Nc();
79  const int Nvol = CommonParameters::Nvol();
80  const int Ndim = CommonParameters::Ndim();
81 
82  const Parameters params_all = ParameterManager::read(filename_input);
83 
84  const Parameters params_test = params_all.lookup("Test_HMC_Wilson");
85  const Parameters params_action_G = params_all.lookup("Action_G");
86  const Parameters params_fopr = params_all.lookup("Fopr");
87  const Parameters params_proj = params_all.lookup("Projection");
88  const Parameters params_smear = params_all.lookup("Smear");
89  const Parameters params_dr_smear = params_all.lookup("Director_Smear");
90  const Parameters params_solver_MD = params_all.lookup("Solver_MD");
91  const Parameters params_solver_H = params_all.lookup("Solver_H");
92  const Parameters params_integrator = params_all.lookup("Builder_Integrator");
93  const Parameters params_hmc = params_all.lookup("HMC_General");
94 
95  const string str_gconf_status = params_test.get_string("gauge_config_status");
96  const string str_gconf_read = params_test.get_string("gauge_config_type_input");
97  const string readfile = params_test.get_string("config_filename_input");
98  const string str_gconf_write = params_test.get_string("gauge_config_type_output");
99  const string writefile = params_test.get_string("config_filename_output");
100  const string str_rand_type = params_test.get_string("random_number_type");
101  const unsigned long seed = params_test.get_unsigned_long("seed_for_random_number");
102  int i_conf = params_test.get_int("trajectory_number");
103  const int Ntraj = params_test.get_int("trajectory_number_step");
104  const int i_save_conf = params_test.get_int("save_config_interval");
105  const string str_vlevel = params_test.get_string("verbose_level");
106 
107  const bool do_check = params_test.is_set("expected_result");
108  const double expected_result = do_check ? params_test.get_double("expected_result") : 0.0;
109 
110  const string str_action_G_type = params_action_G.get_string("action_type");
111  const string str_fopr_type = params_fopr.get_string("fermion_type");
112  const string str_gmset_type = params_fopr.get_string("gamma_matrix_type");
113  const string str_proj_type = params_proj.get_string("projection_type");
114  const string str_smear_type = params_smear.get_string("smear_type");
115  const string str_solver_MD_type = params_solver_MD.get_string("solver_type");
116  const string str_solver_H_type = params_solver_H.get_string("solver_type");
117  const int Nlevels = params_integrator.get_int("number_of_levels");
118  const std::vector<int> level_action = params_integrator.get_int_vector("level_of_actions");
119 
120  const Bridge::VerboseLevel vl = vout.set_verbose_level(str_vlevel);
121 
122  //- print input parameters
123  vout.general(vl, " gconf_status = %s\n", str_gconf_status.c_str());
124  vout.general(vl, " gconf_read = %s\n", str_gconf_read.c_str());
125  vout.general(vl, " readfile = %s\n", readfile.c_str());
126  vout.general(vl, " gconf_write = %s\n", str_gconf_write.c_str());
127  vout.general(vl, " writefile = %s\n", writefile.c_str());
128  vout.general(vl, " rand_type = %s\n", str_rand_type.c_str());
129  vout.general(vl, " seed = %lu\n", seed);
130  vout.general(vl, " i_conf = %d\n", i_conf);
131  vout.general(vl, " Ntraj = %d\n", Ntraj);
132  vout.general(vl, " i_save_conf = %d\n", i_save_conf);
133  vout.general(vl, " vlevel = %s\n", str_vlevel.c_str());
134  vout.general(vl, " gmset_type = %s\n", str_gmset_type.c_str());
135  vout.general(vl, " proj_type = %s\n", str_proj_type.c_str());
136  vout.general(vl, " smear_type = %s\n", str_smear_type.c_str());
137  vout.general(vl, " solver_MD_type = %s\n", str_solver_MD_type.c_str());
138  vout.general(vl, " solver_H_type = %s\n", str_solver_H_type.c_str());
139  vout.general(vl, "\n");
140 
141  //- input parameter check
142  int err = 0;
143  err += ParameterCheck::non_NULL(str_gconf_status);
144  err += ParameterCheck::non_negative(i_conf);
145  err += ParameterCheck::non_negative(Ntraj);
146  err += ParameterCheck::non_negative(i_save_conf);
147 
148  if (err) {
149  vout.crucial(vl, "Error at %s: input parameters have not been set\n", test_name.c_str());
150  exit(EXIT_FAILURE);
151  }
152 
153 
154  RandomNumberManager::initialize(str_rand_type, seed);
155 
156 
157  // ##### object setup #####
158  unique_ptr<Field_G> U(new Field_G(Nvol, Ndim));
159 
160  if (str_gconf_status == "Continue") {
161  GaugeConfig(str_gconf_read).read(U, readfile);
162  } else if (str_gconf_status == "Cold_start") {
163  GaugeConfig("Unit").read(U);
164  } else if (str_gconf_status == "Hot_start") {
165  GaugeConfig("Random").read(U);
166  } else {
167  vout.crucial(vl, "Error at %s: unsupported gconf status \"%s\"\n", test_name.c_str(), str_gconf_status.c_str());
168  exit(EXIT_FAILURE);
169  }
170 
171  const unique_ptr<GaugeConfig> gconf_write(new GaugeConfig(str_gconf_write));
172 
173 
174  unique_ptr<Action> action_G(Action::New(str_action_G_type));
175  action_G->set_parameters(params_action_G);
176 
177  //-- N_f=2 part
178  unique_ptr<Fopr> fopr(Fopr::New(str_fopr_type, str_gmset_type));
179  fopr->set_parameters(params_fopr);
180 
181  unique_ptr<Force> force_fopr(new Force_F_Wilson_Nf2(str_gmset_type));
182  // define fermion force (SA)
183  force_fopr->set_parameters(params_fopr);
184 
185  // define smearing method (SA)
186  unique_ptr<Projection> proj(Projection::New(str_proj_type));
187  proj->set_parameters(params_proj);
188 
189  unique_ptr<Smear> smear(Smear::New(str_smear_type, proj));
190  smear->set_parameters(params_smear);
191 
192  // define force smearing method (SA)
193  unique_ptr<ForceSmear> force_smear(ForceSmear::New(str_smear_type, proj));
194  force_smear->set_parameters(params_smear);
195 
196  unique_ptr<Director> dr_smear(new Director_Smear(smear));
197  dr_smear->set_parameters(params_dr_smear);
198 
199  unique_ptr<Fopr> fopr_smear(Fopr::New("Smeared", fopr, dr_smear));
200  // define smeared fermion operator (SA)
201  unique_ptr<Force> force_fopr_smear(
202  new Force_F_Smeared(force_fopr, force_smear, dr_smear));
203  // define smeared fermion force (SA)
204 
205  unique_ptr<Solver> solver_MD(Solver::New(str_solver_MD_type, fopr_smear));
206  solver_MD->set_parameters(params_solver_MD);
207  unique_ptr<Fprop> fprop_MD(new Fprop_Standard_lex(solver_MD));
208 
209  unique_ptr<Solver> solver_H(Solver::New(str_solver_H_type, fopr_smear));
210  solver_H->set_parameters(params_solver_H);
211  unique_ptr<Fprop> fprop_H(new Fprop_Standard_lex(solver_H));
212 
213  unique_ptr<Action> action_F(
214  new Action_F_Standard_lex(fopr_smear, force_fopr_smear,
215  fprop_MD, fprop_H));
216  // define fermion action (SA)
217 
218 
219  ActionList actions(Nlevels);
220  actions.append(level_action[0], action_F);
221  actions.append(level_action[1], action_G);
222 
223  std::vector<Director *> directors(1);
224  directors[0] = (Director *)dr_smear.get(); // register director[0] (SA)
225 
226  const unique_ptr<Builder_Integrator> builder(new Builder_Integrator(actions, directors));
227  builder->set_parameters(params_integrator);
228  Integrator *integrator = builder->build();
229 
230  //- Random number is initialized with a parameter specified by i_conf
232 
233  HMC_General hmc(actions, directors, integrator, rand); // define hmc_leapfrog (SA)
234  hmc.set_parameters(params_hmc);
235 
236  const unique_ptr<Timer> timer(new Timer(test_name));
237 
238 
239  // #### Execution main part ####
240  timer->start();
241 
242  vout.general(vl, "HMC: Ntraj = %d\n", Ntraj); // a number of trajectory (SA)
243 
244  double result = 0.0;
245  for (int traj = 0; traj < Ntraj; ++traj) {
246  vout.general(vl, "\n");
247  vout.general(vl, "traj = %d\n", traj);
248 
249  result = hmc.update(*U); // hmc update (SA)
250 
251  if ((i_conf + traj + 1) % i_save_conf == 0) {
252  std::string filename = FileUtils::generate_filename("%s-%06d", writefile.c_str(), (i_conf + traj + 1));
253  gconf_write->write_file(U, filename);
254  }
255  }
256 
257  gconf_write->write_file(U, writefile);
258 
259  timer->report();
260 
262 
263 
264  if (do_check) {
265  return Test::verify(result, expected_result);
266  } else {
267  vout.detailed(vl, "check skipped: expected_result not set.\n\n");
268  return EXIT_SKIP;
269  }
270  }
271 } // namespace Test_HMC_Wilson
#define EXIT_SKIP
Definition: test.h:17
Random number generator base on M-series.
BridgeIO vout
Definition: bridgeIO.cpp:503
void detailed(const char *format,...)
Definition: bridgeIO.cpp:216
virtual void set_parameters(const Parameters &param)=0
void set_parameters(const Parameters &params)
void general(const char *format,...)
Definition: bridgeIO.cpp:197
General HMC update class.
Definition: hmc_General.h:45
virtual void set_parameters(const Parameters &)=0
int get_int(const string &key) const
Definition: parameters.cpp:192
Class for parameters.
Definition: parameters.h:46
bool append(const int level, Action *action)
Definition: action_list.cpp:27
virtual void set_parameters(const Parameters &param)=0
Base class of Integrator class family.
Definition: integrator.h:29
virtual void set_parameters(const Parameters &params)=0
Parameters lookup(const string &key) const
Definition: parameters.h:79
static bool RegisterTest(const std::string &key, const Test_function func)
Definition: testManager.h:69
static bool initialize(const std::string &rng_type, unsigned long seed)
SU(N) gauge field.
Definition: field_G.h:38
void set_parameters(const Parameters &params)
Manager of commonly used data object in HMC.
Definition: director.h:37
void read(Field_G *U, const string &filename=string())
unsigned long get_unsigned_long(const string &key) const
Definition: parameters.cpp:209
double get_double(const string &key) const
Definition: parameters.cpp:175
std::string generate_filename(const char *fmt,...)
Definition: file_utils.cpp:17
pointer get() const
virtual void set_parameters(const Parameters &)=0
int non_NULL(const std::string v)
bool is_set(const string &key) const
Definition: parameters.cpp:528
void start()
Definition: timer.cpp:44
Get quark propagator for Fopr with lexical site index.
void crucial(const char *format,...)
Definition: bridgeIO.cpp:178
virtual void set_parameters(const Parameters &params)=0
static void read(const std::string &params_file, Parameters &params)
Force for the standard Wilson fermion operator.
lists of actions at respective integrator levels.
Definition: action_list.h:40
Manager of smeared configurations.
int verify(const double result, const double expected, double eps)
Definition: test.cpp:27
Bridge::VerboseLevel vl
VerboseLevel
Definition: bridgeIO.h:42
double update(Field_G &)
void write_file(Field_G *U, const string &filename)
Definition: gaugeConfig.h:109
Force calculation for smeared fermion operators.
int non_negative(const int v)
virtual void set_parameters(const Parameters &)=0
GaugeConfig class for file I/O of gauge configuration.
Definition: gaugeConfig.h:78
Definition: timer.h:31
string get_string(const string &key) const
Definition: parameters.cpp:221
vector< int > get_int_vector(const string &key) const
Definition: parameters.cpp:267
void report(const Bridge::VerboseLevel vl=Bridge::GENERAL)
Definition: timer.cpp:128
virtual void set_parameters(const Parameters &)=0
Builder of MD integrator for HMC.
Standard fermion action for HMC.
static VerboseLevel set_verbose_level(const std::string &str)
Definition: bridgeIO.cpp:131