Fix threshold index check for photoelectric cross sections

This commit is contained in:
Paul Romano 2019-01-22 06:59:46 -06:00
parent 450359ac87
commit 7650ea42d2

View file

@ -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) +