def invoke(self):
from armi.nuclearDataIO import isotxs
from armi.nuclearDataIO import xsLibraries
runLog.setVerbosity(0)
refIsotxs = isotxs.readBinary(self.args.reference)
for fname in self.args.comparisonFiles:
cmpIsotxs = isotxs.readBinary(fname)
xsLibraries.compare(refIsotxs, cmpIsotxs)
def invoke(self):
from armi.nuclearDataIO import isotxs
from armi.nuclearDataIO import xsLibraries
runLog.setVerbosity(0)
refIsotxs = isotxs.readBinary(self.args.reference)
# flux = read_broad_group_flux(refIsotxs.neutronEnergyGroups)
for fname in self.args.comparisonFiles:
cmpIsotxs = isotxs.readBinary(fname)
xsLibraries.compare(refIsotxs, cmpIsotxs) # , flux)
def test_mergeXSLibrariesWithDifferentDataWorks(self):
if self.xsLibGenerationErrorStack is not None:
print(self.xsLibGenerationErrorStack)
raise Exception("see stdout for stack trace")
# check to make sure they labels overlap... or are actually the same
labels = set(self.xsLib.nuclideLabels)
self.assertEqual(labels, set(self.isotxsAA.nuclideLabels))
self.assertEqual(labels, set(self.gamisoAA.nuclideLabels))
self.assertEqual(labels, set(self.pmatrxAA.nuclideLabels))
# the whole thing is different from the sum of its components
self.assertFalse(xsLibraries.compare(self.xsLib, self.isotxsAA))
self.assertFalse(xsLibraries.compare(self.xsLib, self.gamisoAA))
self.assertFalse(xsLibraries.compare(self.xsLib, self.pmatrxAA))
# individual components are the same
self.assertTrue(isotxs.compare(self.xsLib, self.isotxsAA))
self.assertTrue(gamiso.compare(self.xsLib, self.gamisoAA))
self.assertTrue(pmatrx.compare(self.xsLib, self.pmatrxAA))
def test_mergeTwoXSLibFiles(self):
emptyXSLib = xsLibraries.IsotxsLibrary()
emptyXSLib.merge(self.libAA)
self.__class__.libAA = None
emptyXSLib.merge(self.libAB)
self.__class__.libAB = None
self.assertEqual(set(self.libCombined.nuclideLabels),
set(emptyXSLib.nuclideLabels))
self.assertTrue(xsLibraries.compare(emptyXSLib, self.libCombined))
self.getWriteFunc()(emptyXSLib, self.testFileName)
self.assertTrue(filecmp.cmp(self.getLibAA_ABPath(), self.testFileName))
def test_compareDifferentComponentsOfAnXSLibrary(self):
self.assertTrue(xsLibraries.compare(self.isotxsAA, self.isotxsAA))
self.assertTrue(xsLibraries.compare(self.pmatrxAA, self.pmatrxAA))
aa = isotxs.readBinary(ISOTXS_AA)
del aa[aa.nuclideLabels[0]]
self.assertFalse(xsLibraries.compare(aa, self.isotxsAA))
def test_canPickleAndUnpicklePMATRX(self):
pikAA = cPickle.loads(cPickle.dumps(self.pmatrxAA))
self.assertTrue(xsLibraries.compare(pikAA, self.pmatrxAA))
def test_canPickleAndUnpickleGAMISO(self):
pikAA = cPickle.loads(cPickle.dumps(self.gamisoAA))
self.assertTrue(xsLibraries.compare(pikAA, self.gamisoAA))
def test_canPickleAndUnpickleISOTXS(self):
pikAA = cPickle.loads(cPickle.dumps(self.isotxsAA))
self.assertTrue(xsLibraries.compare(pikAA, self.isotxsAA))
