#ifndef WARPX_PSATD_ALGORITHM_H_ #define WARPX_PSATD_ALGORITHM_H_ #include #include /* TODO: Write documentation */ class PsatdAlgorithm { using SpectralCoefficients = amrex::FabArray< amrex::BaseFab >; public: PsatdAlgorithm( const SpectralKSpace& spectral_kspace, const amrex::DistributionMapping& dm, const int norder_x, const int norder_y, const int norder_z, const amrex::Real dt); PsatdAlgorithm() = default; // Default constructor PsatdAlgorithm& operator=(PsatdAlgorithm&& algorithm) = default; // Default move assignment void pushSpectralFields( SpectralFieldData& f ) const; private: // Modified finite-order vectors KVectorComponent modified_kx_vec, modified_kz_vec; #if (AMREX_SPACEDIM==3) KVectorComponent modified_ky_vec; #endif SpectralCoefficients C_coef, S_ck_coef, X1_coef, X2_coef, X3_coef; }; #endif // WARPX_PSATD_ALGORITHM_H_