41 namespace Test_Rational {
42 const std::string
test_name =
"Rational.Smeared_Rational";
46 const std::string filename_input =
"test_Rational_Smeared.yaml";
52 #ifdef USE_TESTMANAGER_AUTOREGISTER
54 #if defined(USE_GROUP_SU2)
73 int NinF = 2 * Nc * Nd;
85 const string str_gconf_status = params_test.
get_string(
"gauge_config_status");
86 const string str_gconf_read = params_test.
get_string(
"gauge_config_type_input");
87 const string readfile = params_test.
get_string(
"config_filename_input");
88 const string str_rand_type = params_test.
get_string(
"random_number_type");
89 const unsigned long seed = params_test.
get_unsigned_long(
"seed_for_random_number");
90 const string str_vlevel = params_test.
get_string(
"verbose_level");
92 const bool do_check = params_test.
is_set(
"expected_result");
93 const double expected_result = do_check ? params_test.
get_double(
"expected_result") : 0.0;
95 const string str_gmset_type = params_clover.
get_string(
"gamma_matrix_type");
96 const string str_proj_type = params_proj.
get_string(
"projection_type");
97 const string str_smear_type = params_smear.
get_string(
"smear_type");
98 const string str_source_type = params_source.
get_string(
"source_type");
103 vout.
general(vl,
" gconf_status = %s\n", str_gconf_status.c_str());
104 vout.
general(vl,
" gconf_read = %s\n", str_gconf_read.c_str());
105 vout.
general(vl,
" readfile = %s\n", readfile.c_str());
106 vout.
general(vl,
" rand_type = %s\n", str_rand_type.c_str());
108 vout.
general(vl,
" vlevel = %s\n", str_vlevel.c_str());
109 vout.
general(vl,
" gmset_type = %s\n", str_gmset_type.c_str());
110 vout.
general(vl,
" proj_type = %s\n", str_proj_type.c_str());
111 vout.
general(vl,
" smear_type = %s\n", str_smear_type.c_str());
112 vout.
general(vl,
" source_type = %s\n", str_source_type.c_str());
131 if (str_gconf_status ==
"Continue") {
133 }
else if (str_gconf_status ==
"Cold_start") {
135 }
else if (str_gconf_status ==
"Hot_start") {
138 vout.
crucial(vl,
"Error at %s: unsupported gconf status \"%s\"\n",
test_name.c_str(), str_gconf_status.c_str());
179 Field xq(NinF, Nvol, 1), b(NinF, Nvol, 1);
180 Field v(NinF, Nvol, 1), w(NinF, Nvol, 1);
183 vout.
general(vl,
"Consistency check of smeared-rational and rational-smeared.\n");
194 int idx = icolor + Nc * ispin;
195 source->
set(b2, idx);
199 fopr_smear_2.
mult(v, b);
202 fopr_smear_1.
mult(w, v);
204 fopr_smear_2.
mult(v, w);
207 double vv = v.norm2();
223 int idx = icolor + Nc * ispin;
224 source->
set(b2, idx);
228 fopr_r_smear_1.
mult(v, b);
231 fopr_smear_1.
mult(w, v);
233 fopr_r_smear_1.
mult(v, w);
252 vout.
detailed(vl,
"check skipped: expected_result not set.\n\n");
Fermion operator for rational approximation.
const std::string test_name
void detailed(const char *format,...)
int smeared_rational(void)
virtual void set_parameters(const Parameters ¶m)=0
void general(const char *format,...)
Container of Field-type object.
virtual void set_parameters(const Parameters &)=0
void mult(Field &v, const Field &f)
multiplies fermion operator to a given field (2nd argument)
Wilson-type fermion field.
smeared fermion operator.
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
void set_config(Field *U)
setting pointer to the gauge configuration.
virtual void set_parameters(const Parameters &)=0
int non_NULL(const std::string v)
bool is_set(const string &key) const
virtual void set(Field &, int)=0
void axpy(Field &y, const double a, const Field &x)
axpy(y, a, x): y := a * x + y
void crucial(const char *format,...)
virtual void set_parameters(const Parameters ¶ms)=0
static void read(const std::string ¶ms_file, Parameters ¶ms)
void mult(Field &v, const Field &f)
multiply smeared fermion operator
Manager of smeared configurations.
int verify(const double result, const double expected, double eps)
void set_mode(std::string mode)
set the mode of fermion operator
void set_parameters(const Parameters ¶ms)
GaugeConfig class for file I/O of gauge configuration.
Base class of fermion operator family.
string get_string(const string &key) const
void report(const Bridge::VerboseLevel vl=Bridge::GENERAL)
void set_config(Field *U)
set pointer to original thin link variable
virtual void set_parameters(const Parameters &)=0
static VerboseLevel set_verbose_level(const std::string &str)