crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Fix bug in summary.py (assuming root universe is 0)

Browse source
This commit is contained in:
Paul Romano 2017-03-06 19:53:09 -06:00
parent 2d1897a051
commit b19a4c5ec4
1 changed files with 13 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

15
openmc/summary.py
View file

@ -163,6 +163,11 @@ class Summary(object):
def _finalize_geometry(self, cells, cell_fills, universes, lattices):
materials = {m.id: m for m in self.materials}
# Keep track of universes that are used as fills. That way, we can
# determine which universe is NOT used as a fill (and hence is the root
# universe)
fill_univ_ids = set()
# Iterate over all Cells and add fill Materials, Universes and Lattices
for cell_id, (fill_type, fill_id) in cell_fills.items():
# Retrieve the object corresponding to the fill type and ID
@ -174,14 +179,20 @@ class Summary(object):
fill = materials[fill_id] if fill_id > 0 else None
elif fill_type == 'universe':
fill = universes[fill_id]
fill_univ_ids.add(fill_id)
else:
fill = lattices[fill_id]
for idx in fill._natural_indices:
univ = fill.get_universe(idx)
fill_univ_ids.add(univ.id)
# Set the fill for the Cell
cells[cell_id].fill = fill
# Set the root universe for the Geometry
self.geometry.root_universe = universes[0]
# Determine root universe for geometry
non_fill = set(universes.keys()) - fill_univ_ids
self.geometry.root_universe = universes[non_fill.pop()]
def add_volume_information(self, volume_calc):
"""Add volume information to the geometry within the summary file