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 "bridgeIO.h"
using Bridge::vout;


class ShiftField_eo {
 private:
  int                  Nx, Ny, Nz, Nt, Nx2;
  Index_eo             idx;
  Bridge::VerboseLevel m_vl;

 public:
  ShiftField_eo() :
    Nx(CommonParameters::Nx()),
    Ny(CommonParameters::Ny()),
    Nz(CommonParameters::Nz()),
    Nt(CommonParameters::Nt()),
    Nx2(CommonParameters::Nx() / 2),
    m_vl(CommonParameters::Vlevel())
  {
  }

  void forward_h(Field&, const Field&, const int mu, const int ieo);
  void backward_h(Field&, const Field&, const int mu, const int ieo);

  //                                    direction, ieo (0: e<-o, 1: o<-e)

  void forward_h(Field&, const Field&, const int bc, const int mu,
                 const int ieo);
  void backward_h(Field&, const Field&, const int bc, const int mu,
                  const int ieo);

  //                boundary condition, direction, ieo (0: e<-o, 1: o<-e)

  void forward(Field&, const Field&, const int mu);
  void backward(Field&, const Field&, const int mu);

  void forward(Field&, const Field&, const int bc, const int mu);
  void backward(Field&, const Field&, const int bc, const int mu);



 private:
  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