crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Add check on user-specified bounding boxes for volume calculations

Browse source
This commit is contained in:
Paul Romano 2016-08-15 17:50:03 -05:00
parent e0525d9425
commit 6c29c39443
1 changed files with 18 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

19
openmc/volume.py
View file

@ -1,6 +1,7 @@
from collections import Iterable, Mapping
from numbers import Real, Integral
from xml.etree import ElementTree as ET
from warnings import warn
import numpy as np
import pandas as pd
@ -68,10 +69,26 @@ class VolumeCalculation(object):
self.samples = samples
if lower_left is not None:
self.lower_left = lower_left
if upper_right is None:
raise ValueError('Both lower-left and upper-right coordinates '
'should be specified')
# For cell domains, try to compute bounding box and make sure
# user-specified one is valid
if self.domain_type == 'cell':
for c in domains:
if c.region is None:
continue
ll, ur = c.region.bounding_box
if np.any(np.isinf(ll)) or np.any(np.isinf(ur)):
continue
if (np.any(np.asarray(lower_left) > ll) or
np.any(np.asarray(upper_right) < ur)):
warn("Specified bounding box is smaller than computed "
"bounding box for cell {}. Volume calculation may "
"be incorrect!".format(c.id))
self.lower_left = lower_left
self.upper_right = upper_right
else:
if self.domain_type == 'cell':