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

Make sure using multiple S(a,b) tables doesn't break depletion

Browse source
This commit is contained in:
Paul Romano 2020-07-08 15:40:42 -05:00
parent c6c4f3a091
commit d426b3a24e
1 changed files with 8 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

8
src/material.cpp
View file

@ -5,6 +5,7 @@
#include <iterator>
#include <string>
#include <sstream>
#include <unordered_set>
#include "xtensor/xbuilder.hpp"
#include "xtensor/xoperation.hpp"
@ -412,7 +413,14 @@ void Material::init_thermal()
{
std::vector<ThermalTable> tables;
std::unordered_set<int> already_checked;
for (const auto& table : thermal_tables_) {
// Make sure each S(a,b) table only gets checked once
if (already_checked.find(table.index_table) != already_checked.end()) {
continue;
}
already_checked.insert(table.index_table);
// In order to know which nuclide the S(a,b) table applies to, we need
// to search through the list of nuclides for one which has a matching
// name