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 Shannon entropy mesh without dimension can be assigned

Browse source
This commit is contained in:
Paul Romano 2019-10-01 10:31:29 -05:00
parent ae1209e283
commit 9a78836d22
2 changed files with 4 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

3
src/mesh.cpp
View file

@ -155,7 +155,8 @@ RegularMesh::RegularMesh(pugi::xml_node node)
fatal_error("Must specify either <upper_right> and <width> on a mesh.");
}
if (shape_.dimension() > 0) {
// TODO: Change to zero when xtensor is updated
if (shape_.size() > 1) {
if (shape_.size() != lower_left_.size()) {
fatal_error("Number of entries on <lower_left> must be the same "
"as the number of entries on <dimension>.");

3
src/settings.cpp
View file

@ -512,7 +512,8 @@ void read_settings_xml()
if (!m) fatal_error("Only regular meshes can be used as an entropy mesh");
simulation::entropy_mesh = m;
if (m->shape_.dimension() == 0) {
// TODO: Change to zero when xtensor is updated
if (m->shape_.size() == 1) {
// If the user did not specify how many mesh cells are to be used in
// each direction, we automatically determine an appropriate number of
// cells