From 2ae943c37c279337cd5667dcea5ee609f6c703e5 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Wed, 22 Jun 2016 13:07:34 +0700 Subject: [PATCH] Make sure densities are repeated when calculating microscopic cross sections for Legendre moments --- openmc/mgxs/mgxs.py | 1 + 1 file changed, 1 insertion(+) diff --git a/openmc/mgxs/mgxs.py b/openmc/mgxs/mgxs.py index 088db649fa..c12c61bdf4 100644 --- a/openmc/mgxs/mgxs.py +++ b/openmc/mgxs/mgxs.py @@ -1517,6 +1517,7 @@ class MGXS(object): densities = self.get_nuclide_densities(nuclides) else: densities = self.get_nuclide_densities('sum') + densities = np.repeat(densities, len(self.rxn_rate_tally.scores)) tile_factor = df.shape[0] / len(densities) df['mean'] /= np.tile(densities, tile_factor) df['std. dev.'] /= np.tile(densities, tile_factor)