fopr_Wilson_General.cpp

Go to the documentation of this file.

#include "fopr_Wilson_General.h"
#include "fopr_Wilson_General_impl.h"

#ifdef USE_FACTORY
namespace {
  Fopr *create_object()
  {
    return new Fopr_Wilson_General();
  }


  Fopr *create_object_with_repr(const std::string& repr)
  {
    return new Fopr_Wilson_General(repr);
  }


  bool init1 = Fopr::Factory_noarg::Register("Wilson_General", create_object);
  bool init2 = Fopr::Factory_string::Register("Wilson_General", create_object_with_repr);
}
#endif

//- parameter entries
namespace {
  void append_entry(Parameters& param)
  {
    param.Register_string("gamma_matrix_type", "NULL");
    param.Register_double("hopping_parameter_spatial", 0.0);
    param.Register_double("hopping_parameter_temporal", 0.0);
    param.Register_double("dispersion_parameter_spatial", 0.0);
    param.Register_double("Wilson_parameter_spatial", 0.0);
    param.Register_int_vector("boundary_condition", std::vector<int>());

    param.Register_string("verbose_level", "NULL");
  }


#ifdef USE_PARAMETERS_FACTORY
  bool init_param = ParametersFactory::Register("Fopr.Wilson_General", append_entry);
#endif
}
//- end

//- parameters class
Parameters_Fopr_Wilson_General::Parameters_Fopr_Wilson_General() { append_entry(*this); }
//- end

const std::string Fopr_Wilson_General::class_name = "Fopr_Wilson_General";

//====================================================================
void Fopr_Wilson_General::set_parameters(const Parameters& params)
{
  const string str_vlevel = params.get_string("verbose_level");

  m_vl = vout.set_verbose_level(str_vlevel);

  //- fetch and check input parameters
  double           kappa_s, kappa_t;
  double           nu_s, r_s;
  std::vector<int> bc;

  int err = 0;
  err += params.fetch_double("hopping_parameter_spatial", kappa_s);
  err += params.fetch_double("hopping_parameter_temporal", kappa_t);
  err += params.fetch_double("dispersion_parameter_spatial", nu_s);
  err += params.fetch_double("Wilson_parameter_spatial", r_s);
  err += params.fetch_int_vector("boundary_condition", bc);

  if (err) {
    vout.crucial(m_vl, "%s: fetch error, input parameter not found.\n", class_name.c_str());
    exit(EXIT_FAILURE);
  }

  set_parameters(kappa_s, kappa_t, nu_s, r_s, bc);
}


//====================================================================
void Fopr_Wilson_General::set_parameters(const double kappa_s, const double kappa_t,
                                         const double nu_s, const double r_s,
                                         const std::vector<int> bc)
{ m_impl->set_parameters(kappa_s, kappa_t, nu_s, r_s, bc); }

//====================================================================
Fopr_Wilson_General::Fopr_Wilson_General() : Fopr(), m_impl(new Fopr_Wilson_General_impl)
{ m_impl->init("Dirac"); }

//====================================================================
Fopr_Wilson_General::Fopr_Wilson_General(std::string repr) : Fopr(), m_impl(new Fopr_Wilson_General_impl)
{ m_impl->init(repr); }

//====================================================================
Fopr_Wilson_General::~Fopr_Wilson_General() { delete m_impl; }


//====================================================================
void Fopr_Wilson_General::set_config(Field *U)
{ m_impl->set_config(U); }

//====================================================================
void Fopr_Wilson_General::set_mode(std::string mode)
{ m_impl->set_mode(mode); }

//====================================================================
std::string Fopr_Wilson_General::get_mode() const
{ return m_impl->get_mode(); }


//====================================================================
void Fopr_Wilson_General::mult(Field& v, const Field& f)
{ m_impl->mult(v, f); }

//====================================================================
void Fopr_Wilson_General::mult_dag(Field& v, const Field& f)
{ m_impl->mult_dag(v, f); }

//====================================================================
void Fopr_Wilson_General::proj_chiral(Field& w, const int ex1,
                                      const Field& v, const int ex2, const int ipm)
{ m_impl->proj_chiral(w, ex1, v, ex2, ipm); }

//====================================================================
void Fopr_Wilson_General::mult_gm5(Field& w, const Field& v)
{ m_impl->mult_gm5(w, v); }


//====================================================================
void Fopr_Wilson_General::D(Field& w, const Field& v)
{ m_impl->D(w, v); }

//====================================================================
void Fopr_Wilson_General::Ddag(Field& w, const Field& v)
{ m_impl->Ddag(w, v); }

//====================================================================
void Fopr_Wilson_General::DdagD(Field& w, const Field& v)
{ m_impl->DdagD(w, v); }

//====================================================================
void Fopr_Wilson_General::DDdag(Field& w, const Field& v)
{ m_impl->DDdag(w, v); }

//====================================================================
void Fopr_Wilson_General::H(Field& w, const Field& v)
{ m_impl->H(w, v); }

//====================================================================
void Fopr_Wilson_General::mult_gm5p(int mu, Field_F& v, const Field_F& w)
{ m_impl->mult_gm5p(mu, v, w); }

//====================================================================
void Fopr_Wilson_General::D_ex(Field& w, const int ex1,
                               const Field& v, const int ex2)
{ m_impl->D_ex(w, ex1, v, ex2); }

//====================================================================
void Fopr_Wilson_General::mult_up(int mu, Field& w, const Field& v)
{ m_impl->mult_up(mu, w, v); }

//====================================================================
void Fopr_Wilson_General::mult_dn(int mu, Field& w, const Field& v)
{ m_impl->mult_dn(mu, w, v); }

//====================================================================
double Fopr_Wilson_General::flop_count()
{ return m_impl->flop_count(); }

//====================================================================
//============================================================END=====