fopr_Wilson_eo.h

Go to the documentation of this file.

#ifndef FOPR_WILSON_EO_INCLUDED
#define FOPR_WILSON_EO_INCLUDED

#include "fopr_Wilson.h"

#include "fopr_eo.h"
#include "shiftField_eo.h"

#include "gammaMatrixSet.h"

#include "bridgeIO.h"
using Bridge::vout;

//- parameters class
class Parameters_Fopr_Wilson_eo : virtual public Parameters
{
 public:
  Parameters_Fopr_Wilson_eo();
};
//- end


class Fopr_Wilson_eo : public Fopr_eo
{
 public:
  static const std::string class_name;

 private:
  int m_Nc, m_Nd, m_Ndim;
  int m_Nvol, m_Nvol2;

  double             m_kappa;
  std::valarray<int> m_boundary;
  Index_eo           m_index;
  Field_G            *m_Ueo;

  ShiftField_eo shift;
  Field_F       trf, trf2;
  Vec_SU_N      v1, v2;

  std::valarray<GammaMatrix> m_GM;

  std::string m_mode;
  std::string m_repr;

  //  Field_F m_vF1;     //!< working field.

  void (Fopr_Wilson_eo::*m_mult)(Field &, const Field &);
  void (Fopr_Wilson_eo::*m_mult_dag)(Field &, const Field &);
  void (Fopr_Wilson_eo::*m_D)(Field &, const Field &);
  void (Fopr_Wilson_eo::*m_gm5)(Field &, const Field &);
  void (Fopr_Wilson_eo::*m_preProp)(Field &, Field &, const Field &);
  void (Fopr_Wilson_eo::*m_postProp)(Field &, const Field &, const Field &);

 public:

  Fopr_Wilson_eo();

  Fopr_Wilson_eo(std::string repr);

  ~Fopr_Wilson_eo();

  void set_parameters(const Parameters& params);
  void set_parameters(const double kappa, const std::valarray<int> bc);

  void set_config(Field *U);

  // method for even odd fermion operator
  void preProp(Field& Be, Field& bo, const Field& b)
  { (this->*m_preProp)(Be, bo, b); }

  void postProp(Field& x, const Field& xe, const Field& bo)
  { (this->*m_postProp)(x, xe, bo); }

  void prePropD(Field&, Field&, const Field&);
  void postPropD(Field&, const Field&, const Field&);
  void prePropDag(Field&, Field&, const Field&);
  void postPropDag(Field&, const Field&, const Field&);

  const Field mult(const Field& f)
  {
    assert_single_thread();
    Field v(f.nin(), f.nvol(), f.nex());
    mult(v, f);
    return v;
  }

  const Field mult_dag(const Field& f)
  {
    assert_single_thread();
    Field v(f.nin(), f.nvol(), f.nex());
    mult_dag(v, f);
    return v;
  }

  void mult(Field& v, const Field& f)
  { (this->*m_mult)(v, f); }

  void mult_dag(Field& v, const Field& f)
  { (this->*m_mult_dag)(v, f); }

  void set_mode(std::string mode)
  {
    m_mode = mode;

    if (m_mode == "D") {
      m_mult     = &Fopr_Wilson_eo::D;
      m_mult_dag = &Fopr_Wilson_eo::Ddag;
      m_preProp  = &Fopr_Wilson_eo::prePropD;
      m_postProp = &Fopr_Wilson_eo::postPropD;
    } else if (m_mode == "Ddag") {
      m_mult     = &Fopr_Wilson_eo::Ddag;
      m_mult_dag = &Fopr_Wilson_eo::D;
      m_preProp  = &Fopr_Wilson_eo::prePropDag;
      m_postProp = &Fopr_Wilson_eo::postPropDag;
    } else if (m_mode == "DdagD") {
      m_mult     = &Fopr_Wilson_eo::DdagD;
      m_mult_dag = &Fopr_Wilson_eo::DdagD;
    } else if (m_mode == "DDdag") {
      m_mult     = &Fopr_Wilson_eo::DDdag;
      m_mult_dag = &Fopr_Wilson_eo::DDdag;
    } else if (m_mode == "H") {
      m_mult     = &Fopr_Wilson_eo::H;
      m_mult_dag = &Fopr_Wilson_eo::H;
    } else {
      vout.crucial("%s: input mode is undefined.\n", class_name.c_str());
      abort();
    }
  }

  std::string get_mode() const
  { return m_mode; }

  void assert_single_thread();

  double flop_count();

  void mult_undef(Field&, const Field&);
  void D(Field& v, const Field& f);
  void Ddag(Field& v, const Field& f);
  void DdagD(Field& v, const Field& f);
  void DDdag(Field& v, const Field& f);
  void H(Field& v, const Field& f);
  void mult_gm5(Field& v, const Field& f);

  const Field MeoMoe(const Field& f);
  const Field D(const Field& f);
  const Field Ddag(const Field& f);
  const Field DdagD(const Field& f);
  const Field DDdag(const Field& f);
  const Field H(const Field& f);
  const Field mult_gm5(const Field&);

  const Field_F Meo(const Field_F&, const int ieo);
  const Field_F Mdageo(const Field_F&, const int ieo);

  //  void Meo(Field_F&, const Field_F&, const int ieo);
  void Meo(Field&, const Field&, const int ieo);

  const Field_F Meo_gm5(const Field_F&, const int ieo);

  const Field gm5p(const int mu, const Field& v);

  int field_nvol() { return m_Nvol2; }
  int field_nin() { return 2 * m_Nc * m_Nd; }
  int field_nex() { return 1; }

 private:
  //  void init(const std::string);
  void mult_p(int mu, Field_F&, const Field_F&, const int ieo);
  void mult_m(int mu, Field_F&, const Field_F&, const int ieo);

  class Fopr_Wilson_eo_impl;

  Fopr_Wilson_eo_impl *m_impl;
};
#endif