From 6650dcf833f0c3f2975f190e0f234dc28548b8b1 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Wed, 28 Apr 2021 08:02:47 -0500 Subject: [PATCH] Use nu_bank_ for fixed source mode too Right now, if an outgoing energy filter is used with a nu-fission score in fixed source mode, OpenMC will segfault because it tries to pull information on fission sites from nu_bank_, which is not currently populated. --- src/physics.cpp | 12 +++++------- src/physics_mg.cpp | 12 +++++------- 2 files changed, 10 insertions(+), 14 deletions(-) diff --git a/src/physics.cpp b/src/physics.cpp index 4d692abef2..7a8be4ee09 100644 --- a/src/physics.cpp +++ b/src/physics.cpp @@ -217,13 +217,11 @@ create_fission_sites(Particle& p, int i_nuclide, const Reaction& rx) } // Write fission particles to nuBank - if (use_fission_bank) { - p.nu_bank_.emplace_back(); - Particle::NuBank* nu_bank_entry = &p.nu_bank_.back(); - nu_bank_entry->wgt = site.wgt; - nu_bank_entry->E = site.E; - nu_bank_entry->delayed_group = site.delayed_group; - } + p.nu_bank_.emplace_back(); + Particle::NuBank* nu_bank_entry = &p.nu_bank_.back(); + nu_bank_entry->wgt = site.wgt; + nu_bank_entry->E = site.E; + nu_bank_entry->delayed_group = site.delayed_group; } // If shared fission bank was full, and no fissions could be added, diff --git a/src/physics_mg.cpp b/src/physics_mg.cpp index b8e8fea454..b2dd4078ea 100644 --- a/src/physics_mg.cpp +++ b/src/physics_mg.cpp @@ -176,13 +176,11 @@ create_fission_sites(Particle& p) } // Write fission particles to nuBank - if (use_fission_bank) { - p.nu_bank_.emplace_back(); - Particle::NuBank* nu_bank_entry = &p.nu_bank_.back(); - nu_bank_entry->wgt = site.wgt; - nu_bank_entry->E = site.E; - nu_bank_entry->delayed_group = site.delayed_group; - } + p.nu_bank_.emplace_back(); + Particle::NuBank* nu_bank_entry = &p.nu_bank_.back(); + nu_bank_entry->wgt = site.wgt; + nu_bank_entry->E = site.E; + nu_bank_entry->delayed_group = site.delayed_group; } // If shared fission bank was full, and no fissions could be added,