shiftField_eo.h

Go to the documentation of this file.

 
#ifndef SHIFTFIELD_EO_INCLUDED
#define SHIFTFIELD_EO_INCLUDED
 
#include <iostream>
 
#include "field.h"
#include "index_eo.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class ShiftField_eo {
 public:
  static const std::string class_name;
 
 private:
  int m_Nin;
  int m_Nx, m_Ny, m_Nz, m_Nt;
  int m_Nvol;
  int m_Nx2, m_Nvol2;
 
  Index_eo m_index_eo;
 
  Bridge::VerboseLevel m_vl;
 
  std::vector<int> m_yzt_eo;
 
  Field m_wt_x, m_vt_x;  
  Field m_wt_y, m_vt_y;  
  Field m_wt_z, m_vt_z;  
  Field m_wt_t, m_vt_t;  
 
  Field m_we, m_wo;      
  Field m_ve, m_vo;      
  Field m_w1;            
 
 public:
 
  ShiftField_eo(const int Nin) { init(Nin); }
 
  ShiftField_eo() { init(); }
 
 private:
  // non-copyable
  ShiftField_eo(const ShiftField_eo&);
  ShiftField_eo& operator=(const ShiftField_eo&);
 
 public:
  // forward shift (mu: direction, ieo = 0: even<-odd, 1: odd<-even)
  void forward_h(Field&, const Field&, const int mu, const int ieo);
 
  // backward shift (mu: direction, ieo = 0: even<-odd, 1: odd<-even)
  void backward_h(Field&, const Field&, const int mu, const int ieo);
 
  // forward shift with boundary condition
  void forward_h(Field&, const Field&, const int boundary_condition,
                 const int mu, const int ieo);
 
  // backward shift with boundary condition
  void backward_h(Field&, const Field&, const int boundary_condition,
                  const int mu, const int ieo);
 
  void forward(Field&, const Field&, const int mu);
  void backward(Field&, const Field&, const int mu);
 
  void forward(Field&, const Field&,
               const int boundary_condition, const int mu);
  void backward(Field&, const Field&,
                const int boundary_condition, const int mu);
 
 private:
 
  void init(const int Nin = 0);
 
  void up_xh(Field&, const Field&, const int, const int);
  void up_yh(Field&, const Field&, const int, const int);
  void up_zh(Field&, const Field&, const int, const int);
  void up_th(Field&, const Field&, const int, const int);
 
  void dn_xh(Field&, const Field&, const int, const int);
  void dn_yh(Field&, const Field&, const int, const int);
  void dn_zh(Field&, const Field&, const int, const int);
  void dn_th(Field&, const Field&, const int, const int);
};
#endif