mirror of
https://github.com/openmc-dev/openmc.git
synced 2026-07-27 13:45:36 -04:00
34 lines
1.1 KiB
Python
34 lines
1.1 KiB
Python
|
import glob
|
||
|
|
import os
|
||
|
|
|
||
|
|
from openmc import StatePoint
|
||
|
|
|
||
|
|
from tests.testing_harness import TestHarness
|
||
|
|
|
||
|
|
|
||
|
|
class EntropyTestHarness(TestHarness):
|
||
|
|
def _get_results(self):
|
||
|
|
"""Digest info in the statepoint and return as a string."""
|
||
|
|
# Read the statepoint file.
|
||
|
|
statepoint = glob.glob(os.path.join(os.getcwd(), self._sp_name))[0]
|
||
|
|
with StatePoint(statepoint) as sp:
|
||
|
|
# Write out k-combined.
|
||
|
|
outstr = 'k-combined:\n'
|
||
|
|
outstr += '{:12.6E} {:12.6E}\n'.format(sp.keff.n, sp.keff.s)
|
||
|
|
|
||
|
|
# Write out entropy data.
|
||
|
|
outstr += 'entropy:\n'
|
||
|
|
results = ['{:12.6E}'.format(x) for x in sp.entropy]
|
||
|
|
outstr += '\n'.join(results) + '\n'
|
||
|
|
|
||
|
|
return outstr
|
||
|
|
|
||
|
|
'''
|
||
|
|
# This test is adapted from "Monte Carlo power iteration: Entropy and spatial correlations,"
|
||
|
|
M. Nowak et al. The cross sections are defined explicitly so that the value for entropy
|
||
|
|
is exactly 9 and the eigenvalue is exactly 1.
|
||
|
|
'''
|
||
|
|
def test_entropy():
|
||
|
|
harness = EntropyTestHarness('statepoint.10.h5')
|
||
|
|
harness.main()
|