Check for zero fission cross section before sampling fission reaction (…

…openmc-dev#2600)
stchaker · Oct 25, 2023 · 138eb5c · 138eb5c
1 parent 7b5faff
commit 138eb5c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/physics.cpp b/src/physics.cpp
@@ -106,7 +106,7 @@ void sample_neutron_reaction(Particle& p)
 
   const auto& nuc {data::nuclides[i_nuclide]};
 
-  if (nuc->fissionable_) {
+  if (nuc->fissionable_ && p.neutron_xs(i_nuclide).fission > 0.0) {
     auto& rx = sample_fission(i_nuclide, p);
     if (settings::run_mode == RunMode::EIGENVALUE) {
       create_fission_sites(p, i_nuclide, rx);