From d266f7019e81c833f9fbfe766ab6de848e8acd2c Mon Sep 17 00:00:00 2001 From: David Grote Date: Tue, 14 Jul 2020 20:46:48 -0700 Subject: Fix bug in RZ PSATD rho terms (#1172) --- .../SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms') diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp index d85b0a537..b094843ac 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp @@ -120,10 +120,10 @@ PsatdAlgorithmRZ::pushSpectralFields(SpectralFieldDataRZ & f) // Update E (see WarpX online documentation: theory section) fields(i,j,k,Ep_m) = C*Ep_old + S_ck*(-c2*I*kr/2._rt*Bz_old + c2*kz*Bp_old - inv_ep0*Jp) - + kr*(X2*rho_new - X3*rho_old); + + 0.5_rt*kr*(X2*rho_new - X3*rho_old); fields(i,j,k,Em_m) = C*Em_old + S_ck*(-c2*I*kr/2._rt*Bz_old - c2*kz*Bm_old - inv_ep0*Jm) - - kr*(X2*rho_new - X3*rho_old); + - 0.5_rt*kr*(X2*rho_new - X3*rho_old); fields(i,j,k,Ez_m) = C*Ez_old + S_ck*(c2*I*kr*Bp_old + c2*I*kr*Bm_old - inv_ep0*Jz) - I*kz*(X2*rho_new - X3*rho_old); -- cgit v1.2.3