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

removed precision parameter

Browse source
This commit is contained in:
Ethan Peterson 2022-10-24 16:55:48 -04:00
parent fec53aa1fa
commit 33021f7af1
1 changed files with 9 additions and 28 deletions
Download patch file Download diff file Expand all files Collapse all files

37
openmc/geometry.py
View file

@ -436,18 +436,14 @@ class Geometry:
surfaces = cell.region.get_surfaces(surfaces)
return surfaces
def get_redundant_surfaces(self, precision=None):
"""Return all of the topologically redundant surface IDs
def get_redundant_surfaces(self):
"""Return all of the topologically redundant surface IDs.
Uses surface_precision attribute of Geometry instance for rounding and
comparing surface coefficients.
.. versionadded:: 0.12
Parameters
----------
precision : int, optional
The number of decimal places to round to for considering surface
coefficients to be equal. Defaults to the value specified by
Geometry.surface_precision.
Returns
-------
dict
@ -456,12 +452,9 @@ class Geometry:
that should replace it.
"""
precision = self.surface_precision if precision is None else precision
check_less_than('precision', precision, 16)
check_greater_than('precision', precision, 0)
tally = defaultdict(list)
for surf in self.get_all_surfaces().values():
coeffs = tuple(round(surf._coefficients[k], precision)
coeffs = tuple(round(surf._coefficients[k], self.surface_precision)
for k in surf._coeff_keys)
key = (surf._type,) + coeffs
tally[key].append(surf)
@ -615,23 +608,11 @@ class Geometry:
"""
return self._get_domains_by_name(name, case_sensitive, matching, 'lattice')
def remove_redundant_surfaces(self, precision=None):
"""Remove redundant surfaces from the geometry.
Parameters
----------
precision : int, optional
The number of decimal places to round to for considering surface
coefficients to be equal. Defaults to the value specified by
Geometry.surface_precision.
"""
precision = self.surface_precision if precision is None else precision
check_less_than('precision', precision, 16)
check_greater_than('precision', precision, 0)
def remove_redundant_surfaces(self):
"""Remove redundant surfaces from the geometry."""
# Get redundant surfaces
redundant_surfaces = self.get_redundant_surfaces(precision=precision)
redundant_surfaces = self.get_redundant_surfaces()
# Iterate through all cells contained in the geometry
for cell in self.get_all_cells().values():