1 files changed, 29 insertions, 0 deletions
diff --git a/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.H b/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.H
new file mode 100644
index 000000000..a235315ae
--- /dev/null
+++ b/Source/FieldSolver/SpectralSolver/PsatdAlgorithm.H
@@ -0,0 +1,29 @@
+#ifndef WARPX_PSATD_ALGORITHM_H_
+#define WARPX_PSATD_ALGORITHM_H_
+
+#include <SpectralKSpace.H>
+#include <SpectralData.H>
+
+using namespace amrex;
+using namespace Gpu;
+
+/* TODO: Write documentation
+*/
+class PsatdAlgorithm
+{
+    using SpectralCoefficients = FabArray<BaseFab<Real>>;
+
+    public:
+        PsatdAlgorithm(const SpectralKSpace& spectral_kspace,
+                                 const DistributionMapping& dm,
+                                 const int norder_x, const int norder_y,
+                                 const int norder_z, const Real dt);
+        void pushSpectralFields( SpectralData& f ) const;
+
+    private:
+        // Modified finite-order vectors
+        SpectralKVector modified_kx_vec, modified_ky_vec, modified_kz_vec;
+        SpectralCoefficients C_coef, S_ck_coef, X1_coef, X2_coef, X3_coef;
+};
+
+#endif // WARPX_PSATD_ALGORITHM_H_