From f4978e1001494e2b148128380fb37f3c2450209f Mon Sep 17 00:00:00 2001 From: Andrew Myers Date: Mon, 22 Jun 2020 09:28:16 -0700 Subject: 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 * 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 --- Source/Particles/Gather/GetExternalFields.cpp | 46 +++++++++++++++++++++++++++ 1 file changed, 46 insertions(+) create mode 100644 Source/Particles/Gather/GetExternalFields.cpp (limited to 'Source/Particles/Gather/GetExternalFields.cpp') 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(); + } +} -- cgit v1.2.3