60 }
else if (
m_mode ==
"Ddag") {
63 }
else if (
m_mode ==
"DdagD") {
66 }
else if (
m_mode ==
"DDdag") {
69 }
else if (
m_mode ==
"H") {
89 const string str_vlevel = params.
get_string(
"verbose_level");
116 for (
int mu = 0; mu <
m_Ndim; ++mu) {
122 assert(bc.size() ==
m_Ndim);
128 for (
int mu = 0; mu <
m_Ndim; ++mu) {
141 for (
int mu = 0; mu <
m_Ndim; ++mu) {
159 for (
int mu = 0; mu <
m_Ndim; ++mu) {
174 assert(v.
nex() == f.
nex());
175 assert(v.
nin() == f.
nin());
187 assert(ipm == 1 || ipm == -1);
196 }
else if (ipm == -1) {
242 for (
int ex = 0; ex < f.
nex(); ++ex) {
257 for (
int ex = 0; ex < f.
nex(); ++ex) {
278 int flop_per_site = Nc * Nd * 2 * 8 * (4 * Nc - 1);
280 flop_per_site += Nc * Nd * 2 * (4 * 3 + 4 * 2);
281 flop_per_site += Nc * Nd * 2 * 8;
282 flop_per_site += Nc * Nd * 2 * 2;
284 double flop =
static_cast<double>(flop_per_site) *
285 static_cast<double>(Lvol);
void scal(Field &x, const double a)
scal(x, a): x = a * x
void mult_Field_Gd(Field_F &y, const int ex, const Field_G &u, int ex1, const Field_F &x, int ex2)
void set(const int jin, const int site, const int jex, double v)
static const std::string class_name
Bridge::VerboseLevel m_vl
void general(const char *format,...)
static Bridge::VerboseLevel Vlevel()
const Field_F mult_gm5p(int mu, const Field_F &w)
Container of Field-type object.
int fetch_double(const string &key, double &value) const
void proj_chiral(Field &w, const int ex1, const Field &v, const int ex2, const int ipm)
void mult_gm5(Field &v, const Field &f)
gamma_5 multiplication. [31 Mar 2017 H.Matsufuru]
void H(Field &w, const Field &f)
void addpart_ex(int ex, const Field &w, int exw)
std::string get_mode() const
only for Fopr_Overlap
Wilson-type fermion field.
void mult_GMproj2(Field_F &y, const int pm, const GammaMatrix &gm, const Field_F &x)
projection with gamma matrix: (1 gamma)
void set_parameters(const Parameters ¶ms)
std::vector< int > m_boundary
void(Fopr_Wilson::* m_mult)(Field &, const Field &)
void set_mode(std::string mode)
setting the mode of multiplication if necessary. Default implementation here is just to avoid irrelev...
void D(Field &v, const Field &f)
void mult_up(int mu, Field &w, const Field &f)
nearest neighbor hopping term: temporary entry [H.Matsufuru]
void init(std::string repr)
void backward(Field &, const Field &, const int mu)
void mult_Field_Gn(Field_F &y, const int ex, const Field_G &u, int ex1, const Field_F &x, int ex2)
void DdagD(Field &w, const Field &f)
void axpy(Field &y, const double a, const Field &x)
axpy(y, a, x): y := a * x + y
std::vector< GammaMatrix > m_GM
void crucial(const char *format,...)
void mult_GM(Field_F &y, const GammaMatrix &gm, const Field_F &x)
gamma matrix multiplication
void DDdag(Field &w, const Field &f)
void(Fopr_Wilson::* m_mult_dag)(Field &, const Field &)
void mult_dn(int mu, Field &w, const Field &f)
void setpart_ex(int ex, const Field &w, int exw)
string get_string(const string &key) const
int fetch_int_vector(const string &key, vector< int > &value) const
double flop_count()
this returns the number of floating point operations.
void D_ex(Field &v, const int ex1, const Field &f, const int ex2)
void Ddag(Field &w, const Field &f)
static VerboseLevel set_verbose_level(const std::string &str)
void mult_undef(Field &, const Field &f)
void forward(Field &, const Field &, const int mu)