Remove reference to theta

dpgrote · Jul 10, 2024 · b4c30c9 · b4c30c9
1 parent 6121572
commit b4c30c9
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 9 deletions.
diff --git a/Source/FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.H b/Source/FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.H
@@ -67,15 +67,8 @@ public:
                       int              a_nl_iter,
                       bool             a_from_jacobian ) override;
 
-    [[nodiscard]] amrex::Real theta () const { return m_theta; }
-
 private:
 
-    /**
-     * \brief Time-biasing parameter for fields used on RHS to advance system
-     */
-    amrex::Real m_theta = 0.5;
-
     /**
      * \brief Solver vectors to be used in the nonlinear solver to solve for the
      * electric field E. The main logic for determining which variables should be

diff --git a/Source/FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.cpp b/Source/FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.cpp
@@ -46,7 +46,6 @@ void StrangImplicitSpectralEM::PrintParameters () const
     amrex::Print() << "------------------------------------------------------------------------" << std::endl;
     amrex::Print() << "----------- STRANG SPLIT IMPLICIT SPECTRAL EM SOLVER PARAMETERS --------" << std::endl;
     amrex::Print() << "------------------------------------------------------------------------" << std::endl;
-    amrex::Print() << "Time-bias parameter theta:  " << m_theta << std::endl;
     amrex::Print() << "max particle iterations:    " << m_max_particle_iterations << std::endl;
     amrex::Print() << "particle tolerance:         " << m_particle_tolerance << std::endl;
     if (m_nlsolver_type==NonlinearSolverType::Picard) {
@@ -108,7 +107,7 @@ void StrangImplicitSpectralEM::ComputeRHS ( WarpXSolverVec& a_Erhs,
                                             bool a_from_jacobian )
 {
     // update WarpX-owned Efield_fp and Bfield_fp using current state of E from
-    // the nonlinear solver at time n+theta
+    // the nonlinear solver at time n+1/2
     UpdateWarpXFields( a_E, a_time, a_dt );
 
     // Self consistently update particle positions and velocities using the