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

Expanding on eff checks for log-log case.

Browse source
This commit is contained in:
Patrick Shriwise 2022-07-30 00:43:20 -05:00
parent c09f539cd5
commit aeca1173e7
1 changed files with 26 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

35
tests/regression_tests/filter_energyfun/test.py
View file

@ -72,17 +72,34 @@ class FilterEnergyFunHarness(PyAPITestHarness):
# Read the statepoint file.
sp = openmc.StatePoint(self._sp_name)
# check that information is round-tripped correctly from
# the statepoint file
sp_tally = sp.get_tally(id=3)
sp_filt = sp_tally.find_filter(openmc.EnergyFunctionFilter)
# statepoint file round-trip checks
# linear-linear interpolation tally
sp_lin_lin_tally = sp.get_tally(id=2)
sp_lin_lin_filt = sp_lin_lin_tally.find_filter(openmc.EnergyFunctionFilter)
model_tally = self._model.tallies[2]
model_filt = model_tally.find_filter(openmc.EnergyFunctionFilter)
model_lin_lin_tally = self._model.tallies[1]
model_lin_lin_filt = model_lin_lin_tally.find_filter(openmc.EnergyFunctionFilter)
assert sp_filt.interpolation == 'log-log'
assert all(sp_filt.energy == model_filt.energy)
assert all(sp_filt.y == model_filt.y)
assert sp_lin_lin_filt.interpolation == 'linear-linear'
assert all(sp_lin_lin_filt.energy == model_lin_lin_filt.energy)
assert all(sp_lin_lin_filt.y == model_lin_lin_filt.y)
# log-log interpolation tally
sp_log_log_tally = sp.get_tally(id=3)
sp_log_log_filt = sp_log_log_tally.find_filter(openmc.EnergyFunctionFilter)
model_log_log_tally = self._model.tallies[2]
model_log_log_filt = model_log_log_tally.find_filter(openmc.EnergyFunctionFilter)
assert sp_log_log_filt.interpolation == 'log-log'
assert all(sp_log_log_filt.energy == model_log_log_filt.energy)
assert all(sp_log_log_filt.y == model_log_log_filt.y)
# because the values of y are monotonically increasing,
# we expect the log-log tally to have a higher value
assert all(sp_lin_lin_tally.mean < sp_log_log_tally.mean)
def test_filter_energyfun(model):