Remove persistent E+B (#1050)

* add functor for doing the tmp particles copy for the back-transformed diagnosti

* merge the particle push options into one kernel

* EOL

* fix assertion

* add a FieldGatherandPushPX method to PhysicalParticleContainer

* handle offset in copyAttribs

* allow this functor to be constructed even it we aren't doing the back transformed diagnostics

* EOL

* update the overloads of PushPX for the Photon and RigidInjected ParticleContainers

* function for dispatching the right field gather

* init this val to 0.0

* fix some typos

* handle scaling the fields for rigid injection

* EOL

* don't need to get pointers to E and B arrays in PushPX any more.

* actually I can't remove these yet

* EOL

* variable order bug

* move the QED stuff to the proper place

* EOL

* make sure we don't build these functors unless the runtime options are toggled

* EOL

* perform the field gather prior to the photon particle push

* remove E and B components and FieldGather methods. Reimplement PushP for rigid injected and physical particles

* update ionization to do field gather inline

* remove E and B from the particle diagnostics

* don't write E or B in these tests for particles

* add missing files

* remove EB from the Regtest ini file too

* no need to do this twice

* important typo

* also do the gather inline for the QED processes that need to

* move these sources inside ifdef for QED

* fix bug in RZ

* remove some fields from the Python tests.

* remove all particle E and B comps from json benchmarks

* don't assert that Ey is the langmuir output

* remove uy from this output

* update test

* restore the mesh fields I turned off by mistake

* turn off field IO for a few python tests I missed

* fix typo

* reset Langmuir_multi benchmark

* update Langmuir_multi_nodal benchmark

* update single precision langmuir bench

* update psatd single precision languir one too

* also do ionization_lab

* finally, ionizaiton_boost

* update benchmarks_json/Langmuir_multi_psatd.json

* update benchmarks_json/Langmuir_multi_psatd_current_correction.json

* update benchmarks_json/Langmuir_multi_psatd_momentum_conserving.json

* update benchmarks_json/Langmuir_multi_psatd_nodal.json

* remove the particle E and B from the choices in the docs

* fix offset bug

* also add the Gather subdirectory

* Update Source/WarpX.H

Co-authored-by: MaxThevenet <mthevenet@lbl.gov>

* add docstring for LowerCornerWithGalilean

* add new source files to CMakeLists.txt

* also need to update the GPU regression tests

* update the name of the output file for this python test

* remove field gather call from FieldDiagnostics

* fix typo in docstring

* init fields to 0

* add docstring to the CopyParticleAttribs constructor

* some explicit amrex::namepace

Co-authored-by: MaxThevenet <mthevenet@lbl.gov>

1 files changed, 46 insertions, 0 deletions

diff --git a/Source/Particles/Gather/GetExternalFields.cpp b/Source/Particles/Gather/GetExternalFields.cpp
new file mode 100644
index 000000000..6e4eb6d1b
--- /dev/null
+++ b/Source/Particles/Gather/GetExternalFields.cpp
@@ -0,0 +1,46 @@
+#include "WarpX.H"
+#include "Particles/Gather/GetExternalFields.H"
+
+GetExternalEField::GetExternalEField (const WarpXParIter& a_pti, int a_offset) noexcept
+{
+    auto& warpx = WarpX::GetInstance();
+    auto& mypc = warpx.GetPartContainer();
+    if (mypc.m_E_ext_particle_s=="constant" || mypc.m_E_ext_particle_s=="default")
+    {
+        m_type = Constant;
+        m_field_value[0] = mypc.m_E_external_particle[0];
+        m_field_value[1] = mypc.m_E_external_particle[1];
+        m_field_value[2] = mypc.m_E_external_particle[2];
+    }
+    else if (mypc.m_E_ext_particle_s=="parse_e_ext_particle_function")
+    {
+        m_type = Parser;
+        m_time = warpx.gett_new(a_pti.GetLevel());
+        m_get_position = GetParticlePosition(a_pti, a_offset);
+        m_xfield_partparser = mypc.m_Ex_particle_parser.get();
+        m_yfield_partparser = mypc.m_Ey_particle_parser.get();
+        m_zfield_partparser = mypc.m_Ez_particle_parser.get();
+    }
+}
+
+GetExternalBField::GetExternalBField (const WarpXParIter& a_pti, int a_offset) noexcept
+{
+    auto& warpx = WarpX::GetInstance();
+    auto& mypc = warpx.GetPartContainer();
+    if (mypc.m_B_ext_particle_s=="constant" || mypc.m_B_ext_particle_s=="default")
+    {
+        m_type = Constant;
+        m_field_value[0] = mypc.m_B_external_particle[0];
+        m_field_value[1] = mypc.m_B_external_particle[1];
+        m_field_value[2] = mypc.m_B_external_particle[2];
+    }
+    else if (mypc.m_B_ext_particle_s=="parse_e_ext_particle_function")
+    {
+        m_type = Parser;
+        m_time = warpx.gett_new(a_pti.GetLevel());
+        m_get_position = GetParticlePosition(a_pti, a_offset);
+        m_xfield_partparser = mypc.m_Bx_particle_parser.get();
+        m_yfield_partparser = mypc.m_By_particle_parser.get();
+        m_zfield_partparser = mypc.m_Bz_particle_parser.get();
+    }
+}