diff options
| author |  David Grote <grote1@llnl.gov>
| 2020-07-14 20:46:48 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2020-07-14 20:46:48 -0700 |
| commit | d266f7019e81c833f9fbfe766ab6de848e8acd2c (patch) | |
| tree | c70412f7d8fb6dada3d62655c4de218c74ca6993 /Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp | |
| parent | 89b4801d778f58f04162b2c886c6c57ea4cae009 (diff) | |
| download | WarpX-d266f7019e81c833f9fbfe766ab6de848e8acd2c.tar.gz WarpX-d266f7019e81c833f9fbfe766ab6de848e8acd2c.tar.zst WarpX-d266f7019e81c833f9fbfe766ab6de848e8acd2c.zip | |
Fix bug in RZ PSATD rho terms (#1172)
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp')
| -rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
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); |