From b4c30c9f263f8d11c0a6df3a42e56aad5ef38550 Mon Sep 17 00:00:00 2001 From: Dave Grote Date: Wed, 10 Jul 2024 08:54:15 -0700 Subject: [PATCH] Remove reference to theta --- .../FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.H | 7 ------- .../ImplicitSolvers/StrangImplicitSpectralEM.cpp | 3 +-- 2 files changed, 1 insertion(+), 9 deletions(-) diff --git a/Source/FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.H b/Source/FieldSolver/ImplicitSolvers/StrangImplicitSpectralEM.H index 3b4f7e107e4..fde0fb44da4 100644 --- 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 index 79830063f27..70204f83f33 100644 --- 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