Multi-J: J Always Linear in Time (#2679)

* Remove WarpX::J_linear_in_time

* Fix Bug

* Implement Deposition of Rho at Half Time

* Implement New Equations for E

* Fix New Equations for E

* Cleaning

* Fix Limits of X2, X3, X7

* Implement New Equations for F

* Implement New Equations for F

* Revert E Equations to Avoid Roundoff Changes

* Remove Quadratic Equations

* Remove Deposition of Rho at Half Time Step

1 files changed, 3 insertions, 3 deletions

diff --git a/Source/FieldSolver/SpectralSolver/SpectralSolverRZ.cpp b/Source/FieldSolver/SpectralSolver/SpectralSolverRZ.cpp
index b47dfa4ad..e187b9d27 100644
--- a/Source/FieldSolver/SpectralSolver/SpectralSolverRZ.cpp
+++ b/Source/FieldSolver/SpectralSolver/SpectralSolverRZ.cpp
@@ -34,7 +34,7 @@ SpectralSolverRZ::SpectralSolverRZ (const int lev,
                                     amrex::RealVect const dx, amrex::Real const dt,
                                     bool const update_with_rho,
                                     const bool fft_do_time_averaging,
-                                    const bool J_linear_in_time,
+                                    const bool do_multi_J,
                                     const bool dive_cleaning,
                                     const bool divb_cleaning)
     : k_space(realspace_ba, dm, dx)
@@ -47,7 +47,7 @@ SpectralSolverRZ::SpectralSolverRZ (const int lev,
 
     const bool pml = false;
     m_spectral_index = SpectralFieldIndex(update_with_rho, fft_do_time_averaging,
-                                          J_linear_in_time, dive_cleaning, divb_cleaning, pml);
+                                          do_multi_J, dive_cleaning, divb_cleaning, pml);
 
     // - Select the algorithm depending on the input parameters
     //   Initialize the corresponding coefficients over k space
@@ -56,7 +56,7 @@ SpectralSolverRZ::SpectralSolverRZ (const int lev,
          // v_galilean is 0: use standard PSATD algorithm
         algorithm = std::make_unique<PsatdAlgorithmRZ>(
             k_space, dm, m_spectral_index, n_rz_azimuthal_modes, norder_z, nodal, dt,
-            update_with_rho, fft_do_time_averaging, J_linear_in_time, dive_cleaning, divb_cleaning);
+            update_with_rho, fft_do_time_averaging, do_multi_J, dive_cleaning, divb_cleaning);
     } else {
         // Otherwise: use the Galilean algorithm
         algorithm = std::make_unique<GalileanPsatdAlgorithmRZ>(