# Copyright 2019 TerraPower, LLC
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
"""Test rtflux reading and writing."""
import os
import unittest
from armi.nuclearDataIO.cccc import rtflux
from armi.utils.directoryChangers import TemporaryDirectoryChanger
THIS_DIR = os.path.dirname(__file__)
# This rtflux was made by DIF3D 11 in a Cartesian test case.
SIMPLE_RTFLUX = os.path.join(THIS_DIR, "fixtures", "simple_cartesian.rtflux")
[docs]class Testrtflux(unittest.TestCase):
r"""
Tests the rtflux class.
"""
[docs] def test_readrtflux(self):
"""Ensure we can read a rtflux file."""
flux = rtflux.RtfluxStream.readBinary(SIMPLE_RTFLUX)
self.assertEqual(
flux.groupFluxes.shape,
(
flux.metadata["NINTI"],
flux.metadata["NINTJ"],
flux.metadata["NINTK"],
flux.metadata["NGROUP"],
),
)
[docs] def test_writertflux(self):
"""Ensure that we can write a modified rtflux file."""
with TemporaryDirectoryChanger():
flux = rtflux.RtfluxStream.readBinary(SIMPLE_RTFLUX)
# perturb off-diag item to check row/col ordering
flux.groupFluxes[2, 1, 3, 5] *= 1.1
flux.groupFluxes[1, 2, 4, 6] *= 1.2
rtflux.RtfluxStream.writeBinary(flux, "rtflux2")
flux2 = rtflux.RtfluxStream.readBinary("rtflux2")
self.assertAlmostEqual(
flux2.groupFluxes[2, 1, 3, 5], flux.groupFluxes[2, 1, 3, 5]
)
[docs] def test_rwAscii(self):
"""Ensure that we can read/write in ascii format."""
with TemporaryDirectoryChanger():
flux = rtflux.RtfluxStream.readBinary(SIMPLE_RTFLUX)
rtflux.RtfluxStream.writeAscii(flux, "rtflux.ascii")
flux2 = rtflux.RtfluxStream.readAscii("rtflux.ascii")
self.assertTrue((flux2.groupFluxes == flux.groupFluxes).all())
[docs] def test_adjoint(self):
"""Ensure adjoint reads energy groups differently."""
real = rtflux.RtfluxStream.readBinary(SIMPLE_RTFLUX)
adjoint = rtflux.AtfluxStream.readBinary(SIMPLE_RTFLUX)
self.assertFalse((real.groupFluxes == adjoint.groupFluxes).all())
g = 3
self.assertTrue(
(
real.groupFluxes[:, :, :, g]
== adjoint.groupFluxes[:, :, :, real.metadata["NGROUP"] - g - 1]
).all()
)