aboutsummaryrefslogtreecommitdiff
path: root/Source/Particles/ElementaryProcess
diff options
context:
space:
mode:
authorGravatar Edoardo Zoni <59625522+EZoni@users.noreply.github.com> 2021-08-19 18:08:05 -0700
committerGravatar GitHub <noreply@github.com> 2021-08-19 18:08:05 -0700
commit7a0a1afb81dfb2e3bac59a50f7820c00352ebc5c (patch)
treee815699a8c986eb9e46b14c1c98882e235f58587 /Source/Particles/ElementaryProcess
parent263d9627bb8038227ba1a5f40ee937b0affbf8d6 (diff)
downloadWarpX-7a0a1afb81dfb2e3bac59a50f7820c00352ebc5c.tar.gz
WarpX-7a0a1afb81dfb2e3bac59a50f7820c00352ebc5c.tar.zst
WarpX-7a0a1afb81dfb2e3bac59a50f7820c00352ebc5c.zip
Ionization.H: Fix Bug (Division by Zero) (#2214)
Diffstat (limited to 'Source/Particles/ElementaryProcess')
-rw-r--r--Source/Particles/ElementaryProcess/Ionization.H2
1 files changed, 1 insertions, 1 deletions
diff --git a/Source/Particles/ElementaryProcess/Ionization.H b/Source/Particles/ElementaryProcess/Ionization.H
index bb919550a..c970fb4f9 100644
--- a/Source/Particles/ElementaryProcess/Ionization.H
+++ b/Source/Particles/ElementaryProcess/Ionization.H
@@ -124,7 +124,7 @@ struct IonizationFilterFunc
);
// Compute probability of ionization p
- amrex::Real w_dtau = 1._rt/ ga * m_adk_prefactor[ion_lev] *
+ amrex::Real w_dtau = (E == 0._rt) ? 0._rt : 1._rt/ ga * m_adk_prefactor[ion_lev] *
std::pow(E, m_adk_power[ion_lev]) *
std::exp( m_adk_exp_prefactor[ion_lev]/E );
amrex::Real p = 1._rt - std::exp( - w_dtau );