26 const string str_vlevel = params.
get_string(
"verbose_level");
30 double c_plaq, c_rect;
67 double E_plaq = 36.0 * (1.0 - plaq);
71 std::ofstream log_file;
98 double E_clover_1x1 = 0.0;
99 double E_clover_1x2 = 0.0;
104 std::vector<Field_G> Fmunu_1x1(6);
108 for (
int mu = 0; mu < Ndim; ++mu) {
109 for (
int nu = mu + 1; nu < Ndim; ++nu) {
117 for (
int site = 0; site < Nvol; ++site) {
118 for (
int i_munu = 0; i_munu < 6; i_munu++) {
119 F2_1x1 +=
ReTr(Fmunu_1x1[i_munu].mat(site) * Fmunu_1x1[i_munu].mat(site));
130 std::vector<Field_G> Fmunu_1x2(6);
133 for (
int mu = 0; mu < Ndim; ++mu) {
134 for (
int nu = mu + 1; nu < Ndim; ++nu) {
142 for (
int site = 0; site < Nvol; ++site) {
143 for (
int i_munu = 0; i_munu < 6; i_munu++) {
144 F2_1x2 +=
ReTr(Fmunu_1x2[i_munu].mat(site) * Fmunu_1x2[i_munu].mat(site));
160 std::ofstream log_file;
void construct_Fmunu_1x2(Field_G &Fmunu, const int mu, const int nu, const Field_G &U)
FieldStrength m_field_strength
std::string m_filename_output
static double epsilon_criterion()
double E_plaq(const Field_G &U)
void general(const char *format,...)
int fetch_double(const string &key, double &value) const
double E_clover(const Field_G &U)
double plaquette(const Field_G &)
calculates plaquette value.
void init(const std::string &filename)
virtual void set_parameters(const Parameters ¶ms)
setting parameters.
void construct_Fmunu_1x1(Field_G &Fmunu, const int mu, const int nu, const Field_G &U)
Bridge::VerboseLevel m_vl
void crucial(const char *format,...)
std::ostream & getStream()
static int reduce_sum(int count, double *recv_buf, double *send_buf, int pattern=0)
make a global sum of an array of double over the communicator. pattern specifies the dimensions to be...
static const std::string class_name
string get_string(const string &key) const
static VerboseLevel set_verbose_level(const std::string &str)
double ReTr(const Mat_SU_N &m)