33 namespace Test_GradientFlow {
44 const std::string filename_input_RK1 =
"test_GradientFlow_Nf0_RK1.yaml";
45 const std::string filename_input_RK2 =
"test_GradientFlow_Nf0_RK2.yaml";
46 const std::string filename_input_RK3 =
"test_GradientFlow_Nf0_RK3.yaml";
47 const std::string filename_input_RK4 =
"test_GradientFlow_Nf0_RK4.yaml";
48 const std::string filename_input_RK_adaptive =
"test_GradientFlow_Nf0_RK_adaptive.yaml";
52 int update(
const std::string& filename_input);
55 {
return update(filename_input_RK1); }
58 {
return update(filename_input_RK2); }
61 {
return update(filename_input_RK3); }
64 {
return update(filename_input_RK4); }
67 {
return update(filename_input_RK_adaptive); }
69 #ifdef USE_TESTMANAGER_AUTOREGISTER
71 #if defined(USE_GROUP_SU2)
99 int update(
const std::string& filename_input)
113 const string str_gconf_status = params_test.
get_string(
"gauge_config_status");
114 const string str_gconf_read = params_test.
get_string(
"gauge_config_type_input");
115 const string readfile = params_test.
get_string(
"config_filename_input");
116 const string str_gconf_write = params_test.
get_string(
"gauge_config_type_output");
117 const string writefile = params_test.
get_string(
"config_filename_output");
118 const string str_rand_type = params_test.
get_string(
"random_number_type");
119 const unsigned long seed = params_test.
get_unsigned_long(
"seed_for_random_number");
120 const int Nstep = params_test.
get_int(
"number_of_steps");
121 const double t_flow_max = params_test.
get_double(
"max_flow_time");
122 const string str_vlevel = params_test.
get_string(
"verbose_level");
124 const bool do_check = params_test.
is_set(
"expected_result");
125 const double expected_result = do_check ? params_test.
get_double(
"expected_result") : 0.0;
127 const string str_action_G_type = params_action_G.
get_string(
"action_type");
132 vout.
general(vl,
" gconf_status = %s\n", str_gconf_status.c_str());
133 vout.
general(vl,
" gconf_read = %s\n", str_gconf_read.c_str());
134 vout.
general(vl,
" readfile = %s\n", readfile.c_str());
135 vout.
general(vl,
" gconf_write = %s\n", str_gconf_write.c_str());
136 vout.
general(vl,
" writefile = %s\n", writefile.c_str());
137 vout.
general(vl,
" rand_type = %s\n", str_rand_type.c_str());
140 vout.
general(vl,
" t_flow_max = %10.6f\n", t_flow_max);
141 vout.
general(vl,
" vlevel = %s\n", str_vlevel.c_str());
160 if (str_gconf_status ==
"Continue") {
162 }
else if (str_gconf_status ==
"Cold_start") {
164 }
else if (str_gconf_status ==
"Hot_start") {
167 vout.
crucial(vl,
"Error at %s: unsupported gconf status \"%s\"\n",
test_name.c_str(), str_gconf_status.c_str());
173 params_action_G.
set_double(
"beta", static_cast<double>(Nc));
191 for (
int i = 0; i < Nstep; ++i) {
192 result = g_flow->
evolve(t_flow, *U);
194 double t2 = t_flow * t_flow;
195 double E_plaq = energy_density->
E_plaq(*U);
196 double E_clover = energy_density->
E_clover(*U);
198 vout.
general(vl,
" (t, t^2 E_plaq, t^2 E_clover) = %.8f %.16f %.16f\n",
199 t_flow, t2 * E_plaq, t2 * E_clover);
202 if (t_flow > t_flow_max)
break;
213 vout.
detailed(vl,
"check skipped: expected_result not set.\n\n");
void detailed(const char *format,...)
double E_plaq(const Field_G &U)
void general(const char *format,...)
const std::string test_name
double evolve(double &t, Field_G &U)
double E_clover(const Field_G &U)
int get_int(const string &key) const
const std::string test_name_RK1
const std::string test_name_RK_adaptive
virtual void set_parameters(const Parameters ¶m)=0
virtual void set_parameters(const Parameters ¶ms)
setting parameters.
int update(const std::string &filename_input)
const std::string test_name_RK2
Parameters lookup(const string &key) const
static bool RegisterTest(const std::string &key, const Test_function func)
static bool initialize(const std::string &rng_type, unsigned long seed)
void read(Field_G *U, const string &filename=string())
unsigned long get_unsigned_long(const string &key) const
double get_double(const string &key) const
energy density measurement.
const std::string test_name_RK3
int non_NULL(const std::string v)
bool is_set(const string &key) const
void crucial(const char *format,...)
static void read(const std::string ¶ms_file, Parameters ¶ms)
GradientFlow construction.
int verify(const double result, const double expected, double eps)
void set_double(const string &key, const double value)
const std::string test_name_RK4
GaugeConfig class for file I/O of gauge configuration.
string get_string(const string &key) const
int run_test_RK_adaptive()
void report(const Bridge::VerboseLevel vl=Bridge::GENERAL)
void set_parameters(const Parameters ¶ms)
static VerboseLevel set_verbose_level(const std::string &str)