# 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.
"""Tests the Interface."""
import unittest
from armi import interfaces
from armi import settings
[docs]class DummyInterface(interfaces.Interface):
name = "Dummy"
function = "dummyAction"
[docs]class TestCodeInterface(unittest.TestCase):
"""Test Code interface."""
def setUp(self):
self.cs = settings.Settings()
[docs] def test_isRequestedDetailPoint(self):
"""Tests notification of detail points."""
newSettings = {"dumpSnapshot": ["000001", "995190"]}
cs = self.cs.modified(newSettings=newSettings)
i = DummyInterface(None, cs)
self.assertEqual(i.isRequestedDetailPoint(0, 1), True)
self.assertEqual(i.isRequestedDetailPoint(995, 190), True)
self.assertEqual(i.isRequestedDetailPoint(5, 10), False)
[docs] def test_enabled(self):
"""Test turning interfaces on and off."""
i = DummyInterface(None, self.cs)
self.assertEqual(i.enabled(), True)
i.enabled(False)
self.assertEqual(i.enabled(), False)
i.enabled(True)
self.assertEqual(i.enabled(), True)
[docs] def test_nameContains(self):
i = DummyInterface(None, self.cs)
self.assertFalse(i.nameContains("nope"))
self.assertTrue(i.nameContains("Dum"))
[docs] def test_distributable(self):
i = DummyInterface(None, self.cs)
self.assertEqual(i.distributable(), 1)
[docs] def test_preDistributeState(self):
i = DummyInterface(None, self.cs)
self.assertEqual(i.preDistributeState(), {})
[docs] def test_duplicate(self):
i = DummyInterface(None, self.cs)
iDup = i.duplicate()
self.assertEqual(type(i), type(iDup))
self.assertEqual(i.enabled(), iDup.enabled())
[docs]class TestTightCoupler(unittest.TestCase):
"""Test the tight coupler class."""
def setUp(self):
cs = settings.Settings()
cs["tightCoupling"] = True
cs["tightCouplingSettings"] = {
"dummyAction": {"parameter": "nothing", "convergence": 1.0e-5}
}
self.interface = DummyInterface(None, cs)
[docs] def test_couplerActive(self):
self.assertIsNotNone(self.interface.coupler)
[docs] def test_storePreviousIterationValue(self):
self.interface.coupler.storePreviousIterationValue(1.0)
self.assertEqual(self.interface.coupler._previousIterationValue, 1.0)
[docs] def test_storePreviousIterationValueException(self):
with self.assertRaises(TypeError) as cm:
self.interface.coupler.storePreviousIterationValue({5.0})
the_exception = cm.exception
self.assertEqual(the_exception.error_code, 3)
[docs] def test_isConvergedValueError(self):
with self.assertRaises(ValueError) as cm:
self.interface.coupler.isConverged(1.0)
the_exception = cm.exception
self.assertEqual(the_exception.error_code, 3)
[docs] def test_isConverged(self):
"""Ensure TightCoupler.isConverged() works with float, 1D list, and ragged 2D list.
.. test:: The tight coupling logic is based around a convergence criteria.
:id: T_ARMI_OPERATOR_PHYSICS1
:tests: R_ARMI_OPERATOR_PHYSICS
Notes
-----
2D lists can end up being ragged as assemblies can have different number of blocks.
Ragged lists are easier to manage with lists as opposed to numpy.arrays,
namely, their dimension is preserved.
"""
# show a situation where it doesn't converge
previousValues = {
"float": 1.0,
"list1D": [1.0, 2.0],
"list2D": [[1, 2, 3], [1, 2]],
}
updatedValues = {
"float": 5.0,
"list1D": [5.0, 6.0],
"list2D": [[5, 6, 7], [5, 6]],
}
for previous, current in zip(previousValues.values(), updatedValues.values()):
self.interface.coupler.storePreviousIterationValue(previous)
self.assertFalse(self.interface.coupler.isConverged(current))
# show a situation where it DOES converge
previousValues = updatedValues
for previous, current in zip(previousValues.values(), updatedValues.values()):
self.interface.coupler.storePreviousIterationValue(previous)
self.assertTrue(self.interface.coupler.isConverged(current))
[docs] def test_isConvergedRuntimeError(self):
"""Test to ensure 3D arrays do not work."""
previous = [[[1, 2, 3]], [[1, 2, 3]], [[1, 2, 3]]]
updatedValues = [[[5, 6, 7]], [[5, 6, 7]], [[5, 6, 7]]]
self.interface.coupler.storePreviousIterationValue(previous)
with self.assertRaises(RuntimeError) as cm:
self.interface.coupler.isConverged(updatedValues)
the_exception = cm.exception
self.assertEqual(the_exception.error_code, 3)
[docs] def test_getListDimension(self):
a = [1, 2, 3]
self.assertEqual(interfaces.TightCoupler.getListDimension(a), 1)
a = [[1, 2, 3]]
self.assertEqual(interfaces.TightCoupler.getListDimension(a), 2)
a = [[[1, 2, 3]]]
self.assertEqual(interfaces.TightCoupler.getListDimension(a), 3)