26 bool init = Source::Factory::Register(
"Wall", create_object);
37 const string str_vlevel = params.
get_string(
"verbose_level");
42 std::vector<int> source_position;
64 std::vector<int> Lsize(Ndim);
72 std::vector<int> Nsize(Ndim);
78 const int t_dir = Ndim - 1;
82 vout.
general(
m_vl,
" source_position[t] = %d\n", source_position[t_dir]);
96 for (
int mu = 0; mu < Ndim; ++mu) {
99 m_source_position[t_dir] = (source_position[t_dir] + Lsize[t_dir]) % Lsize[t_dir];
125 const int Lvol3 = Lx * Ly * Lz;
128 std::vector<int> Nsize(Ndim);
141 for (
int z = 0; z < Nsize[2]; ++z) {
142 for (
int y = 0; y < Nsize[1]; ++y) {
143 for (
int x = 0; x < Nsize[0]; ++x) {
147 src.
set(2 * j, isite, 0, 1.0 / Lvol3);
void set(const int jin, const int site, const int jex, double v)
int site(const int &x, const int &y, const int &z, const int &t) const
void general(const char *format,...)
void set(Field &v, int j)
static const std::string class_name
Container of Field-type object.
static int ipe(const int dir)
logical coordinate of current proc.
void set_parameters(const Parameters ¶ms)
void crucial(const char *format,...)
Bridge::VerboseLevel m_vl
int non_negative(const int v)
std::vector< int > m_source_position
Base class of source for a linear solver.
string get_string(const string &key) const
int fetch_int_vector(const string &key, vector< int > &value) const
Wall source for 4-spinor fermion.
static VerboseLevel set_verbose_level(const std::string &str)