Galilean PSATD: current correction and rho-free formulation (#1151)

* Introduce option to update E with/without rho

* Clean up

* Implement current correction for Galilean PSATD (needs bug fix)

* Include equations in docs

* Fix EOL whitespaces error

* Small clean-up

* Implement Galilean PSATD update without rho

* Clean up

* Fix bug in current correction

* Fix EOL whitespaces

* Clean up

* Fix unused import

* Remove unused variable

* [skip CI] Improve docs

* Clean up style

* Fix EOL whitespaces

* Fix EOL whitespaces

* Clean up style

* Revert analysis script to old status

* [skip CI] Clean up style

* Make equations more human-readable and improve comments

* 2D test with current correction works

* Temporary build fix as in #1197

* 3D test with current correction works

* Rename th and th_star as theta and theta_star

* Fix a couple of wrong comments

* Add vertical spaces to improve readability

* Improve documentation

* Function CurrentCorrection is now pure

* 2D benchmark fields data are now correct

* Add limits of coefficients for nu=-1

* Change default of update_with_rho for Galilean PSATD

1 files changed, 16 insertions, 0 deletions

diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.H b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.H
index 347012d5a..348282dce 100644
--- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.H
+++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.H
@@ -35,6 +35,22 @@ class PMLPsatdAlgorithm : public SpectralBaseAlgorithm
         }
 
         /**
+         * \brief Virtual function for current correction in Fourier space
+         * (<a href="https://doi.org/10.1016/j.jcp.2013.03.010"> Vay et al, 2013</a>).
+         * This function overrides the virtual function \c CurrentCorrection in the
+         * base class \c SpectralBaseAlgorithm and cannot be overridden by further
+         * derived classes.
+         *
+         * \param[in,out] field_data All fields in Fourier space
+         * \param[in,out] current    Array of unique pointers to \c MultiFab storing
+         *                           the three components of the current density
+         * \param[in]     rho        Unique pointer to \c MultiFab storing the charge density
+         */
+        virtual void CurrentCorrection (SpectralFieldData& field_data,
+                                        std::array<std::unique_ptr<amrex::MultiFab>,3>& current,
+                                        const std::unique_ptr<amrex::MultiFab>& rho) override final;
+
+        /**
          * \brief Virtual function for Vay current deposition in Fourier space
          * (<a href="https://doi.org/10.1016/j.jcp.2013.03.010"> Vay et al, 2013</a>).
          * This function overrides the virtual function \c VayDeposition in the