import os import openmc from openmc.examples import random_ray_three_region_cube from tests.testing_harness import WeightWindowPyAPITestHarness class MGXSTestHarness(WeightWindowPyAPITestHarness): def _cleanup(self): super()._cleanup() f = 'mgxs.h5' if os.path.exists(f): os.remove(f) def test_random_ray_adjoint_fixed_source(): model = random_ray_three_region_cube() ww_mesh = openmc.RegularMesh() n = 6 width = 30.0 ww_mesh.dimension = (n, n, n) ww_mesh.lower_left = (0.0, 0.0, 0.0) ww_mesh.upper_right = (width, width, width) wwg = openmc.WeightWindowGenerator( method="fw_cadis", mesh=ww_mesh, max_realizations=model.settings.batches) model.settings.weight_window_generators = wwg model.settings.random_ray['volume_estimator'] = 'naive' harness = MGXSTestHarness('statepoint.10.h5', model) harness.main()