fopr_Wilson_eo_impl.h

Go to the documentation of this file.

#ifndef FOPR_WILSON_EO_ORG_INCLUDED
#define FOPR_WILSON_EO_ORG_INCLUDED

#include "Fopr/fopr_eo.h"
#include "Field/shiftField_eo.h"
#include "Tools/gammaMatrixSet.h"

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


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

   private:
    int              m_Nvol, m_Nvol2, m_Ndim;
    int              m_Nc, m_Nd;
    double           m_kappa;
    std::vector<int> m_boundary;
    Index_eo         m_index;
    Field_G          *m_Ueo;

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

    std::vector<GammaMatrix> m_GM;

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

   public:

    Fopr_Wilson_eo()
      : m_Nvol(CommonParameters::Nvol()),
        m_Nvol2(CommonParameters::Nvol() / 2),
        m_Ndim(CommonParameters::Ndim()),
        m_Nc(CommonParameters::Nc()),
        m_Nd(CommonParameters::Nd()),
        trf(m_Nvol2, 1), trf2(m_Nvol2, 1)
    {
      init("Dirac");
    }

    Fopr_Wilson_eo(std::string repr)
      : m_Nvol(CommonParameters::Nvol()),
        m_Nvol2(CommonParameters::Nvol() / 2),
        m_Ndim(CommonParameters::Ndim()),
        m_Nc(CommonParameters::Nc()),
        m_Nd(CommonParameters::Nd()),
        trf(m_Nvol2, 1), trf2(m_Nvol2, 1)
    {
      init(repr);
    }

    ~Fopr_Wilson_eo()
    {
      delete m_Ueo;
    }

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

    void set_config(Field *U);

    void set_config(unique_ptr<Field_G>& U)
    { return set_config(U.get()); }

    void set_mode(std::string mode);
    std::string get_mode() const;

    // 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&);

    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 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);

    // ieo=0: even <-- odd
    // ieo=1: odd  <-- even

    void Meo(Field&, const Field&, const int ieo);
    void Mdageo(Field&, const Field&, const int ieo);
    void MeoMoe(Field& v, const Field& f);
    void Meo_gm5(Field&, const Field&, const int ieo);

    void mult_gm5(Field&, const Field&);

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

    int field_nvol() { return CommonParameters::Nvol() / 2; }
    int field_nin() { return 2 * CommonParameters::Nc() * CommonParameters::Nd(); }
    int field_nex() { return 1; }

    double flop_count();

   private:
    void init(const std::string);

    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&);

    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);

    /*
    void mult_xp(Field_F&, const Field_F&, const int ieo);
    void mult_xm(Field_F&, const Field_F&, const int ieo);
    void mult_yp(Field_F&, const Field_F&, const int ieo);
    void mult_ym(Field_F&, const Field_F&, const int ieo);
    void mult_zp(Field_F&, const Field_F&, const int ieo);
    void mult_zm(Field_F&, const Field_F&, const int ieo);
    void mult_tp(Field_F&, const Field_F&, const int ieo);
    void mult_tm(Field_F&, const Field_F&, const int ieo);
    */
  };
}
#endif