fopr_Wilson.h

Go to the documentation of this file.

#ifndef FOPR_WILSON_INCLUDED
#define FOPR_WILSON_INCLUDED

#include <valarray>
#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::valarray<int> bc);

  void set_config(Field *U);

  void set_mode(std::string mode);

  std::string get_mode() const;

  inline const Field mult(const Field& f)
  {
    Field v(f.nin(), f.nvol(), f.nex());

    mult(v, f);
    return v;
  }

  inline const Field mult_dag(const Field& f)
  {
    Field v(f.nin(), f.nvol(), f.nex());

    mult_dag(v, f);
    return v;
  }

  inline const Field mult_gm5(const Field& f)
  {
    Field v(f.nin(), f.nvol(), f.nex());

    mult_gm5(v, f);
    return v;
  }

  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 gm5p(Field_F& v, const int mu, const Field_F& w);

  const Field_F mult_gm5p(int mu, const Field_F& w);
  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);

  //  void fprop_normalize(Field& v);
  //  void fopr_normalize(Field& v);

  //  const double get_fprop_normfactor();
  //  const double get_fopr_normfactor();

  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