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