diff options
author | 2021-08-19 18:08:05 -0700 | |
---|---|---|
committer | 2021-08-19 18:08:05 -0700 | |
commit | 7a0a1afb81dfb2e3bac59a50f7820c00352ebc5c (patch) | |
tree | e815699a8c986eb9e46b14c1c98882e235f58587 /Source/Particles/ElementaryProcess | |
parent | 263d9627bb8038227ba1a5f40ee937b0affbf8d6 (diff) | |
download | WarpX-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.H | 2 |
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 ); |