fopr_Clover_eo.h

Go to the documentation of this file.

 
#ifndef FOPR_CLOVER_EO_INCLUDED
#define FOPR_CLOVER_EO_INCLUDED
 
#include "Fopr/fopr_Wilson_eo.h"
#include "Fopr/fopr_CloverTerm_eo.h"
#include "Field/field_F.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
 
class Fopr_Clover_eo : public Fopr_eo
{
 public:
  static const std::string class_name;
 
 private:
  // input parameters
  double m_kappa;               
  double m_cSW;                 
  std::vector<int> m_boundary;  
  std::string m_repr;           
  Bridge::VerboseLevel m_vl;    
 
  std::string m_mode;
 
  // internal data members
  int m_Nc, m_Nd, m_NinF, m_Ndim;
  int m_Nvol, m_Nvol2;
 
  Fopr_Wilson_eo *m_fopr_w;
  Fopr_CloverTerm_eo *m_fopr_csw;
 
  Index_eo m_idx;
  Field_G *m_U;
 
  Field m_w1, m_w2;       
  Field m_v1;             
 
  bool m_is_initial_step; 
 
 public:
  DEPRECATED
  Fopr_Clover_eo(const std::string repr) { init(repr); }
 
  Fopr_Clover_eo(const Parameters& params) { init(params); }
 
  ~Fopr_Clover_eo() { tidyup(); }
 
  void set_parameters(const Parameters& params);
 
  void set_parameters(const double kappa,
                      const double cSW,
                      const std::vector<int> bc);
 
  void get_parameters(Parameters& params) const;
 
  void set_config(Field *U);
 
  void set_mode(const std::string mode);
 
  std::string get_mode() const { return m_mode; }
 
  void mult(Field& v, const Field& f);
 
  void mult_dag(Field& v, const Field& f);
 
  //- method for even odd fermion operator
  void preProp(Field& Be, Field& bo, const Field& b);
 
  void postProp(Field& x, const Field& xe, const Field& bo);
 
  //  const Field_F mult_csw_inv(const Field_F&, const int ieo);
 
  void mult_gm5(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);
 
  // Meo: 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 mult_isigma(Field_F& w, const Field_F& f,
                   const int mu, const int nu);
 
  int field_nin()  { return 2 * m_Nc * m_Nd; }
  int field_nvol() { return m_Nvol2; }
  int field_nex()  { return 1; }
 
  double flop_count();
 
 private:
  void init(const std::string repr);
 
  void init(const Parameters& params);
 
  void setup();
 
  void tidyup();
 
#ifdef USE_FACTORY
 private:
  static Fopr *create_object_with_repr(const std::string& repr)
  {
    return new Fopr_Clover_eo(repr);
  }
 
  static Fopr *create_object_with_params(const Parameters& params)
  {
    return new Fopr_Clover_eo(params);
  }
 
 public:
  static bool register_factory()
  {
    bool init = true;
    init &= Fopr::Factory_string::Register("Clover_eo",
                                           create_object_with_repr);
    init &= Fopr::Factory_params::Register("Clover_eo",
                                           create_object_with_params);
    return init;
  }
#endif
};
#endif