crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Make sure E_out is not undefined in ContinuousTabular::sample

Browse source
This commit is contained in:
Paul Romano 2020-07-02 21:18:20 -05:00
parent 28df499464
commit 1af42dbb7d
1 changed files with 1 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/distribution_energy.cpp
View file

@ -224,13 +224,11 @@ double ContinuousTabular::sample(double E, uint64_t* seed) const
double E_l_k = distribution_[l].e_out[k];
double p_l_k = distribution_[l].p[k];
double E_out;
double E_out = E_l_k;
if (distribution_[l].interpolation == Interpolation::histogram) {
// Histogram interpolation
if (p_l_k > 0.0 && k >= n_discrete) {
E_out = E_l_k + (r1 - c_k)/p_l_k;
} else {
E_out = E_l_k;
}
} else if (distribution_[l].interpolation == Interpolation::lin_lin) {