diff options
| author |  Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
| 2022-01-09 07:53:02 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2022-01-09 07:53:02 -0800 |
| commit | bd91b3d5ba42cdf867c132bfc333fac7c887de73 (patch) | |
| tree | 8dd430ee991860e5c19b7a21321da2b9a04f6d74 /Source/FieldSolver/WarpXPushFieldsEM.cpp | |
| parent | b673c598713a8dba4e2477caecabe7a720e13045 (diff) | |
| download | WarpX-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/WarpXPushFieldsEM.cpp')
| -rw-r--r-- | Source/FieldSolver/WarpXPushFieldsEM.cpp | 9 |
1 files changed, 3 insertions, 6 deletions
diff --git a/Source/FieldSolver/WarpXPushFieldsEM.cpp b/Source/FieldSolver/WarpXPushFieldsEM.cpp index ea112bbff..a096401bc 100644 --- a/Source/FieldSolver/WarpXPushFieldsEM.cpp +++ b/Source/FieldSolver/WarpXPushFieldsEM.cpp @@ -236,12 +236,9 @@ WarpX::PSATDForwardTransformJ () { const SpectralFieldIndex& Idx = spectral_solver_fp[0]->m_spectral_index; - const int idx_jx = (WarpX::J_linear_in_time) ? static_cast<int>(Idx.Jx_new) - : static_cast<int>(Idx.Jx); - const int idx_jy = (WarpX::J_linear_in_time) ? static_cast<int>(Idx.Jy_new) - : static_cast<int>(Idx.Jy); - const int idx_jz = (WarpX::J_linear_in_time) ? static_cast<int>(Idx.Jz_new) - : static_cast<int>(Idx.Jz); + const int idx_jx = (WarpX::do_multi_J) ? static_cast<int>(Idx.Jx_new) : static_cast<int>(Idx.Jx); + const int idx_jy = (WarpX::do_multi_J) ? static_cast<int>(Idx.Jy_new) : static_cast<int>(Idx.Jy); + const int idx_jz = (WarpX::do_multi_J) ? static_cast<int>(Idx.Jz_new) : static_cast<int>(Idx.Jz); for (int lev = 0; lev <= finest_level; ++lev) { |