From 7650ea42d2bd253a20ae05c07eb64ce085836a20 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Tue, 22 Jan 2019 06:59:46 -0600 Subject: [PATCH] Fix threshold index check for photoelectric cross sections --- src/physics.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/physics.cpp b/src/physics.cpp index 817577cab..90f0708f3 100644 --- a/src/physics.cpp +++ b/src/physics.cpp @@ -303,7 +303,7 @@ void sample_photon_reaction(Particle* p) // Check threshold of reaction int i_start = shell.threshold; - if (i_grid <= i_start) continue; + if (i_grid < i_start) continue; // Evaluation subshell photoionization cross section double xs = std::exp(shell.cross_section(i_grid - i_start) +