diff options
| author |  Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
| 2023-03-04 16:10:51 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2023-03-04 16:10:51 -0800 |
| commit | 860a5ebbbcc80f2f9dd5de7f178a331e9712930d (patch) | |
| tree | 3dbf441e233f46815fc429751babe98e88187caa /Source/FieldSolver/FiniteDifferenceSolver | |
| parent | dec1e18adbf3abb8e071df158b22e859ebd73807 (diff) | |
| download | WarpX-860a5ebbbcc80f2f9dd5de7f178a331e9712930d.tar.gz WarpX-860a5ebbbcc80f2f9dd5de7f178a331e9712930d.tar.zst WarpX-860a5ebbbcc80f2f9dd5de7f178a331e9712930d.zip | |
RZ FDTD: fix bug w/ div(E) cleaning (extra c2) (#3726)
Diffstat (limited to 'Source/FieldSolver/FiniteDifferenceSolver')
| -rw-r--r-- | Source/FieldSolver/FiniteDifferenceSolver/EvolveF.cpp | 3 |
1 files changed, 1 insertions, 2 deletions
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/EvolveF.cpp b/Source/FieldSolver/FiniteDifferenceSolver/EvolveF.cpp index c29c0cdb6..27d7fb4e5 100644 --- a/Source/FieldSolver/FiniteDifferenceSolver/EvolveF.cpp +++ b/Source/FieldSolver/FiniteDifferenceSolver/EvolveF.cpp @@ -169,7 +169,6 @@ void FiniteDifferenceSolver::EvolveFCylindrical ( Box const& tf = mfi.tilebox(Ffield->ixType().toIntVect()); Real constexpr inv_epsilon0 = 1./PhysConst::ep0; - Real constexpr c2 = PhysConst::c * PhysConst::c; // Use the right shift in components: // - the first WarpX::ncomps (2*n_rz_azimuthal_modes-1) components correspond to rho old (i.e. rhocomp=0) @@ -195,7 +194,7 @@ void FiniteDifferenceSolver::EvolveFCylindrical ( + T_Algo::DownwardDrr_over_r(Er, r, dr, coefs_r, n_coefs_r, i, j, 0, 2*m-1) + m * Et( i, j, 0, 2*m )/r + T_Algo::DownwardDz(Ez, coefs_z, n_coefs_z, i, j, 0, 2*m-1) ); // Real part - F(i, j, 0, 2*m ) += c2 * dt *( + F(i, j, 0, 2*m ) += dt *( - rho(i, j, 0, rho_shift + 2*m-1) * inv_epsilon0 + T_Algo::DownwardDrr_over_r(Er, r, dr, coefs_r, n_coefs_r, i, j, 0, 2*m-1) - m * Et( i, j, 0, 2*m-1 )/r |