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