fopr_Wilson.h

Go to the documentation of this file.

#ifndef FOPR_WILSON_INCLUDED
#define FOPR_WILSON_INCLUDED

#include <string>
#include <cassert>

#include "fopr.h"
#include "field_F.h"

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


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

class Fopr_Wilson : public Fopr
{
 public:
  static const std::string class_name;

 public:
  // interface definitions

  Fopr_Wilson();
  Fopr_Wilson(std::string repr);

  ~Fopr_Wilson();

  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)
  {
    set_config(U.get());
  }

  void set_mode(std::string mode);

  std::string get_mode() const;

  void mult(Field& v, const Field& f);
  void mult_dag(Field& v, const Field& f);
  void mult_gm5(Field& w, const Field& v);
  void proj_chiral(Field& w, const int ex1,
                   const Field& v, const int ex2, const int ipm);

  void D(Field& w, const Field& v);
  void Ddag(Field& w, const Field& v);
  void DdagD(Field& w, const Field& v);
  void DDdag(Field& w, const Field& v);
  void H(Field& w, const Field& v);

  void D_ex(Field& w, const int ex1, const Field& v, const int ex2);

  void mult_gm5p(int mu, Field_F& v, const Field_F& w);

  void mult_up(int mu, Field& w, const Field& v);
  void mult_dn(int mu, Field& w, const Field& v);


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

  double flop_count();

 private:
  // pimple prescription for Fopr_Wilson class implementation.
  class Fopr_Wilson_impl;

  Fopr_Wilson_impl *m_impl;
};
#endif