Merge pull request #1333 from compas-dev/use-tolerance

Use tolerance

Author: tomvanmele

CHANGELOG.md

@@ -29,6 +29,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 * Fixed broken import in `copas.geometry.trimesh_smoothing_numpy`.
 * Changed `RhinoBrep.trimmed` to return single result or raise `BrepTrimmingError` instead of returning a list.
 * Changed order of imports according to `isort` and changed line length to `179`.
+* Changed use of `compas.geometry.allclose` to `compas.tolerance.TOL.is_allclose`.
+* Changed use of `compas.geometry.close` to `compas.tolerance.TOL.is_close`.
 
 ### Removed
 

src/compas/geometry/_core/tangent.py

@@ -22,13 +22,13 @@ def tangent_points_to_circle_xy(circle, point):
 
     Examples
     --------
-    >>> from compas.geometry import allclose
+    >>> from compas.tolerance import TOL
     >>> circle = ((0, 0, 0), (0, 0, 1)), 1.0
     >>> point = (2, 4, 0)
     >>> t1, t2 = tangent_points_to_circle_xy(circle, point)
-    >>> allclose(t1, [-0.772, 0.636, 0.000], 1e-3)
+    >>> TOL.is_allclose(t1, [-0.772, 0.636, 0.000], atol=1e-3)
     True
-    >>> allclose(t2, [0.972, -0.236, 0.000], 1e-3)
+    >>> TOL.is_allclose(t2, [0.972, -0.236, 0.000], atol=1e-3)
     True
 
     """

src/compas/geometry/_core/transformations.py

@@ -62,13 +62,13 @@ def orthonormalize_axes(xaxis, yaxis):
 
     Examples
     --------
-    >>> from compas.geometry import allclose
+    >>> from compas.tolerance import TOL
     >>> xaxis = [1, 4, 5]
     >>> yaxis = [1, 0, -2]
     >>> xaxis, yaxis = orthonormalize_axes(xaxis, yaxis)
-    >>> allclose(xaxis, [0.1543, 0.6172, 0.7715], tol=0.001)
+    >>> TOL.is_allclose(xaxis, [0.1543, 0.6172, 0.7715], atol=0.001)
     True
-    >>> allclose(yaxis, [0.6929, 0.4891, -0.5298], tol=0.001)
+    >>> TOL.is_allclose(yaxis, [0.6929, 0.4891, -0.5298], atol=0.001)
     True
 
     """

src/compas/geometry/curves/arc.py

@@ -10,7 +10,7 @@
 from compas.geometry import Frame
 from compas.geometry import Point
 from compas.geometry import Vector
-from compas.geometry import close
+from compas.tolerance import TOL
 
 from .curve import Curve
 
@@ -244,7 +244,7 @@ def circumference(self):
 
     @property
     def is_circle(self):
-        return close(abs(abs(self.angle) - PI2), 0.0)
+        return TOL.is_close(abs(self.angle), PI2)
 
     @property
     def is_closed(self):

src/compas/geometry/curves/polyline.py

@@ -5,11 +5,11 @@
 from compas.geometry import Frame
 from compas.geometry import Line
 from compas.geometry import Point
-from compas.geometry import allclose
 from compas.geometry import is_point_on_line
 from compas.geometry import is_point_on_polyline
 from compas.geometry import transform_points
 from compas.itertools import pairwise
+from compas.tolerance import TOL
 
 from .curve import Curve
 
@@ -118,7 +118,7 @@ def __len__(self):
     def __eq__(self, other):
         if not hasattr(other, "__iter__") or not hasattr(other, "__len__") or len(self) != len(other):
             return False
-        return allclose(self, other)
+        return TOL.is_allclose(self, other)
 
     # ==========================================================================
     # properties

src/compas/geometry/intersections.py

@@ -609,15 +609,15 @@ def intersection_sphere_line(sphere, line):
 
     Examples
     --------
-    >>> from compas.geometry import allclose
+    >>> from compas.tolerance import TOL
 
     >>> sphere = (3.0, 7.0, 4.0), 10.0
     >>> line = (1.0, 0, 0.5), (2.0, 1.0, 0.5)
     >>> x1, x2 = intersection_sphere_line(sphere, line)
 
-    >>> allclose(x1, [11.634, 10.634, 0.500], 1e-3)
+    >>> TOL.is_allclose(x1, [11.634, 10.634, 0.500], atol=1e-3)
     True
-    >>> allclose(x2, [-0.634, -1.634, 0.50], 1e-3)
+    >>> TOL.is_allclose(x2, [-0.634, -1.634, 0.50], atol=1e-3)
     True
 
     """

src/compas/geometry/polygon.py

@@ -10,7 +10,6 @@
 from compas.geometry import Plane
 from compas.geometry import Point
 from compas.geometry import Transformation
-from compas.geometry import allclose
 from compas.geometry import area_polygon
 from compas.geometry import bounding_box
 from compas.geometry import centroid_polygon
@@ -19,6 +18,7 @@
 from compas.geometry import is_polygon_convex
 from compas.geometry import transform_points
 from compas.itertools import pairwise
+from compas.tolerance import TOL
 
 
 class Polygon(Geometry):
@@ -119,7 +119,7 @@ def __iter__(self):
     def __eq__(self, other):
         if not hasattr(other, "__iter__") or not hasattr(other, "__len__") or len(self) != len(other):
             return False
-        return allclose(self, other)
+        return TOL.is_allclose(self, other)
 
     # ==========================================================================
     # Properties

src/compas/geometry/projection.py

@@ -12,12 +12,12 @@
 """
 
 from compas.geometry import Transformation
-from compas.geometry import allclose
 from compas.geometry import decompose_matrix
 from compas.geometry import matrix_from_orthogonal_projection
 from compas.geometry import matrix_from_parallel_projection
 from compas.geometry import matrix_from_perspective_entries
 from compas.geometry import matrix_from_perspective_projection
+from compas.tolerance import TOL
 from compas.utilities import flatten
 
 
@@ -47,7 +47,7 @@ class Projection(Transformation):
     def __init__(self, matrix=None, check=False, name=None):
         if matrix and check:
             _, _, _, _, perspective = decompose_matrix(matrix)
-            if not allclose(flatten(matrix), flatten(matrix_from_perspective_entries(perspective))):
+            if not TOL.is_allclose(flatten(matrix), flatten(matrix_from_perspective_entries(perspective))):
                 raise ValueError("This is not a proper projection matrix.")
         super(Projection, self).__init__(matrix=matrix, name=name)
 

src/compas/geometry/quaternion.py

@@ -292,12 +292,12 @@ def from_frame(cls, frame):
 
         Examples
         --------
-        >>> from compas.geometry import allclose
+        >>> from compas.tolerance import TOL
         >>> from compas.geometry import Frame
         >>> q = [1., -2., 3., -4.]
         >>> F = Frame.from_quaternion(q)
         >>> Q = Quaternion.from_frame(F)
-        >>> allclose(Q.canonized(), quaternion_canonize(quaternion_unitize(q)))
+        >>> TOL.is_allclose(Q.canonized(), quaternion_canonize(quaternion_unitize(q)))
         True
 
         """

src/compas/geometry/reflection.py

@@ -12,13 +12,13 @@
 """
 
 from compas.geometry import Transformation
-from compas.geometry import allclose
 from compas.geometry import cross_vectors
 from compas.geometry import decompose_matrix
 from compas.geometry import dot_vectors
 from compas.geometry import identity_matrix
 from compas.geometry import matrix_from_perspective_entries
 from compas.geometry import normalize_vector
+from compas.tolerance import TOL
 from compas.utilities import flatten
 
 
@@ -50,7 +50,7 @@ class Reflection(Transformation):
     def __init__(self, matrix=None, check=False, name=None):
         if matrix and check:
             _, _, _, _, perspective = decompose_matrix(matrix)
-            if not allclose(flatten(matrix), flatten(matrix_from_perspective_entries(perspective))):
+            if not TOL.is_allclose(flatten(matrix), flatten(matrix_from_perspective_entries(perspective))):
                 raise ValueError("This is not a proper reflection matrix.")
         super(Reflection, self).__init__(matrix=matrix, name=name)