From 0bfbd2c823f2c574fd666f8bacfd77142a1b198f Mon Sep 17 00:00:00 2001 From: Shikhar Kumar Date: Tue, 4 Jun 2019 20:13:21 -0400 Subject: [PATCH] numpy errstate required for serr --- openmc/cmfd.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/openmc/cmfd.py b/openmc/cmfd.py index c7dd33155..76f6fe3b1 100644 --- a/openmc/cmfd.py +++ b/openmc/cmfd.py @@ -1855,8 +1855,9 @@ class CMFDRun(object): kerr = abs(k_o - k_n) / k_n # Calculate max error in source - serr = np.sqrt(np.sum(np.where(s_n > 0, ((s_n-s_o) / s_n)**2, 0)) - / len(s_n)) + with np.errstate(divide='ignore', invalid='ignore'): + serr = np.sqrt(np.sum(np.where(s_n > 0, ((s_n-s_o) / s_n)**2, 0)) + / len(s_n)) # Check for convergence iconv = kerr < self._cmfd_ktol and serr < self._stol