aboutsummaryrefslogtreecommitdiff
path: root/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
diff options
context:
space:
mode:
authorGravatar Edoardo Zoni <59625522+EZoni@users.noreply.github.com> 2022-01-09 07:53:02 -0800
committerGravatar GitHub <noreply@github.com> 2022-01-09 07:53:02 -0800
commitbd91b3d5ba42cdf867c132bfc333fac7c887de73 (patch)
tree8dd430ee991860e5c19b7a21321da2b9a04f6d74 /Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
parentb673c598713a8dba4e2477caecabe7a720e13045 (diff)
downloadWarpX-bd91b3d5ba42cdf867c132bfc333fac7c887de73.tar.gz
WarpX-bd91b3d5ba42cdf867c132bfc333fac7c887de73.tar.zst
WarpX-bd91b3d5ba42cdf867c132bfc333fac7c887de73.zip
Multi-J: J Always Linear in Time (#2679)
* Remove WarpX::J_linear_in_time * Fix Bug * Implement Deposition of Rho at Half Time * Implement New Equations for E * Fix New Equations for E * Cleaning * Fix Limits of X2, X3, X7 * Implement New Equations for F * Implement New Equations for F * Revert E Equations to Avoid Roundoff Changes * Remove Quadratic Equations * Remove Deposition of Rho at Half Time Step
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp')
-rw-r--r--Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp8
1 files changed, 5 insertions, 3 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
index 9d8563b6b..86fe8f7bd 100644
--- a/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
+++ b/Source/FieldSolver/SpectralSolver/SpectralFieldData.cpp
@@ -34,7 +34,7 @@ using namespace amrex;
SpectralFieldIndex::SpectralFieldIndex (const bool update_with_rho,
const bool time_averaging,
- const bool J_linear_in_time,
+ const bool do_multi_J,
const bool dive_cleaning,
const bool divb_cleaning,
const bool pml)
@@ -66,9 +66,11 @@ SpectralFieldIndex::SpectralFieldIndex (const bool update_with_rho,
if (divb_cleaning) G = c++;
- if (J_linear_in_time)
+ if (do_multi_J)
{
- Jx_new = c++; Jy_new = c++; Jz_new = c++;
+ Jx_new = c++;
+ Jy_new = c++;
+ Jz_new = c++;
}
}
else // PML