10 #ifndef WARPX_PhysicalParticleContainer_H_
11 #define WARPX_PhysicalParticleContainer_H_
52 const std::string&
name);
139 bool skip_deposition=
false,
152 int lev,
int gather_lev,
166 int lev,
int gather_lev,
170 void PushP (
int lev, amrex::Real
dt,
226 amrex::ParticleReal& ux, amrex::ParticleReal& uy, amrex::ParticleReal& uz,
227 amrex::Real t_lab = 0.)
const;
237 amrex::ParticleReal q_tot,
241 amrex::ParticleReal x, amrex::ParticleReal y, amrex::ParticleReal z,
242 amrex::ParticleReal ux, amrex::ParticleReal uy, amrex::ParticleReal uz,
243 amrex::ParticleReal weight,
251 amrex::Real t_lab= 0.)
const;
267 int n_external_attr_real,
268 int n_external_attr_int)
final;
311 amrex::Elixir& exeli, amrex::Elixir& eyeli, amrex::Elixir& ezeli,
312 amrex::Elixir& bxeli, amrex::Elixir& byeli, amrex::Elixir& bzeli,
329 void resample (
int timestep,
bool verbose=
true) final;
355 (
const std::shared_ptr<BreitWheelerEngine>& ptr)
override;
362 (
const std::shared_ptr<QuantumSynchrotronEngine>& ptr)
override;
DtType
Definition: WarpXDtType.H:11
PushType
Definition: WarpXPushType.H:12
Definition: BreitWheelerEngineWrapper.H:294
Filter functor for the Breit Wheeler process.
Definition: QEDPairGeneration.H:38
Filter functor for the QED photon emission process.
Definition: QEDPhotonEmission.H:44
Definition: PhysicalParticleContainer.H:47
void PushP(int lev, amrex::Real dt, const amrex::MultiFab &Ex, const amrex::MultiFab &Ey, const amrex::MultiFab &Ez, const amrex::MultiFab &Bx, const amrex::MultiFab &By, const amrex::MultiFab &Bz) override
Definition: PhysicalParticleContainer.cpp:2545
void set_breit_wheeler_engine_ptr(const std::shared_ptr< BreitWheelerEngine > &ptr) override
Definition: PhysicalParticleContainer.cpp:3401
std::vector< std::string > getUserIntAttribs() const override
Definition: PhysicalParticleContainer.H:378
std::vector< std::string > m_user_int_attribs
Definition: PhysicalParticleContainer.H:440
PhysicalParticleContainer & operator=(PhysicalParticleContainer &&)=default
void ImplicitPushXP(WarpXParIter &pti, amrex::FArrayBox const *exfab, amrex::FArrayBox const *eyfab, amrex::FArrayBox const *ezfab, amrex::FArrayBox const *bxfab, amrex::FArrayBox const *byfab, amrex::FArrayBox const *bzfab, amrex::IntVect ngEB, int, long offset, long np_to_push, int lev, int gather_lev, amrex::Real dt, ScaleFields scaleFields, DtType a_dt_type=DtType::Full)
Definition: PhysicalParticleContainer.cpp:2982
IonizationFilterFunc getIonizationFunc(const WarpXParIter &pti, int lev, amrex::IntVect ngEB, const amrex::FArrayBox &Ex, const amrex::FArrayBox &Ey, const amrex::FArrayBox &Ez, const amrex::FArrayBox &Bx, const amrex::FArrayBox &By, const amrex::FArrayBox &Bz)
Definition: PhysicalParticleContainer.cpp:3314
std::shared_ptr< BreitWheelerEngine > m_shr_p_bw_engine
Definition: PhysicalParticleContainer.H:437
amrex::Vector< amrex::Parser * > getUserIntAttribParser() const override
Definition: PhysicalParticleContainer.H:386
Resampling m_resampler
Definition: PhysicalParticleContainer.H:408
amrex::Vector< std::unique_ptr< amrex::Parser > > m_user_int_attrib_parser
Definition: PhysicalParticleContainer.H:444
PhysicalParticleContainer(PhysicalParticleContainer &&)=default
void Evolve(int lev, const amrex::MultiFab &Ex, const amrex::MultiFab &Ey, const amrex::MultiFab &Ez, const amrex::MultiFab &Bx, const amrex::MultiFab &By, const amrex::MultiFab &Bz, amrex::MultiFab &jx, amrex::MultiFab &jy, amrex::MultiFab &jz, amrex::MultiFab *cjx, amrex::MultiFab *cjy, amrex::MultiFab *cjz, amrex::MultiFab *rho, amrex::MultiFab *crho, const amrex::MultiFab *cEx, const amrex::MultiFab *cEy, const amrex::MultiFab *cEz, const amrex::MultiFab *cBx, const amrex::MultiFab *cBy, const amrex::MultiFab *cBz, amrex::Real t, amrex::Real dt, DtType a_dt_type=DtType::Full, bool skip_deposition=false, PushType push_type=PushType::Explicit) override
Evolve is the central function PhysicalParticleContainer that advances plasma particles for a time dt...
Definition: PhysicalParticleContainer.cpp:2007
void set_quantum_sync_engine_ptr(const std::shared_ptr< QuantumSynchrotronEngine > &ptr) override
Definition: PhysicalParticleContainer.cpp:3408
void SplitParticles(int lev)
Definition: PhysicalParticleContainer.cpp:2345
bool m_do_qed_quantum_sync
Definition: PhysicalParticleContainer.H:428
BreitWheelerEngine * get_breit_wheeler_engine_ptr() const override
Definition: PhysicalParticleContainer.H:365
void ContinuousInjection(const amrex::RealBox &injection_box) override
Definition: PhysicalParticleContainer.cpp:2698
bool has_breit_wheeler() const override
Definition: PhysicalParticleContainer.cpp:3394
PhysicalParticleContainer(PhysicalParticleContainer const &)=delete
bool m_do_qed_breit_wheeler
Definition: PhysicalParticleContainer.H:431
PairGenerationFilterFunc getPairGenerationFilterFunc()
Definition: PhysicalParticleContainer.cpp:3421
bool do_classical_radiation_reaction
Definition: PhysicalParticleContainer.H:421
void MapParticletoBoostedFrame(amrex::ParticleReal &x, amrex::ParticleReal &y, amrex::ParticleReal &z, amrex::ParticleReal &ux, amrex::ParticleReal &uy, amrex::ParticleReal &uz, amrex::Real t_lab=0.) const
Definition: PhysicalParticleContainer.cpp:487
std::shared_ptr< QuantumSynchrotronEngine > m_shr_p_qs_engine
Definition: PhysicalParticleContainer.H:434
void PartitionParticlesInBuffers(long &nfine_current, long &nfine_gather, long np, WarpXParIter &pti, int lev, amrex::iMultiFab const *current_masks, amrex::iMultiFab const *gather_masks)
Definition: Partition.cpp:52
void DefaultInitializeRuntimeAttributes(typename ContainerLike< amrex::PinnedArenaAllocator >::ParticleTileType &pinned_tile, int n_external_attr_real, int n_external_attr_int) final
Default initialize runtime attributes in a tile. This routine does not initialize the first n_externa...
Definition: PhysicalParticleContainer.cpp:838
bool impose_t_lab_from_file
Definition: PhysicalParticleContainer.H:406
void InitData() override
Definition: PhysicalParticleContainer.cpp:481
void ReadHeader(std::istream &is) override
Definition: ParticleIO.cpp:110
virtual void AddParticles(int lev)
Definition: PhysicalParticleContainer.cpp:886
void BackwardCompatibility()
Definition: PhysicalParticleContainer.cpp:465
void resample(int timestep, bool verbose=true) final
This function determines if resampling should be done for the current species, and if so,...
Definition: PhysicalParticleContainer.cpp:3347
amrex::Vector< std::unique_ptr< amrex::Parser > > m_user_real_attrib_parser
Definition: PhysicalParticleContainer.H:446
void AddPlasma(PlasmaInjector const &plasma_injector, int lev, amrex::RealBox part_realbox=amrex::RealBox())
Definition: PhysicalParticleContainer.cpp:955
void CheckAndAddParticle(amrex::ParticleReal x, amrex::ParticleReal y, amrex::ParticleReal z, amrex::ParticleReal ux, amrex::ParticleReal uy, amrex::ParticleReal uz, amrex::ParticleReal weight, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_x, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_y, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_z, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_ux, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_uy, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_uz, amrex::Gpu::HostVector< amrex::ParticleReal > &particle_w, amrex::Real t_lab=0.) const
Definition: PhysicalParticleContainer.cpp:860
QuantumSynchrotronEngine * get_quantum_sync_engine_ptr() const override
Definition: PhysicalParticleContainer.H:369
bool m_save_previous_position
Definition: PhysicalParticleContainer.H:424
void InitIonizationModule() override
Definition: PhysicalParticleContainer.cpp:3229
void PostRestart() final
Definition: PhysicalParticleContainer.H:187
std::vector< std::string > m_user_real_attribs
Definition: PhysicalParticleContainer.H:442
void ContinuousFluxInjection(amrex::Real t, amrex::Real dt) override
Definition: PhysicalParticleContainer.cpp:2710
~PhysicalParticleContainer() override=default
void applyNCIFilter(int lev, const amrex::Box &box, amrex::Elixir &exeli, amrex::Elixir &eyeli, amrex::Elixir &ezeli, amrex::Elixir &bxeli, amrex::Elixir &byeli, amrex::Elixir &bzeli, amrex::FArrayBox &filtered_Ex, amrex::FArrayBox &filtered_Ey, amrex::FArrayBox &filtered_Ez, amrex::FArrayBox &filtered_Bx, amrex::FArrayBox &filtered_By, amrex::FArrayBox &filtered_Bz, const amrex::FArrayBox &Ex, const amrex::FArrayBox &Ey, const amrex::FArrayBox &Ez, const amrex::FArrayBox &Bx, const amrex::FArrayBox &By, const amrex::FArrayBox &Bz, amrex::FArrayBox const *&ex_ptr, amrex::FArrayBox const *&ey_ptr, amrex::FArrayBox const *&ez_ptr, amrex::FArrayBox const *&bx_ptr, amrex::FArrayBox const *&by_ptr, amrex::FArrayBox const *&bz_ptr)
Apply NCI Godfrey filter to all components of E and B before gather.
Definition: PhysicalParticleContainer.cpp:2268
void AddPlasmaFromFile(PlasmaInjector &plasma_injector, amrex::ParticleReal q_tot, amrex::ParticleReal z_shift)
Definition: PhysicalParticleContainer.cpp:711
void WriteHeader(std::ostream &os) const override
Definition: ParticleIO.cpp:117
void AddGaussianBeam(PlasmaInjector const &plasma_injector)
Definition: PhysicalParticleContainer.cpp:534
virtual void PushPX(WarpXParIter &pti, amrex::FArrayBox const *exfab, amrex::FArrayBox const *eyfab, amrex::FArrayBox const *ezfab, amrex::FArrayBox const *bxfab, amrex::FArrayBox const *byfab, amrex::FArrayBox const *bzfab, amrex::IntVect ngEB, int, long offset, long np_to_push, int lev, int gather_lev, amrex::Real dt, ScaleFields scaleFields, DtType a_dt_type=DtType::Full)
Definition: PhysicalParticleContainer.cpp:2729
PlasmaInjector * GetPlasmaInjector(int i) override
Definition: PhysicalParticleContainer.cpp:3338
bool boost_adjust_transverse_positions
Definition: PhysicalParticleContainer.H:401
PhysicalParticleContainer & operator=(PhysicalParticleContainer const &)=delete
std::vector< std::unique_ptr< PlasmaInjector > > plasma_injectors
Definition: PhysicalParticleContainer.H:396
bool has_quantum_sync() const override
Definition: PhysicalParticleContainer.cpp:3389
std::vector< std::string > getUserRealAttribs() const override
Definition: PhysicalParticleContainer.H:382
PhotonEmissionFilterFunc getPhotonEmissionFilterFunc()
Definition: PhysicalParticleContainer.cpp:3414
std::string species_name
Definition: PhysicalParticleContainer.H:395
bool m_rz_random_theta
Definition: PhysicalParticleContainer.H:403
bool do_backward_propagation
Definition: PhysicalParticleContainer.H:402
PhysicalParticleContainer(amrex::AmrCore *amr_core, int ispecies, const std::string &name)
Definition: PhysicalParticleContainer.cpp:247
amrex::Vector< amrex::Parser * > getUserRealAttribParser() const override
Definition: PhysicalParticleContainer.H:390
void AddPlasmaFlux(PlasmaInjector const &plasma_injector, amrex::Real dt)
Definition: PhysicalParticleContainer.cpp:1486
Definition: PlasmaInjector.H:42
Definition: QuantumSyncEngineWrapper.H:273
This is a general class used for resampling that is instantiated as a member of MultiParticleContaine...
Definition: Resampling.H:51
Definition: WarpXParticleContainer.H:53
Definition: WarpXParticleContainer.H:111
ParticleTile< ParticleType, NArrayReal, NArrayInt, Allocator > ParticleTileType
Vector< T * > GetVecOfPtrs(Vector< T > &a)
i
Definition: check_interp_points_and_weights.py:174
name
Definition: run_automated.py:229
float dt
Definition: stencil.py:442
int z_shift
Definition: yt3d_mpi.py:45
Definition: Ionization.H:31
Functor that scales E and B by a factor before pushing the particles. This is used for rigid injectio...
Definition: ScaleFields.H:14