mirror of
https://github.com/openmc-dev/openmc.git
synced 2026-07-29 06:35:48 -04:00
37 lines
1,007 B
Python
37 lines
1,007 B
Python
"""Regression tests for openmc.deplete.integrator.cecm algorithm.
|
|
|
|
These tests integrate a simple test problem described in dummy_geometry.py.
|
|
"""
|
|
|
|
from pytest import approx
|
|
import openmc.deplete
|
|
|
|
from tests import dummy_operator
|
|
|
|
|
|
def test_cecm(run_in_tmpdir):
|
|
"""Integral regression test of integrator algorithm using CE/CM."""
|
|
|
|
op = dummy_operator.DummyOperator()
|
|
op.output_dir = "test_integrator_regression"
|
|
|
|
# Perform simulation using the MCNPX/MCNP6 algorithm
|
|
dt = [0.75, 0.75]
|
|
power = 1.0
|
|
openmc.deplete.cecm(op, dt, power, print_out=False)
|
|
|
|
# Load the files
|
|
res = openmc.deplete.ResultsList(op.output_dir / "depletion_results.h5")
|
|
|
|
_, y1 = res.get_atoms("1", "1")
|
|
_, y2 = res.get_atoms("1", "2")
|
|
|
|
# Mathematica solution
|
|
s1 = [1.86872629872102, 1.395525772416039]
|
|
s2 = [2.18097439443550, 2.69429754646747]
|
|
|
|
assert y1[1] == approx(s1[0])
|
|
assert y2[1] == approx(s1[1])
|
|
|
|
assert y1[2] == approx(s2[0])
|
|
assert y2[2] == approx(s2[1])
|