From 58d10d6a3036c14e4fc57fde2c30fe49a44e9f68 Mon Sep 17 00:00:00 2001 From: John Tramm Date: Mon, 16 Dec 2019 20:49:16 +0000 Subject: [PATCH] increase nubank size from 10 to 15. Some nu values were getting be 10 or 11. Also added an assert to enforce this. --- include/openmc/particle.h | 2 +- src/physics.cpp | 2 ++ 2 files changed, 3 insertions(+), 1 deletion(-) diff --git a/include/openmc/particle.h b/include/openmc/particle.h index 134ebc105..00d37162b 100644 --- a/include/openmc/particle.h +++ b/include/openmc/particle.h @@ -329,7 +329,7 @@ public: std::vector> tracks_; - NuBank nu_bank_[10]; + NuBank nu_bank_[15]; }; } // namespace openmc diff --git a/src/physics.cpp b/src/physics.cpp index 8d0af575c..132c45bc9 100644 --- a/src/physics.cpp +++ b/src/physics.cpp @@ -169,6 +169,8 @@ create_fission_sites(Particle* p, int i_nuclide, const Reaction* rx, // First, if our bank is full then don't continue if (nu == 0) return; + assert( nu < 15 ); + // Initialize the counter of delayed neutrons encountered for each delayed // group. double nu_d[MAX_DELAYED_GROUPS] = {0.};