From 596eab4bb69743a85e1f03830229eedb1091957d Mon Sep 17 00:00:00 2001 From: Shikhar Kumar Date: Mon, 22 Jul 2019 12:08:04 -0400 Subject: [PATCH] Define default value for resnb when flux equals 0 --- openmc/cmfd.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/openmc/cmfd.py b/openmc/cmfd.py index 742cdeabfd..ca52502c3a 100644 --- a/openmc/cmfd.py +++ b/openmc/cmfd.py @@ -1141,7 +1141,7 @@ class CMFDRun(object): self._dhat = np.zeros((nx, ny, nz, ng, 6)) # Set reference diffusion parameters - if self._ref_d: + if list(self._ref_d): self._set_reference_params = True # Check length of reference diffusion parameters equal to number of # energy groups @@ -2214,7 +2214,8 @@ class CMFDRun(object): res = leakage + interactions - scattering - (1.0 / keff) * fission # Normalize res by flux and bank res - self._resnb = np.divide(res, self._flux, where=self._flux > 0) + self._resnb = np.divide(res, self._flux, where=self._flux > 0, + out=np.zeros_like(self._flux)) # Calculate RMS and record for this batch self._balance.append(np.sqrt(