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

Avoid out-of-bounds access on index_inelastic_scatter_

Browse source
This commit is contained in:
Paul Romano 2023-02-06 21:49:43 -06:00
parent 6218becb17
commit 5470ec75b4
1 changed files with 2 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/physics.cpp
View file

@ -705,17 +705,10 @@ void scatter(Particle& p, int i_nuclide)
// =======================================================================
// INELASTIC SCATTERING
int j = 0;
int n = nuc->index_inelastic_scatter_.size();
int i = 0;
while (prob < cutoff) {
for (int j = 0; j < n && prob < cutoff; ++j) {
i = nuc->index_inelastic_scatter_[j];
++j;
// Check to make sure inelastic scattering reaction sampled
if (i >= nuc->reactions_.size()) {
p.write_restart();
fatal_error("Did not sample any reaction for nuclide " + nuc->name_);
}
// add to cumulative probability
prob += nuc->reactions_[i]->xs(micro);