8 #ifndef WARPX_NUCLEAR_FUSION_FUNC_H_
9 #define WARPX_NUCLEAR_FUSION_FUNC_H_
60 const bool isSameSpecies):
68 #ifdef AMREX_SINGLE_PRECISION_PARTICLES
78 pp_collision_name,
"fusion_probability_target_value",
134 amrex::ParticleReal
const , amrex::ParticleReal
const ,
135 amrex::ParticleReal
const m1, amrex::ParticleReal
const m2,
136 amrex::Real
const dt, amrex::Real
const dV,
index_type coll_idx,
139 amrex::ParticleReal* AMREX_RESTRICT p_pair_reaction_weight,
158 index_type pair_index = cell_start_pair + coll_idx;
170 const auto multiplier_ratio =
static_cast<int>(
173 #if (defined WARPX_DIM_RZ)
181 for (
index_type k = coll_idx; k < max_N; k += min_N)
185 const index_type c1k = (k%NI1 < max_N%NI1) ? c1 + 1: c1;
186 const index_type c2k = (k%NI2 < max_N%NI2) ? c2 + 1: c2;
188 #if (defined WARPX_DIM_RZ)
197 amrex::ParticleReal
const theta = theta2[I2[i2]]-theta1[I1[i1]];
198 amrex::ParticleReal
const u1xbuf = u1x[I1[i1]];
199 u1x[I1[i1]] = u1xbuf*std::cos(theta) - u1y[I1[i1]]*std::sin(theta);
200 u1y[I1[i1]] = u1xbuf*std::sin(theta) + u1y[I1[i1]]*std::cos(theta);
204 u1x[ I1[i1] ], u1y[ I1[i1] ], u1z[ I1[i1] ],
205 u2x[ I2[i2] ], u2y[ I2[i2] ], u2z[ I2[i2] ],
206 m1, m2, w1[ I1[i1] ]/c1k, w2[ I2[i2] ]/c2k,
207 dt, dV,
static_cast<int>(pair_index), p_mask, p_pair_reaction_weight,
213 #if (defined WARPX_DIM_RZ)
214 amrex::ParticleReal
const u1xbuf_new = u1x[I1[i1]];
215 u1x[I1[i1]] = u1xbuf_new*std::cos(-theta) - u1y[I1[i1]]*std::sin(-theta);
216 u1y[I1[i1]] = u1xbuf_new*std::sin(-theta) + u1y[I1[i1]]*std::cos(-theta);
219 p_pair_indices_1[pair_index] = I1[i1];
220 p_pair_indices_2[pair_index] = I2[i2];
#define AMREX_GPU_HOST_DEVICE
NuclearFusionType
Definition: BinaryCollisionUtils.H:26
AMREX_GPU_HOST_DEVICE AMREX_INLINE void SingleNuclearFusionEvent(const amrex::ParticleReal &u1x, const amrex::ParticleReal &u1y, const amrex::ParticleReal &u1z, const amrex::ParticleReal &u2x, const amrex::ParticleReal &u2y, const amrex::ParticleReal &u2z, const amrex::ParticleReal &m1, const amrex::ParticleReal &m2, amrex::ParticleReal w1, amrex::ParticleReal w2, const amrex::Real &dt, const amrex::ParticleReal &dV, const int &pair_index, index_type *AMREX_RESTRICT p_mask, amrex::ParticleReal *AMREX_RESTRICT p_pair_reaction_weight, const amrex::ParticleReal &fusion_multiplier, const int &multiplier_ratio, const amrex::ParticleReal &probability_threshold, const amrex::ParticleReal &probability_target_value, const NuclearFusionType &fusion_type, const amrex::RandomEngine &engine)
This function computes whether the collision between two particles result in a nuclear fusion event,...
Definition: SingleNuclearFusionEvent.H:54
#define WARPX_ABORT_WITH_MESSAGE(MSG)
Definition: TextMsg.H:15
Definition: MultiParticleContainer.H:66
This functor does binary nuclear fusions on a single cell. Particles of the two reacting species are ...
Definition: NuclearFusionFunc.H:37
NuclearFusionFunc(const std::string &collision_name, MultiParticleContainer const *const mypc, const bool isSameSpecies)
Constructor of the NuclearFusionFunc class.
Definition: NuclearFusionFunc.H:59
NuclearFusionFunc()=default
Default constructor of the NuclearFusionFunc class.
Executor const & executor() const
Definition: NuclearFusionFunc.H:238
amrex::ParticleReal m_fusion_multiplier
Definition: NuclearFusionFunc.H:243
bool m_isSameSpecies
Definition: NuclearFusionFunc.H:253
WarpXParticleContainer::ParticleType ParticleType
Definition: NuclearFusionFunc.H:39
ParticleBins::index_type index_type
Definition: NuclearFusionFunc.H:43
amrex::ParticleReal m_probability_threshold
Definition: NuclearFusionFunc.H:250
Executor m_exe
Definition: NuclearFusionFunc.H:255
amrex::ParticleReal m_probability_target_value
Definition: NuclearFusionFunc.H:251
NuclearFusionType m_fusion_type
Definition: NuclearFusionFunc.H:252
ParticleTile< ParticleType, NArrayReal, NArrayInt, Allocator > ParticleTileType
T_ParticleType ParticleType
Definition: BinaryCollisionUtils.cpp:18
NuclearFusionType get_nuclear_fusion_type(const std::string &collision_name, MultiParticleContainer const *const mypc)
Definition: BinaryCollisionUtils.cpp:40
AMREX_GPU_HOST_DEVICE constexpr AMREX_FORCE_INLINE const T & max(const T &a, const T &b) noexcept
AMREX_GPU_HOST_DEVICE constexpr AMREX_FORCE_INLINE const T & min(const T &a, const T &b) noexcept
float dt
Definition: stencil.py:442
int queryWithParser(const amrex::ParmParse &a_pp, char const *const str, T &val)
Definition: ParserUtils.H:137
Definition: NuclearFusionFunc.H:88
NuclearFusionType m_fusion_type
Definition: NuclearFusionFunc.H:234
amrex::ParticleReal m_probability_target_value
Definition: NuclearFusionFunc.H:233
bool m_isSameSpecies
Definition: NuclearFusionFunc.H:235
amrex::ParticleReal m_probability_threshold
Definition: NuclearFusionFunc.H:232
AMREX_GPU_HOST_DEVICE AMREX_INLINE void operator()(index_type const I1s, index_type const I1e, index_type const I2s, index_type const I2e, index_type const *AMREX_RESTRICT I1, index_type const *AMREX_RESTRICT I2, const SoaData_type &soa_1, const SoaData_type &soa_2, GetParticlePosition< PIdx >, GetParticlePosition< PIdx >, amrex::ParticleReal const, amrex::ParticleReal const, amrex::ParticleReal const m1, amrex::ParticleReal const m2, amrex::Real const dt, amrex::Real const dV, index_type coll_idx, index_type const cell_start_pair, index_type *AMREX_RESTRICT p_mask, index_type *AMREX_RESTRICT p_pair_indices_1, index_type *AMREX_RESTRICT p_pair_indices_2, amrex::ParticleReal *AMREX_RESTRICT p_pair_reaction_weight, amrex::RandomEngine const &engine) const
Executor of the NuclearFusionFunc class. Performs nuclear fusions at the cell level using the algorit...
Definition: NuclearFusionFunc.H:127
amrex::ParticleReal m_fusion_multiplier
Definition: NuclearFusionFunc.H:231
@ theta
RZ needs all three position components.
Definition: NamedComponentParticleContainer.H:36
@ uz
Definition: NamedComponentParticleContainer.H:34
@ w
weight
Definition: NamedComponentParticleContainer.H:33
@ uy
Definition: NamedComponentParticleContainer.H:34
@ ux
Definition: NamedComponentParticleContainer.H:34
GpuArray< ParticleReal *, NAR > m_rdata
ParticleTileData< StorageParticleType, NArrayReal, NArrayInt > ParticleTileDataType