Shapely preprocessor

sectionproperties/pre/pre.py

@@ -96,8 +96,6 @@ def create_mesh(
     :return: Object containing generated mesh data
     :rtype: :class:`meshpy.triangle.MeshInfo`
     """
-    print("Pre.mesh beginning")
-
     # check_geometry(points, facets, holes, control_points, atol=atol)
     if not isinstance(mesh_sizes, list):
         mesh_sizes = [mesh_sizes]
@@ -114,7 +112,6 @@ def create_mesh(
         mesh.regions[i] = [cp[0], cp[1], region_id, mesh_sizes[i]]
         region_id += 1
 
-    print("Just before mesh")
     mesh = triangle.build(
         mesh,
         min_angle=30,

sectionproperties/pre/sections.py

@@ -56,7 +56,7 @@ def __init__(
         self.holes = []
         self.perimeter = []
         self._recovery_points = []
-        self.tol = 12 # Represents num of decimal places of precision
+        self.tol = 12 # Represents num of decimal places of precision for point locations
         self.compile_geometry()
 
     def _repr_svg_(self):
@@ -205,7 +205,6 @@ def create_mesh(self, mesh_sizes: Union[float, List[float]]):
 
             Mesh generated from the above geometry.
         """
-        print("Geometry create mesh method begin")
         if isinstance(mesh_sizes, (list, tuple)) and len(mesh_sizes) == 1:
             mesh_size = mesh_sizes[0]
         elif isinstance(mesh_sizes, (float, int)):
@@ -242,8 +241,9 @@ def align_to(
         :return: Geometry object translated to alignment location
         :rtype: :class:`sections.pre.sections.Geometry`
         """
-        # Mappings are indexes in the list of bbox extents of both
-        # self and 'align_to'
+        # Mappings are for indexes in the list of bbox extents of both
+        # 'self' and 'align_to'. i.e. a mapping of which "word" corresponds
+        # to which bounding box coordinate
         align_self_map = {
             "left": 0,
             "right": 1,
@@ -305,12 +305,17 @@ def align_center(self, align_to: Optional[Geometry] = None):
         :rtype: :class:`sections.pre.sections.Geometry`
         """
         cx, cy = list(self.geom.centroid.coords)[0]
+        # Suggested by Agent 6-6-6: Hard-rounding of cx and cy allows
+        # for greater precision in placing geometry with its centroid
+        # near [0, 0]. True [0, 0] placement will not be possible due
+        # to floating point errors.
         if align_to is None:
-            shift_x, shift_y = -cx, -cy
+            shift_x, shift_y = round(-cx, self.tol), round(-cy, self.tol)
+
         else:
             align_cx, align_cy = list(align_to.geom.centroid.coords)[0]
-            shift_x = align_cx - cx
-            shift_y = align_cy - cy
+            shift_x = round(align_cx - cx, self.tol)
+            shift_y = round(align_cy - cy, self.tol)
         new_geom = self.shift_section(x_offset=shift_x, y_offset=shift_y)
         return new_geom
 
@@ -854,7 +859,7 @@ def __init__(self, geoms: Union[MultiPolygon, List[Geometry]]):
         self.holes = []
         self.perimeter = []
         self.compile_geometry()
-
+        self.tol = 12
         self.mesh = None
 
     def _repr_svg_(self):
@@ -975,10 +980,8 @@ def create_mesh(self, mesh_sizes: List[float]):
 
             Mesh generated from the above geometry.
         """
-        print("Begin create mesh method")
         if len(mesh_sizes) == 1:
             mesh_sizes = mesh_sizes * len(self.control_points)
-        print("Call pre.create_mesh")
         self.mesh = pre.create_mesh(
             self.points, self.facets, self.holes, self.control_points, mesh_sizes
         )

sectionproperties/tests/test_validation.py

@@ -45,14 +45,14 @@ def test_angle_all_properties():
     check.almost_equal(angle_section.section_props.sf_xx_plus, 1.788366)
     check.almost_equal(angle_section.section_props.sf_xx_minus, 9.207672e-1)
     check.almost_equal(angle_section.section_props.sf_yy_plus, 1.829943) # Altered from 1.829944
-    check.almost_equal(angle_section.section_props.sf_yy_minus, 5.646721e-1) # Altered from 5.646723
+    check.almost_equal(angle_section.section_props.sf_yy_minus, 5.646721e-1) # Altered from 5.646723e-1
     check.almost_equal(angle_section.section_props.s11, 1.210275e5)
     check.almost_equal(angle_section.section_props.s22, 4.376054e4)
     check.almost_equal(angle_section.section_props.sf_11_plus, 1.238049)
     check.almost_equal(angle_section.section_props.sf_11_minus, 1.744103)
     check.almost_equal(angle_section.section_props.sf_22_plus, 1.564994)
     check.almost_equal(angle_section.section_props.sf_22_minus, 2.107303)
-    check.almost_equal(angle_section.section_props.j, 1.354663e5)
+    check.almost_equal(angle_section.section_props.j, 1.354663e5) # Altered from 1.354663e5
     check.almost_equal(angle_section.section_props.gamma, 162220735.49)
     check.almost_equal(angle_section.section_props.A_s11, 8.855951e2)
     check.almost_equal(angle_section.section_props.A_s22, 1.460240e3)