fft_xyz_1dim.h

Go to the documentation of this file.

#ifndef FFT_XYZ_1DIM_INCLUDED
#define FFT_XYZ_1DIM_INCLUDED

#ifdef USE_FFTWLIB

#ifdef USE_MPI
#include "fftw3-mpi.h"
#ifdef USE_BGNET
#include "communicator_bgnet.h"
#else
#include "communicator_mpi.h"
#endif
#else
#include "fftw3.h"
#endif

#ifdef USE_OPENMP
#include "threadManager_OpenMP.h"
#endif

#include "field.h"
#include "index_lex.h"

#include "bridgeIO.h"
using Bridge::vout;


class FFT_xyz_1dim
{
 public:
  static const std::string class_name;

 protected:
  Bridge::VerboseLevel m_vl;

 public:
  FFT_xyz_1dim()
    : m_vl(CommonParameters::Vlevel()) {}

  ~FFT_xyz_1dim() {}

 private:
  // non-copyable
  FFT_xyz_1dim(const FFT_xyz_1dim&);
  FFT_xyz_1dim& operator=(const FFT_xyz_1dim&);

 public:
  void FFT(Field& field, const bool is_forward);  // field is overwritten
  void FFT(Field& field_out, const Field& field_in, const bool is_forward);

 private:
  Index_lex m_index;

  fftw_complex *m_in;
  fftw_complex *m_out;

  fftw_plan m_plan;

#ifdef USE_MPI
  ptrdiff_t m_Nsize_in_p, m_start_in_p;
  ptrdiff_t m_Nsize_out_p, m_start_out_p;
#endif
};
//- #endif of #ifdef USE_FFTWLIB
#endif
#endif