fopr_NonRelativistic.h

Go to the documentation of this file.

 
#ifndef FOPR_NONRELATIVISTIC_INCLUDED
#define FOPR_NONRELATIVISTIC_INCLUDED
 
#include "Fopr/fopr.h"
#include "Field/field_G.h"
 
 
class Fopr_NonRelativistic : public Fopr
{
 public:
  static const std::string class_name;
 
 private:
  double m_MQ;                    
  int m_nstab;                    
  int m_num_correct;              
  double m_u0;                    
 
  Bridge::VerboseLevel m_vl;      
 
  std::vector<double> m_coeff;    
  std::vector<int> m_boundary;    
  std::string m_evolution_type;   
  std::string m_correction_terms; 
 
  std::string m_repr;             
 
  std::string m_mode;             
 
  int m_Nc, m_Nvc, m_Nd, m_Nd2;
  int m_Nx, m_Ny, m_Nz, m_Nt;
  int m_Nvol, m_Ndim;
  int m_Nspc;                         
  int m_Ltime;                        
 
  int m_itime_src;                    
 
  Field_G m_U;                        
 
  std::vector<Field_G> m_Fstr;        
 
  std::vector<double> buf1_x, buf2_x; 
  std::vector<double> buf1_y, buf2_y; 
  std::vector<double> buf1_z, buf2_z; 
 
  std::vector<Field> m_Xt;            
  Field m_Wt1, m_Wt2, m_Wt3;          
  Field m_Yt1, m_Yt2, m_Yt3;          
  Field m_Zt1, m_Zt2;                 
 
 public:
 
 
  /*
  Fopr_NonRelativistic(): Fopr()
  { init(); }
  */
 
  Fopr_NonRelativistic(const Parameters& params) : Fopr()
  { init(params); }
 
  ~Fopr_NonRelativistic() { tidyup(); }
 
  void set_parameters(const Parameters& params);
 
  void set_parameters(const double MQ,
                      const int nstab,
                      const double m_u0,
                      const std::vector<double> coeff,
                      const std::vector<int> bc,
                      const std::string evolution_type,
                      const std::string correction_terms);
 
  void get_parameters(Parameters& params) const;
 
  void set_config(Field *U);
 
  void set_config(unique_ptr<Field_G>& U) { set_config(U.get()); }
 
  void set_mode(const std::string mode);
 
  std::string get_mode() const { return m_mode; }
 
  void mult(Field& v, const Field& w);
 
  void mult_dag(Field& v, const Field& w);
 
  void mult(Field& v, const Field& w, const std::string mode);
 
  void mult_dag(Field& v, const Field& w, const std::string mode);
 
  void mult_gm5(Field& v, const Field& w);
 
  void mult_up(const int mu, Field& v, const Field& w);
 
  void mult_dn(const int mu, Field& v, const Field& w);
 
  int field_nvol() { return m_Nvol; }
 
  int field_nin() { return 2 * m_Nc * m_Nd; }
 
  int field_nex() { return 1; }
 
  double flop_count();
 
  double flop_count(const std::string mode);
 
 private:
 
  void init(const Parameters& params);
 
  void tidyup();
 
  void evolve(Field&, const Field&);
 
  void set_source(Field& Xt, const Field& b);
 
  void evolve_impl(Field&, const Field&, const int itime);
 
  void evolve_typeA(Field&, const Field&, const int itime);
 
  void evolve_typeB(Field&, const Field&, const int itime);
 
  void evolve_H0(Field&, const Field&, const int itime);
 
  void evolve_U4(Field&, const Field&, const int itime);
 
  void evolve_deltaH(Field&, const Field&, const int itime);
 
  void add_deltaH1(Field&, const Field&, const int itime);
 
  void add_deltaH2(Field&, const Field&, const int itime);
 
  void add_deltaH3(Field&, const Field&, const int itime);
 
  void add_deltaH5(Field&, const Field&, const int itime);
 
  void add_deltaH46(Field&, const Field&, const int itime);
 
  void calc_Delta2(Field&, const Field&, const int itime);
 
  void calc_Delta1(Field&, const Field&, const int idir, const int itime);
 
  void calc_Delta4(Field&, const Field&, const int itime);
 
  void rotation(Field& v, const Field& w, const int jd);
 
  void add_R2(Field& v, const Field& w, const int jd, const int itime);
 
  void add_R3(Field& v, const Field& w, const int itime);
 
  void add_R4(Field& v, const Field& w, const int itime);
 
  void add_R5(Field& v, const Field& w, const int itime);
 
  void mult_sigma1(Field& Xt, const Field& Yt);
 
  void mult_sigma2(Field& Xt, const Field& Yt);
 
  void mult_sigma3(Field& Xt, const Field& Yt);
 
  void shift_forward(Field& Xt, const Field& Yt,
                     const int idir, const int itime);
 
  void shift_backward(Field& Xt, const Field& Yt,
                      const int idir, const int itime);
 
  void shift_xup(Field& Xt, const Field& Yt, const int itime);
  void shift_xdn(Field& Xt, const Field& Yt, const int itime);
  void shift_yup(Field& Xt, const Field& Yt, const int itime);
  void shift_ydn(Field& Xt, const Field& Yt, const int itime);
  void shift_zup(Field& Xt, const Field& Yt, const int itime);
  void shift_zdn(Field& Xt, const Field& Yt, const int itime);
 
  void mult_Gn(Field& Xt, const Field& Yt, const int idir, const int itime);
 
  void mult_Gd(Field& Xt, const Field& Yt, const int idir, const int itime);
 
  void mult_F(Field& Xt, const Field& Yt, const int icomp, const int itime);
 
 
#ifdef USE_FACTORY
 private:
  static Fopr *create_object_with_params(const Parameters& params)
  { return new Fopr_NonRelativistic(params); }
 
 public:
  static bool register_factory()
  {
    bool init = Fopr::Factory_params::Register("NonRelativistic",
                                               create_object_with_params);
    return init;
  }
#endif
};
#endif