Bridge++  Ver. 1.1.x
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
test_HMC_Clover_Leapfrog_Nf2_eo.cpp
Go to the documentation of this file.
1 
14 #include "parameterManager_YAML.h"
15 #include "parameters_factory.h"
16 
17 #include "bridgeIO.h"
18 using Bridge::vout;
19 
20 #include "gaugeConfig.h"
21 
22 #include "action_G_Plaq.h"
23 #include "action_G_Rectangle.h"
24 
25 #include "fopr_Clover.h"
26 #include "force_F_Clover_Nf2.h"
27 
28 #include "randomNumbers_Mseries.h"
29 
30 #include "director_Smear.h"
31 #include "fopr_Smeared.h"
32 #include "force_F_Smeared.h"
33 
34 #include "forceSmear.h"
35 #include "projection.h"
36 #include "smear.h"
37 #include "solver.h"
38 
39 #include "action_F_Standard_eo.h"
40 #include "fopr_Clover_eo.h"
41 #include "fopr_Smeared_eo.h"
42 #include "fprop_Standard_eo.h"
43 
44 #include "hmc_Leapfrog.h"
45 
46 #ifdef USE_TEST
47 #include "test.h"
48 #endif
49 
50 #ifdef USE_TESTMANAGER_AUTOREGISTER
51 #include "testManager.h"
52 #endif
53 
54 //====================================================================
56 
71 namespace Test_HMC_Clover {
72  //- test-private parameters
73  namespace {
74  const std::string filename_input = "test_HMC_Clover_Leapfrog_Nf2.yaml";
75  const std::string filename_output = "stdout";
76 
77  class Parameters_Test_HMC_Clover : public Parameters {
78  public:
79  Parameters_Test_HMC_Clover()
80  {
81  Register_string("gauge_config_status", "NULL");
82 
83  Register_string("gauge_config_type_input", "NULL");
84  Register_string("config_filename_input", "NULL");
85 
86  Register_string("gauge_config_type_output", "NULL");
87  Register_string("config_filename_output", "NULL");
88 
89  Register_int("trajectory_number", 0);
90  Register_int("trajectory_number_step", 0);
91  Register_int("save_config_interval", 0);
92 
93  Register_string("verbose_level", "NULL");
94 
95  Register_double("expected_result", 0.0);
96  }
97  };
98  }
99 
100  //- prototype declaration
101  int leapfrog_Nf2_eo(void);
102 
103 #ifdef USE_TESTMANAGER_AUTOREGISTER
104  namespace {
105  static const bool is_registered = TestManager::RegisterTest(
106  "HMC.Clover.Leapfrog_Nf2_eo",
107  leapfrog_Nf2_eo
108  );
109  }
110 #endif
111 
112  //====================================================================
113  int leapfrog_Nf2_eo(void)
114  {
115  // ##### parameter setup #####
116  int Nc = CommonParameters::Nc();
117  int Nvol = CommonParameters::Nvol();
118  int Ndim = CommonParameters::Ndim();
119  int NinG = 2 * Nc * Nc;
120 
121  Parameters *params_test = new Parameters_Test_HMC_Clover;
122  Parameters *params_action_G = ParametersFactory::New("Action.G_Rectangle");
123  Parameters *params_clover = ParametersFactory::New("Fopr.Clover");
124  Parameters *params_proj = ParametersFactory::New("Projection");
125  Parameters *params_smear = ParametersFactory::New("Smear");
126  Parameters *params_dr_smear = ParametersFactory::New("Director_Smear");
127  Parameters *params_solver_MD = ParametersFactory::New("Solver");
128  Parameters *params_solver_H = ParametersFactory::New("Solver");
129  Parameters *params_hmc = ParametersFactory::New("HMC.Leapfrog");
130 
131  Parameters *params_all = new Parameters;
132 
133  params_all->Register_Parameters("Test_HMC_Clover", params_test);
134  params_all->Register_Parameters("Action_G_Rectangle", params_action_G);
135  params_all->Register_Parameters("Fopr_Clover", params_clover);
136  params_all->Register_Parameters("Projection", params_proj);
137  params_all->Register_Parameters("Smear", params_smear);
138  params_all->Register_Parameters("Director_Smear", params_dr_smear);
139  params_all->Register_Parameters("Solver_MD", params_solver_MD);
140  params_all->Register_Parameters("Solver_H", params_solver_H);
141  params_all->Register_Parameters("HMC_Leapfrog", params_hmc);
142 
143  ParameterManager_YAML params_manager;
144  params_manager.read_params(filename_input, params_all);
145 
146  const string str_gconf_status = params_test->get_string("gauge_config_status");
147  const string str_gconf_read = params_test->get_string("gauge_config_type_input");
148  const string readfile = params_test->get_string("config_filename_input");
149  const string str_gconf_write = params_test->get_string("gauge_config_type_output");
150  const string writefile = params_test->get_string("config_filename_output");
151  int iconf = params_test->get_int("trajectory_number");
152  const int Ntraj = params_test->get_int("trajectory_number_step");
153  const int i_save_conf = params_test->get_int("save_config_interval");
154  const string str_vlevel = params_test->get_string("verbose_level");
155 #ifdef USE_TEST
156  const double expected_result = params_test->get_double("expected_result");
157 #endif
158 
159  const string str_gmset_type = params_clover->get_string("gamma_matrix_type");
160  const string str_proj_type = params_proj->get_string("projection_type");
161  const string str_smear_type = params_smear->get_string("smear_type");
162  const string str_solver_MD_type = params_solver_MD->get_string("solver_type");
163  const string str_solver_H_type = params_solver_H->get_string("solver_type");
164 
165  Bridge::VerboseLevel vl = vout.set_verbose_level(str_vlevel);
166 
167  //- print input parameters
168  vout.general(vl, " gconf_status = %s\n", str_gconf_status.c_str());
169  vout.general(vl, " gconf_read = %s\n", str_gconf_read.c_str());
170  vout.general(vl, " readfile = %s\n", readfile.c_str());
171  vout.general(vl, " gconf_write = %s\n", str_gconf_write.c_str());
172  vout.general(vl, " writefile = %s\n", writefile.c_str());
173  vout.general(vl, " iconf = %d\n", iconf);
174  vout.general(vl, " Ntraj = %d\n", Ntraj);
175  vout.general(vl, " i_save_conf = %d\n", i_save_conf);
176  vout.general(vl, " vlevel = %s\n", str_vlevel.c_str());
177 
178  vout.general(vl, " gmset_type = %s\n", str_gmset_type.c_str());
179  vout.general(vl, " proj_type = %s\n", str_proj_type.c_str());
180  vout.general(vl, " smear_type = %s\n", str_smear_type.c_str());
181  vout.general(vl, " solver_MD_type = %s\n", str_solver_MD_type.c_str());
182  vout.general(vl, " solver_H_type = %s\n", str_solver_H_type.c_str());
183  vout.general(vl, "\n");
184 
185  //- input parameter check
186  int err = 0;
187  err += ParameterCheck::non_NULL(str_gconf_read);
188  err += ParameterCheck::non_NULL(readfile);
189  err += ParameterCheck::non_NULL(str_gconf_write);
190  err += ParameterCheck::non_NULL(writefile);
191  err += ParameterCheck::non_zero(iconf);
192  err += ParameterCheck::non_zero(Ntraj);
193  err += ParameterCheck::non_zero(i_save_conf);
194 
195  if (err) {
196  vout.crucial(vl, "test_HMC_Clover: Input parameters have not been set.\n");
197  abort();
198  }
199 
200 
201  // ##### object setup #####
202  Field_G *U = new Field_G(Nvol, Ndim);
203  GaugeConfig *gconf_read = new GaugeConfig(str_gconf_read);
204 
205  if (str_gconf_status == "Continue") {
206  gconf_read->read_file((Field *)U, readfile);
207  } else if (str_gconf_status == "Start_cold") {
208  gconf_read->set_cold((Field *)U);
209  } else {
210  vout.crucial(vl, "Test_HMC_Clover: unsupported gconf status \"%s\".\n", str_gconf_status.c_str());
211  abort();
212  }
213 
214 
215  Action_G_Rectangle *action_G = new Action_G_Rectangle; // plaq.+rectangle (SA)
216  action_G->set_parameters(*params_action_G);
217 
218 
219  Fopr_Clover *fopr_c = new Fopr_Clover(str_gmset_type); // define fermion operator (SA)
220  fopr_c->set_parameters(*params_clover);
221 
222  Force_F_Clover_Nf2 *force_fopr_c = new Force_F_Clover_Nf2(str_gmset_type);
223  // define fermion force (SA)
224  force_fopr_c->set_parameters(*params_clover);
225 
226  // define smearing method (SA)
227  Projection *proj = Projection::New(str_proj_type);
228  Smear *smear = Smear::New(str_smear_type, proj);
229  smear->set_parameters(*params_smear);
230 
231  // define force smearing method (SA)
232  ForceSmear *force_smear = ForceSmear::New(str_smear_type, proj);
233  force_smear->set_parameters(*params_smear);
234 
235  Director_Smear *dr_smear = new Director_Smear((Smear *)smear);
236  dr_smear->set_parameters(*params_dr_smear);
237 
238  Fopr_Smeared *fopr_smear = new Fopr_Smeared((Fopr *)fopr_c, dr_smear);
239  // define smeared fermion operator (SA)
240  Force_F_Smeared *force_fopr_smear
241  = new Force_F_Smeared((Force *)force_fopr_c, (ForceSmear *)force_smear, dr_smear);
242  // define smeared fermion force (SA)
243 
244  //- NB1 fopr_c_eo->set_config is performed in Action_F_Standard_eo.
245  // NB2 fopr_c_eo->set_mode is performed in fprop.
246  Fopr_Clover_eo *fopr_c_eo = new Fopr_Clover_eo(str_gmset_type);
247  fopr_c_eo->set_parameters(*params_clover);
248 
249  Fopr_Smeared_eo *fopr_smear_eo = new Fopr_Smeared_eo((Fopr_eo *)fopr_c_eo, dr_smear);
250 
251  Solver *solver_eo_MD = Solver::New(str_solver_MD_type, fopr_smear_eo);
252  solver_eo_MD->set_parameters(*params_solver_MD);
253  Fprop *fprop_eo_MD = new Fprop_Standard_eo(solver_eo_MD);
254 
255  Solver *solver_eo_H = Solver::New(str_solver_H_type, fopr_smear_eo);
256  solver_eo_H->set_parameters(*params_solver_H);
257  Fprop *fprop_eo_H = new Fprop_Standard_eo(solver_eo_H);
258 
259  Action_F_Standard_eo *action_F
260  = new Action_F_Standard_eo((Fopr *)fopr_smear, (Force *)force_fopr_smear,
261  fprop_eo_MD, fprop_eo_H);
262  // define fermion action (SA)
263 
264 
265  valarray<Action *> actions(2);
266  actions[0] = (Action *)action_G; // register action[0] (SA)
267  actions[1] = (Action *)action_F; // register action[1] (SA)
268 
269  valarray<Director *> directors(1);
270  directors[0] = (Director *)dr_smear; // register director[0] (SA)
271 
272 
273  //- Random number is initialized with a parameter specified by iconf
274  RandomNumbers *rand = new RandomNumbers_Mseries(iconf);
275 
276  HMC_Leapfrog hmc(actions, directors, rand); // define hmc_leapfrog (SA)
277  hmc.set_parameters(*params_hmc);
278 
279 
280  // #### Execution main part ####
281  vout.general(vl, "HMC: Ntraj = %d\n", Ntraj); // a number of trajectory (SA)
282 
283  double result = 0.0;
284  for (int traj = 0; traj < Ntraj; ++traj) {
285  vout.general(vl, "\n");
286  vout.general(vl, "traj = %d\n", traj);
287 
288  result = hmc.update(*U); // hmc update (SA)
289  }
290 
291 
292  // ##### tidy up #####
293  delete U;
294  delete gconf_read;
295 
296  delete rand;
297 
298  delete action_G;
299  delete action_F;
300 
301  delete fopr_c;
302  delete force_fopr_c;
303 
304  delete proj;
305  delete smear;
306  delete force_smear;
307  delete dr_smear;
308  delete fopr_smear;
309  delete force_fopr_smear;
310 
311  delete fopr_c_eo;
312  delete fopr_smear_eo;
313 
314  delete solver_eo_MD;
315  delete fprop_eo_MD;
316  delete solver_eo_H;
317  delete fprop_eo_H;
318 
319  delete params_test;
320  delete params_action_G;
321  delete params_clover;
322  delete params_proj;
323  delete params_smear;
324  delete params_dr_smear;
325  delete params_solver_MD;
326  delete params_solver_H;
327  delete params_hmc;
328  delete params_all;
329 
330 
331 #ifdef USE_TEST
332  return Test::verify(expected_result, result);
333 
334 #else
335  return EXIT_SUCCESS;
336 #endif
337  }
338 } // namespace Test_HMC_Clover