SpectralFieldDataRZ.H

Go to the documentation of this file.

/* Copyright 2019-2020 David Grote
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef WARPX_SPECTRAL_FIELD_DATA_RZ_H_
#define WARPX_SPECTRAL_FIELD_DATA_RZ_H_
 
#include "SpectralBinomialFilter.H"
#include "SpectralFieldData.H"
#include "SpectralHankelTransform/SpectralHankelTransformer.H"
#include "SpectralKSpaceRZ.H"
 
#include <ablastr/math/fft/AnyFFT.H>
 
#include <AMReX_MultiFab.H>
 
/* \brief Class that stores the fields in spectral space, and performs the
 *  Fourier transforms between real space and spectral space
 */
class SpectralFieldDataRZ
{
 
    public:
 
        // Define the FFTplans type, which holds one fft plan per box
        // (plans are only initialized for the boxes that are owned by
        // the local MPI rank)
        using FFTplans = amrex::LayoutData<ablastr::math::anyfft::VendorFFTPlan>;
 
        // Similarly, define the Hankel transformers and filter for each box.
        using MultiSpectralHankelTransformer = amrex::LayoutData<SpectralHankelTransformer>;
 
        using BinomialFilter = amrex::LayoutData<SpectralBinomialFilter>;
 
        SpectralFieldDataRZ (int lev,
                             const amrex::BoxArray& realspace_ba,
                             const SpectralKSpaceRZ& k_space,
                             const amrex::DistributionMapping& dm,
                             int n_field_required,
                             int n_modes);
        SpectralFieldDataRZ () = default; // Default constructor
        ~SpectralFieldDataRZ ();

        SpectralFieldDataRZ ( SpectralFieldDataRZ const &) = delete;
        SpectralFieldDataRZ& operator= ( SpectralFieldDataRZ const & ) = delete;
        SpectralFieldDataRZ ( SpectralFieldDataRZ&& ) = default;
        SpectralFieldDataRZ& operator=(SpectralFieldDataRZ&& field_data) = default;
 
 
        void ForwardTransform (int lev, const amrex::MultiFab& mf, int field_index,
                               int i_comp=0);
        void ForwardTransform (int lev, const amrex::MultiFab& mf_r, int field_index_r,
                               const amrex::MultiFab& mf_t, int field_index_t);
        void BackwardTransform (int lev, amrex::MultiFab& mf, int field_index,
                                int i_comp=0);
        void BackwardTransform (int lev, amrex::MultiFab& mf_r, int field_index_r,
                                amrex::MultiFab& mf_t, int field_index_t);
 
        void FABZForwardTransform (amrex::MFIter const & mfi, amrex::Box const & realspace_bx,
                                   amrex::MultiFab const & tempHTransformedSplit,
                                   int field_index, bool is_nodal_z);
        void FABZBackwardTransform (amrex::MFIter const & mfi, amrex::Box const & realspace_bx,
                                    int field_index,
                                    amrex::MultiFab & tempHTransformedSplit,
                                    bool is_nodal_z);

        void CopySpectralDataComp (int src_comp, int dest_comp);

        void ZeroOutDataComp(int icomp);

        void ScaleDataComp(int icomp, amrex::Real scale_factor);
 
        void InitFilter (amrex::IntVect const & filter_npass_each_dir, bool compensation,
                         SpectralKSpaceRZ const & k_space);
 
        void ApplyFilter (int lev, int field_index);
        void ApplyFilter (int lev, int field_index1,
                          int field_index2, int field_index3);
 
        // Returns an array that holds the kr for all of the modes
        HankelTransform::RealVector const & getKrArray (amrex::MFIter const & mfi) const {
            return multi_spectral_hankel_transformer[mfi].getKrArray();
        }

        SpectralField fields;

        int n_rz_azimuthal_modes;
        int m_ncomps;
 
    private:
 
        SpectralFieldIndex m_spectral_index;
        int m_n_fields;
 
        // tempHTransformed and tmpSpectralField store fields
        // right before/after the z Fourier transform
        SpectralField tempHTransformed; // contains Complexes
        SpectralField tmpSpectralField; // contains Complexes
        FFTplans forward_plan, backward_plan;
        // Correcting "shift" factors when performing FFT from/to
        // a cell-centered grid in real space, instead of a nodal grid
        SpectralShiftFactor zshift_FFTfromCell, zshift_FFTtoCell;
        MultiSpectralHankelTransformer multi_spectral_hankel_transformer;
        BinomialFilter binomialfilter;
 
};
 
#endif // WARPX_SPECTRAL_FIELD_DATA_RZ_H_