crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

numpy errstate required for serr

Browse source
This commit is contained in:
Shikhar Kumar 2019-06-04 20:13:21 -04:00
parent 5e8f50dd2e
commit 0bfbd2c823
1 changed files with 3 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

5
openmc/cmfd.py
View file

@ -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