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

updated path for dagmc file

Browse source
This commit is contained in:
Jonathan Shimwell 2022-08-12 10:43:40 +01:00
parent b1ea00b6cb
commit 623b867033
1 changed files with 31 additions and 26 deletions
Download patch file Download diff file Expand all files Collapse all files

57
tests/unit_tests/test_dagmc_universe.py
View file

@ -1,50 +1,58 @@
import openmc
import pytest
from pathlib import Path
def test_bounding_box():
"""Checks that the DAGMCUniverse.bounding_box returns the correct values"""
u = openmc.DAGMCUniverse("tests/regression_tests/dagmc/universes/dagmc.h5m")
u = openmc.DAGMCUniverse(str(Path(__file__).parent.resolve() / "dagmc/dagmc.h5m"))
ll, ur = u.bounding_box
assert ll == pytest.approx((-25., -25., -25))
assert ur == pytest.approx((25., 25., 25))
assert ll == pytest.approx((-25.0, -25.0, -25))
assert ur == pytest.approx((25.0, 25.0, 25))
def test_bounding_region():
"""Checks that the DAGMCUniverse.bounding_region() returns a region with
correct surfaces and boundary types"""
u = openmc.DAGMCUniverse("tests/regression_tests/dagmc/universes/dagmc.h5m")
u = openmc.DAGMCUniverse(str(Path(__file__).parent.resolve() / "dagmc/dagmc.h5m"))
region = u.bounding_region() # should default to bounded_type='box'
assert isinstance(region, openmc.Region)
assert len(region) == 6
assert region[0].surface.type == 'x-plane'
assert region[1].surface.type == 'x-plane'
assert region[2].surface.type == 'y-plane'
assert region[3].surface.type == 'y-plane'
assert region[4].surface.type == 'z-plane'
assert region[5].surface.type == 'z-plane'
assert region[0].surface.boundary_type == 'vacuum'
assert region[1].surface.boundary_type == 'vacuum'
assert region[2].surface.boundary_type == 'vacuum'
assert region[3].surface.boundary_type == 'vacuum'
assert region[4].surface.boundary_type == 'vacuum'
assert region[5].surface.boundary_type == 'vacuum'
assert region[0].surface.type == "x-plane"
assert region[1].surface.type == "x-plane"
assert region[2].surface.type == "y-plane"
assert region[3].surface.type == "y-plane"
assert region[4].surface.type == "z-plane"
assert region[5].surface.type == "z-plane"
assert region[0].surface.boundary_type == "vacuum"
assert region[1].surface.boundary_type == "vacuum"
assert region[2].surface.boundary_type == "vacuum"
assert region[3].surface.boundary_type == "vacuum"
assert region[4].surface.boundary_type == "vacuum"
assert region[5].surface.boundary_type == "vacuum"
region = u.bounding_region(bounded_type='sphere', boundary_type='reflective')
region = u.bounding_region(bounded_type="sphere", boundary_type="reflective")
assert isinstance(region, openmc.Region)
assert isinstance(region, openmc.Halfspace)
assert region.surface.type == 'sphere'
assert region.surface.boundary_type == 'reflective'
assert region.surface.type == "sphere"
assert region.surface.boundary_type == "reflective"
def test_bounded_universe():
"""Checks that the DAGMCUniverse.bounded_universe() returns a
openmc.Universe with correct surface ids and cell ids"""
u = openmc.DAGMCUniverse("tests/regression_tests/dagmc/universes/dagmc.h5m")
u = openmc.DAGMCUniverse(str(Path(__file__).parent.resolve() / "dagmc/dagmc.h5m"))
# bounded with defaults
bu = u.bounded_universe()
cells = list(bu.get_all_cells().items())
assert isinstance(bu, openmc.Universe)
assert len(cells) == 1
assert cells[0][0] == 10000 # default bounding_cell_id is 10000
assert cells[0][1].id == 10000 # default bounding_cell_id is 10000
@ -53,16 +61,13 @@ def test_bounded_universe():
assert surfaces[0][1].id == 10000
# bounded with non defaults
bu = u.bounded_universe(
bounding_cell_id=42,
bounded_type='sphere',
starting_id=43
)
bu = u.bounded_universe(bounding_cell_id=42, bounded_type="sphere", starting_id=43)
cells = list(bu.get_all_cells().items())
assert isinstance(bu, openmc.Universe)
assert len(cells) == 1
assert cells[0][0] == 42 # default bounding_cell_id is 10000
assert cells[0][1].id == 42 # default bounding_cell_id is 10000
surfaces = list(cells[0][1].region.get_surfaces().items())
assert surfaces[0][1].type == 'sphere'
assert surfaces[0][1].type == "sphere"
assert surfaces[0][1].id == 43