fprop_alt_Standard_eo_Richardson.h

Go to the documentation of this file.

 
#ifndef FPROP_ALT_STANDARD_EO_RICHARDSON_INCLUDED
#define FPROP_ALT_STANDARD_EO_RICHARDSON_INCLUDED
 
#include "lib_alt/Measurements/Fermion/fprop_alt.h"
 
#include "lib/Tools/timer.h"
 
#include "lib/Fopr/afopr.h"
 
#include "lib_alt/Solver/asolver.h"
#include "lib_alt/Solver/aprecond.h"
//#include "Field/afield.h"
 
 
template<typename AFIELD_d, typename AFIELD_f>
class Fprop_alt_Standard_eo_Richardson : public Fprop_alt<AFIELD_d>
{
 public:
  static const std::string class_name;
  typedef typename AFIELD_d::real_t real_t;
  using Fprop_alt<AFIELD_d>::m_vl;
  using Fprop_alt<AFIELD_d>::m_mode;
 
 private:
  AFopr<AFIELD_d> *m_kernel_d;
  AFopr<AFIELD_d> *m_fopr_d;
  AFopr<AFIELD_f> *m_kernel_f;
  AFopr<AFIELD_f> *m_fopr_f;
  ASolver<AFIELD_f> *m_solver_prec;
  APrecond<AFIELD_d> *m_precond;
  ASolver<AFIELD_d> *m_solver;
 
  AFIELD_d m_axq, m_abq;
  AFIELD_d m_be, m_bo, m_xe, m_xo, m_y1, m_y2;
  AFIELD_f m_y1f, m_y2f;
 
  Timer m_timer;
  double m_flop_count;
  double m_elapsed_time;
 
  Director_Smear *m_dr_smear;
 
  bool m_use_DdagD;  // must be true for domainwall
 
 public:
  Fprop_alt_Standard_eo_Richardson(const Parameters& params_fopr,
                                   const Parameters& params_solver,
                                   Director_Smear *dr_smear)
    : Fprop_alt<AFIELD_d>()
  { init(params_fopr, params_solver, dr_smear); }
 
  Fprop_alt_Standard_eo_Richardson(const Parameters& params_fopr,
                                   const Parameters& params_solver)
    : Fprop_alt<AFIELD_d>()
  { init(params_fopr, params_solver); }
 
  ~Fprop_alt_Standard_eo_Richardson()
  { tidyup(); }
 
  void set_config(Field *);
 
  void invert(Field&, const Field&, int&, double&);
 
  void invert_D(Field&, const Field&, int&, double&);
 
  void invert_DdagD(Field&, const Field&, int&, double&);
 
  void invert(AFIELD_d&, const AFIELD_d&, int&, double&);
 
  void invert_D(AFIELD_d&, const AFIELD_d&, int&, double&);
 
  void invert_DdagD(AFIELD_d&, const AFIELD_d&, int&, double&);
 
  double flop_count();
 
  void reset_performance();
 
  void report_performance();
 
  void mult_performance(const std::string mode, const int Nrepeat);
 
 private:
 
  void init(const Parameters& params_fopr,
            const Parameters& params_solver);
 
  void init(const Parameters& params_fopr,
            const Parameters& params_solver,
            Director_Smear *dr_smear);
 
  void tidyup();
 
  void invert_De(AFIELD_d&, AFIELD_d&, AFIELD_d&, AFIELD_d&,
                 int&, double&);
};
#endif