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

Merge pull request #2310 from joshmay1/uni_nuclide_densities

Browse source 
Fix bug in `openmc.Universe.get_nuclide_densities()`
This commit is contained in:
Paul Romano 2022-11-29 19:49:45 -08:00 • committed by GitHub
commit d7fc8e59a2
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
2 changed files with 12 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

2
openmc/universe.py
View file

@ -469,7 +469,7 @@ class Universe(UniverseBase):
"""
nuclides = OrderedDict()
if self._atoms is not None:
if self._atoms:
volume = self.volume
for name, atoms in self._atoms.items():
nuclide = openmc.Nuclide(name)

11
tests/unit_tests/test_universe.py
View file

@ -132,3 +132,14 @@ def test_create_xml(cell_with_lattice):
assert all(c.get('universe') == str(u.id) for c in cell_elems)
assert not (set(c.get('id') for c in cell_elems) ^
set(str(c.id) for c in cells))
def test_get_nuclide_densities():
surf = openmc.Sphere()
material = openmc.Material()
material.add_elements_from_formula("H2O")
material.set_density("g/cm3", 1)
cell = openmc.Cell(region=-surf, fill=material)
universe = openmc.Universe(cells=[cell])
with pytest.raises(RuntimeError):
universe.get_nuclide_densities()