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

Update trigger_statepoint_restart test results

Browse source
This commit is contained in:
Paul Romano 2021-04-19 15:20:31 -05:00
parent 9552be7316
commit bfac737001
3 changed files with 16 additions and 17 deletions
Download patch file Download diff file Expand all files Collapse all files

8
tests/regression_tests/trigger_statepoint_restart/inputs_true.dat
View file

@ -13,11 +13,11 @@
<?xml version='1.0' encoding='utf-8'?>
<settings>
<run_mode>eigenvalue</run_mode>
<particles>200</particles>
<batches>10</batches>
<inactive>5</inactive>
<particles>400</particles>
<batches>15</batches>
<inactive>10</inactive>
<keff_trigger>
<threshold>0.004</threshold>
<threshold>0.003</threshold>
<type>std_dev</type>
</keff_trigger>
<trigger>

6
tests/regression_tests/trigger_statepoint_restart/results_true.dat
View file

@ -1,5 +1,5 @@
k-combined:
2.916922E-01 3.293799E-03
2.963805E-01 2.986351E-03
tally 1:
6.184423E+01
4.789617E+02
8.017826E+01
6.436316E+02

19
tests/regression_tests/trigger_statepoint_restart/test.py
View file

@ -25,22 +25,21 @@ def model():
# Settings
settings = openmc.Settings()
settings.run_mode = 'eigenvalue'
settings.batches = 10
settings.inactive = 5
settings.particles = 200
# Choose a sufficiently low threshold to trigger after more than 10 batches.
# 0.004 seems to take 13 batches.
settings.keff_trigger = {'type': 'std_dev', 'threshold': 0.004}
settings.batches = 15
settings.inactive = 10
settings.particles = 400
# Choose a sufficiently low threshold to enable use of trigger
settings.keff_trigger = {'type': 'std_dev', 'threshold': 0.003}
settings.trigger_max_batches = 1000
settings.trigger_batch_interval = 1
settings.trigger_active = True
settings.verbosity = 1 # to test that this works even with no output
# Tallies
t = openmc.Tally()
t.scores = ['flux']
tallies = openmc.Tallies([t])
# Put it all together
model = openmc.model.Model(materials=materials,
geometry=geometry,
@ -92,7 +91,7 @@ class TriggerStatepointRestartTestHarness(PyAPITestHarness):
with openmc.StatePoint(spfile) as sp:
sp_batchno_1 = sp.current_batch
k_combined_1 = sp.k_combined
assert sp_batchno_1 > 10
assert sp_batchno_1 > 5
print('Last batch no = %d' % sp_batchno_1)
self._write_inputs(self._get_inputs())
self._compare_inputs()
@ -110,7 +109,7 @@ class TriggerStatepointRestartTestHarness(PyAPITestHarness):
with openmc.StatePoint(spfile) as sp:
sp_batchno_2 = sp.current_batch
k_combined_2 = sp.k_combined
assert sp_batchno_2 > 10
assert sp_batchno_2 > 5
assert sp_batchno_1 == sp_batchno_2, \
'Different final batch number after restart'
# need str() here as uncertainties.ufloat instances are always different