topologicpy 0.8.9__py3-none-any.whl → 0.8.11__py3-none-any.whl

topologicpy/BVH.py

@@ -282,7 +282,7 @@ class BVH:
             # Add an edge from the parent to this vertex (if a parent exists)
             if parent_vertex is not None:
                 d = Vertex.Distance(parent_vertex, current_vertex)
-                if d < tolerance:
+                if d <= tolerance:
                     current_vertex = Topology.Translate(current_vertex, tolerance*random.uniform(2,50), tolerance*random.uniform(2,50), tolerance*random.uniform(2,50))
                 edge = Edge.ByVertices(parent_vertex, current_vertex, tolerance=tolerance)
                 graph = Graph.AddEdge(graph, edge, silent=True)

topologicpy/Cell.py

@@ -998,14 +998,14 @@ class Cell():
                 verticalFaces.append(aFace)
                 verticalApertures += getApertures(aFace)
             elif aCode == 1:
-                if abs(Vertex.Z(Topology.Centroid(aFace)) - zMin) < tolerance:
+                if abs(Vertex.Z(Topology.Centroid(aFace)) - zMin) <= tolerance:
                     bottomHorizontalFaces.append(aFace)
                     bottomHorizontalApertures += getApertures(aFace)
                 else:
                     topHorizontalFaces.append(aFace)
                     topHorizontalApertures += getApertures(aFace)
             elif aCode == 2:
-                if abs(Vertex.Z(Topology.Centroid(aFace)) - zMax) < tolerance:
+                if abs(Vertex.Z(Topology.Centroid(aFace)) - zMax) <= tolerance:
                     topHorizontalFaces.append(aFace)
                     topHorizontalApertures += getApertures(aFace)
                 else:

topologicpy/CellComplex.py

@@ -546,7 +546,7 @@ class CellComplex():
                     internalVerticalApertures += getApertures(aFace)
             elif aCode == 1:
                 if n == 1:
-                    if abs(Vertex.Z(Topology.Centroid(aFace)) - zMin) < tolerance:
+                    if abs(Vertex.Z(Topology.Centroid(aFace)) - zMin) <= tolerance:
                         bottomHorizontalFaces.append(aFace)
                         bottomHorizontalApertures += getApertures(aFace)
                     else:
@@ -557,7 +557,7 @@ class CellComplex():
                     internalHorizontalApertures += getApertures(aFace)
             elif aCode == 2:
                 if n == 1:
-                    if abs(Vertex.Z(Topology.Centroid(aFace)) - zMax) < tolerance:
+                    if abs(Vertex.Z(Topology.Centroid(aFace)) - zMax) <= tolerance:
                         topHorizontalFaces.append(aFace)
                         topHorizontalApertures += getApertures(aFace)
                     else:

topologicpy/Edge.py

@@ -153,15 +153,15 @@ class Edge():
             
             shared_vertex = None
             
-            if Vertex.Distance(start1, start2) < tolerance:
+            if Vertex.Distance(start1, start2) <= tolerance:
                 shared_vertex = start1
-            elif Vertex.Distance(start1, end2) < tolerance:
+            elif Vertex.Distance(start1, end2) <= tolerance:
                 shared_vertex = start1
                 edge2 = Edge.Reverse(edge2)
-            elif Vertex.Distance(end1, start2) < tolerance:
+            elif Vertex.Distance(end1, start2) <= tolerance:
                 shared_vertex = start2
                 edge1 = Edge.Reverse(edge1)
-            elif Vertex.Distance(end1, end2) < tolerance:
+            elif Vertex.Distance(end1, end2) <= tolerance:
                 shared_vertex = end1
                 edge1 = Edge.Reverse(edge1)
                 edge2 = Edge.Reverse(edge2)
@@ -176,10 +176,10 @@ class Edge():
         if not Topology.IsInstance(edgeB, "Edge"):
             print("Edge.Bisect - Error: The input edgeB parameter is not a valid topologic edge. Returning None.")
             return None
-        if Edge.Length(edgeA) < tolerance:
+        if Edge.Length(edgeA) <= tolerance:
             print("Edge.Bisect - Error: The input edgeA parameter is shorter than the input tolerance parameter. Returning None.")
             return None
-        if Edge.Length(edgeB) < tolerance:
+        if Edge.Length(edgeB) <= tolerance:
             print("Edge.Bisect - Error: The input edgeB parameter is shorter than the input tolerance parameter. Returning None.")
             return None
         
@@ -313,7 +313,7 @@ class Edge():
             if not silent:
                 print("Edge.ByStartVertexEndVertex - Error: The input vertexA and vertexB parameters are the same vertex. Returning None.")
             return None
-        if Vertex.Distance(vertexA, vertexB) < tolerance:
+        if Vertex.Distance(vertexA, vertexB) <= tolerance:
             if not silent:
                 print("Edge.ByStartVertexEndVertex - Error: The distance between the input vertexA and vertexB parameters is less than the input tolerance. Returning None.")
             return None
@@ -360,7 +360,7 @@ class Edge():
                 print("Edge.ByVertexDirectionLength - Error: The input vertex parameter is not a valid vertex. Returning None.")
             return None
         
-        if length < tolerance:
+        if length <= tolerance:
             if not silent:
                 print("Edge.ByVertexDirectionLength - Error: The input edge parameter must not be less than the input tolerance parameter. Returning None.")
             return None
@@ -646,7 +646,7 @@ class Edge():
             print("Edge.Extend - Error: The input edge parameter is not a valid topologic edge. Returning None.")
             return None
         distance = abs(distance)
-        if distance < tolerance:
+        if distance <= tolerance:
             return edge
         sv = Edge.StartVertex(edge)
         ev = Edge.EndVertex(edge)
@@ -810,11 +810,11 @@ class Edge():
                 evb = Edge.EndVertex(edges[i])
                 dsvsv = Vertex.Distance(sva, svb)
                 devev = Vertex.Distance(eva, evb)
-                if dsvsv < tolerance and devev < tolerance:
+                if dsvsv <= tolerance and devev <= tolerance:
                     return i
                 dsvev = Vertex.Distance(sva, evb)
                 devsv = Vertex.Distance(eva, svb)
-                if dsvev < tolerance and devsv < tolerance:
+                if dsvev <= tolerance and devsv <= tolerance:
                     return i
         return None
 
@@ -856,7 +856,7 @@ class Edge():
             distances.append(Vertex.Distance(pair[0], pair[1]))
         v_list = Helper.Sort(v_list, distances)
         closest_pair = v_list[0]
-        if Vertex.Distance(closest_pair[0], closest_pair[1]) < tolerance:
+        if Vertex.Distance(closest_pair[0], closest_pair[1]) <= tolerance:
             return Topology.Centroid(Cluster.ByTopologies(closest_pair))
         
         if Edge.IsCollinear(edgeA, edgeB, tolerance=tolerance):
@@ -915,11 +915,11 @@ class Edge():
         if not Topology.IsInstance(edgeB, "Edge"):
             print("Edge.IsCollinear - Error: The input parameter edgeB is not a valid edge. Returning None")
             return None
-        if Edge.Length(edgeA) < tolerance:
-            print("Edge.IsCollinear - Error: The length of edgeA is less than the tolerance. Returning None")
+        if Edge.Length(edgeA) <= tolerance:
+            print("Edge.IsCollinear - Error: The length of edgeA is less than or equal the tolerance. Returning None")
             return None
-        if Edge.Length(edgeB) < tolerance:
-            print("Edge.IsCollinear - Error: The length of edgeB is less than the tolerance. Returning None")
+        if Edge.Length(edgeB) <= tolerance:
+            print("Edge.IsCollinear - Error: The length of edgeB is less than or equal to the tolerance. Returning None")
             return None
         
         # Get start and end points of the first edge
@@ -962,7 +962,7 @@ class Edge():
         distance_end = distance_from_line(end_b, start_a_coords, direction_a)
 
         # Check if both distances are within tolerance
-        return bool(distance_start < tolerance) and bool(distance_end < tolerance)
+        return bool(distance_start <= tolerance) and bool(distance_end <= tolerance)
     
     @staticmethod
     def IsCoplanar(edgeA, edgeB, mantissa: int = 6, tolerance: float = 0.0001):
@@ -996,11 +996,11 @@ class Edge():
         if not Topology.IsInstance(edgeB, "Edge"):
             print("Edge.IsCoplanar - Error: The input parameter edgeB is not a valid edge. Returning None")
             return None
-        if Edge.Length(edgeA) < tolerance:
-            print("Edge.IsCoplanar - Error: The length of edgeA is less than the tolerance. Returning None")
+        if Edge.Length(edgeA) <= tolerance:
+            print("Edge.IsCoplanar - Error: The length of edgeA is less than or equal to the tolerance. Returning None")
             return None
-        if Edge.Length(edgeB) < tolerance:
-            print("Edge.IsCoplanar - Error: The length of edgeB is less than the tolerance. Returning None")
+        if Edge.Length(edgeB) <= tolerance:
+            print("Edge.IsCoplanar - Error: The length of edgeB is less than or equal to the tolerance. Returning None")
             return None
         
         # Extract points
@@ -1513,8 +1513,8 @@ class Edge():
         distance = abs(distance)
         if distance == 0:
             return edge
-        if distance < tolerance:
-            print("Edge.Trim - Warning: The input distance parameter is less than the input tolerance parameter. Returning the input edge.")
+        if distance <= tolerance:
+            print("Edge.Trim - Warning: The input distance parameter is less than or equal to the input tolerance parameter. Returning the input edge.")
             return edge
         sv = Edge.StartVertex(edge)
         ev = Edge.EndVertex(edge)

topologicpy/Face.py

@@ -325,7 +325,7 @@ class Face():
                 print("Face.ByOffset - Warning: The input face parameter is not a valid face. Returning None.")
             return None
         
-        if abs(Face.Normal(face)[2] + 1) < tolerance:
+        if abs(Face.Normal(face)[2] + 1) <= tolerance:
             reverse = not(reverse)
         eb = Face.Wire(face)
         
@@ -494,7 +494,7 @@ class Face():
                 # We want this difference to be as close to 0 as possible
                 loss = (new_area - area) ** 2
                 # If the loss is less than the tolerance, accept the result and return a loss of 0.
-                if loss < tolerance:
+                if loss <= tolerance:
                     return 0
                 # Otherwise, return the actual loss value.
                 return loss 
@@ -1622,7 +1622,7 @@ class Face():
             dV = direction
         uV = Vector.Normalize(dV)
         dot = sum([i*j for (i, j) in zip(uV, faceNormal)])
-        if dot < tolerance:
+        if dot <= tolerance:
             return False
         return True
 
@@ -2278,7 +2278,7 @@ class Face():
         def vertex_part_of_face(vertex, face, tolerance=0.0001):
             vertices = Topology.Vertices(face)
             for v in vertices:
-                if Vertex.Distance(vertex, v) < tolerance:
+                if Vertex.Distance(vertex, v) <= tolerance:
                     return True
             return False
 
@@ -3040,9 +3040,9 @@ class Face():
             if not silent:
                 print("Face.NormalEdge - Error: The input face parameter is not a valid face. Retuning None.")
             return None
-        if length < tolerance:
+        if length <= tolerance:
             if not silent:
-                print("Face.NormalEdge - Error: The input length parameter is less than the input tolerance. Retuning None.")
+                print("Face.NormalEdge - Error: The input length parameter is less than or equal to the input tolerance. Retuning None.")
             return None
         iv = Face.InternalVertex(face)
         u, v = Face.VertexParameters(face, iv)
@@ -3411,6 +3411,61 @@ class Face():
             return None
         return Face.ByWire(wire, tolerance=tolerance)
 
+
+    @staticmethod
+    def ThirdVertex(face, tolerance: float = 0.0001, silent: bool = False):
+        """
+        Returns a third vertex on the input face to enable rotation matrix creation.
+
+        Parameters
+        ----------
+        face : topologic_core.Face
+            The input face.
+        tolerance : float , optional
+            The desired tolerance. The default is 0.0001.
+        silent : bool , optional
+            If set to True, no error and warning messages are printed. Otherwise, they are. The default is False.
+
+        Returns
+        -------
+        topologic_core.Face
+            The created face.
+
+        """
+        from topologicpy.Vector import Vector
+        from topologicpy.Vertex import Vertex
+        from topologicpy.Topology import Topology
+
+        if not Topology.IsInstance(face, "Face"):
+            if not silent:
+                print("Face.ThirdVertex - Error: The input face parameter is not a valid face. Returning None.")
+            return None
+        # Retrieve all vertices of the face
+        vertices = Face.Vertices(face)
+        centroid = Topology.Centroid(face)
+        normal = Face.Normal(face)
+        for vertex in vertices:
+            # Skip the centroid itself
+            if Vertex.Distance(centroid, vertex) <= tolerance:
+                continue
+            
+            # Vector from the centroid to the current vertex
+            vector_to_vertex = Vector.ByVertices(centroid, vertex)
+            vector_to_vertex_normalized = Vector.Normalize(vector_to_vertex)
+
+
+            # Check if the vector_to_vertex is collinear with the normal direction
+            if Vector.IsCollinear(vector_to_vertex_normalized, normal, tolerance):
+                continue
+
+            # If not collinear, return this vertex
+            return vertex
+
+        # No valid third vertex found
+        if not silent:
+                print("Face.ThirdVertex - Warning: No valid third vertex could be found. Returning None.")
+        return None
+
     @staticmethod
     def Trapezoid(origin= None, widthA: float = 1.0, widthB: float = 0.75, offsetA: float = 0.0, offsetB: float = 0.0, length: float = 1.0, direction: list = [0, 0, 1], placement: str = "center", tolerance: float = 0.0001):
         """