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

Fix get_decay_modes ('unknown' wasn't recognized)

Browse source
This commit is contained in:
Paul Romano 2021-08-27 11:08:42 -05:00
parent 209fdd86ed
commit c08d9bc142
2 changed files with 17 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

7
openmc/data/decay.py
View file

@ -57,8 +57,11 @@ def get_decay_modes(value):
List of successive decays, e.g. ('beta-', 'neutron')
"""
return [_DECAY_MODES[int(x)][0] for x in
str(value).strip('0').replace('.', '')]
if int(value) == 10:
return ['unknown']
else:
return [_DECAY_MODES[int(x)][0] for x in
str(value).strip('0').replace('.', '')]
class FissionProductYields(EqualityMixin):

12
tests/unit_tests/test_data_decay.py
View file

@ -31,6 +31,18 @@ def u235_yields():
return openmc.data.FissionProductYields.from_endf(filename)
def test_get_decay_modes():
assert openmc.data.get_decay_modes(1.0) == ['beta-']
assert openmc.data.get_decay_modes(6.0) == ['sf']
assert openmc.data.get_decay_modes(10.0) == ['unknown']
assert openmc.data.get_decay_modes(1.5) == ['beta-', 'n']
assert openmc.data.get_decay_modes(1.4) == ['beta-', 'alpha']
assert openmc.data.get_decay_modes(1.55) == ['beta-', 'n', 'n']
assert openmc.data.get_decay_modes(1.555) == ['beta-', 'n', 'n', 'n']
assert openmc.data.get_decay_modes(2.4) == ['ec/beta+', 'alpha']
def test_nb90_halflife(nb90):
ufloat_close(nb90.half_life, ufloat(52560.0, 180.0))
ufloat_close(nb90.decay_constant, log(2.)/nb90.half_life)