AnyFFT.H

Go to the documentation of this file.

/* Copyright 2019-2020
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */

#ifndef ANYFFT_H_
#define ANYFFT_H_

#ifdef AMREX_USE_GPU
#  include <cufft.h>
#else
#  include <fftw3.h>
#endif

#include <AMReX_LayoutData.H>

namespace AnyFFT
{
    // First, define library-dependent types (complex, FFT plan)

#ifdef AMREX_USE_GPU
#  ifdef AMREX_USE_FLOAT
    using Complex = cuComplex;
#  else
    using Complex = cuDoubleComplex;
#  endif
#else
#  ifdef AMREX_USE_FLOAT
    using Complex = fftwf_complex;
#  else
    using Complex = fftw_complex;
#  endif
#endif

#ifdef AMREX_USE_GPU
    using VendorFFTPlan = cufftHandle;
#else
#  ifdef AMREX_USE_FLOAT
    using VendorFFTPlan = fftwf_plan;
#  else
    using VendorFFTPlan = fftw_plan;
#  endif
#endif

    // Second, define library-independent API

    enum struct direction {R2C, C2R};

    struct FFTplan
    {
        amrex::Real* m_real_array; 
        Complex* m_complex_array; 
        VendorFFTPlan m_plan; 
        direction m_dir;  
        int m_dim; 
    };

    using FFTplans = amrex::LayoutData<FFTplan>;

    FFTplan CreatePlan(const amrex::IntVect& real_size, amrex::Real * const real_array,
                       Complex * const complex_array, const direction dir, const int dim);

    void DestroyPlan(FFTplan& fft_plan);

    void Execute(FFTplan& fft_plan);
}

#endif // ANYFFT_H_