topologicpy 0.7.0__py3-none-any.whl → 0.7.3__py3-none-any.whl

topologicpy/Edge.py

@@ -597,9 +597,10 @@ class Edge():
 
 
     @staticmethod
-    def IsCollinear(edgeA, edgeB, mantissa: int = 6, angTolerance: float = 0.1, tolerance: float = 0.0001) -> bool:
+    def IsCollinear(edgeA, edgeB, mantissa: int = 6, tolerance: float = 0.0001) -> bool:
         """
         Return True if the two input edges are collinear. Returns False otherwise.
+        This code is based on a contribution by https://github.com/gaoxipeng
 
         Parameters
         ----------
@@ -609,8 +610,6 @@ class Edge():
             The second input edge.
         mantissa : int , optional
             The desired length of the mantissa. The default is 6.
-        angTolerance : float , optional
-            The angular tolerance used for the test. The default is 0.1.
         tolerance : float , optional
             The desired tolerance. The default is 0.0001.
 
@@ -622,25 +621,62 @@ class Edge():
         """
         from topologicpy.Vertex import Vertex
         from topologicpy.Topology import Topology
+        import numpy as np
 
         if not Topology.IsInstance(edgeA, "Edge"):
-            print("Edge.IsCollinear - Error: The input edgeA parameter is not a valid topologic edge. Returning None.")
+            print("Edge.IsCollinear - Error: The input parameter edgeA is not a valid edge. Returning None")
             return None
         if not Topology.IsInstance(edgeB, "Edge"):
-            print("Edge.IsCollinear - Error: The input edgeB parameter is not a valid topologic edge. Returning None.")
+            print("Edge.IsCollinear - Error: The input parameter edgeB is not a valid edge. Returning None")
             return None
-        ang = Edge.Angle(edgeA, edgeB, mantissa=mantissa, bracket=True)
-        svA = Edge.StartVertex(edgeA)
-        evA = Edge.EndVertex(edgeA)
-        svB = Edge.StartVertex(edgeB)
-        evB = Edge.EndVertex(edgeB)
-        d1 = Vertex.Distance(svA, svB)
-        d2 = Vertex.Distance(svA, evB)
-        d3 = Vertex.Distance(evA, svB)
-        d4 = Vertex.Distance(evA, evB)
-        if (d1 < tolerance or d2 < tolerance or d3 < tolerance or d4 < tolerance) and (abs(ang) < angTolerance or (abs(180 - ang) < angTolerance)):
-            return True
-        return False
+        if Edge.Length(edgeA) < tolerance:
+            print("Edge.IsCollinear - Error: The length of edgeA is less than 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")
+            return None
+        
+        # Get start and end points of the first edge
+        start_a = Edge.StartVertex(edgeA)
+        end_a = Edge.EndVertex(edgeA)
+        start_a_coords = np.array([Vertex.X(start_a), Vertex.Y(start_a), Vertex.Z(start_a)])
+        end_a_coords = np.array(
+            [Vertex.X(end_a, mantissa=mantissa), Vertex.Y(end_a, mantissa=mantissa), Vertex.Z(end_a, mantissa=mantissa)])
+
+        # Calculate the direction vector of the first edge
+        direction_a = end_a_coords - start_a_coords
+
+        # Normalize the direction vector
+        norm_a = np.linalg.norm(direction_a)
+        if norm_a == 0:
+            print("Edge.IsCollinear - Error: Division by zero. Returning None.")
+            return None
+        direction_a /= norm_a
+
+        # Function to calculate perpendicular distance from a point to the line defined by a point and direction vector
+        def distance_from_line(point, line_point, line_dir):
+            point = np.array([Vertex.X(point, mantissa=mantissa), Vertex.Y(point, mantissa=mantissa),
+                            Vertex.Z(point, mantissa=mantissa)])
+            line_point = np.array(line_point)
+            diff = point - line_point
+            cross_product = np.cross(diff, line_dir)
+            line_dir_norm = np.linalg.norm(line_dir)
+            if line_dir_norm == 0:
+                print("Edge.IsCollinear - Error: Division by zero. Returning None.")
+                return None
+            distance = np.linalg.norm(cross_product) / np.linalg.norm(line_dir)
+            return distance
+
+        # Get start and end points of the second edge
+        start_b = Edge.StartVertex(edgeB)
+        end_b = Edge.EndVertex(edgeB)
+
+        # Calculate distances for start and end vertices of the second edge to the line defined by the first edge
+        distance_start = distance_from_line(start_b, start_a_coords, direction_a)
+        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)
     
     @staticmethod
     def IsParallel(edgeA, edgeB, mantissa: int = 6, angTolerance: float = 0.1) -> bool:

topologicpy/Face.py

@@ -927,7 +927,19 @@ class Face():
             The external boundary of the input face.
 
         """
-        return face.ExternalBoundary()
+        from topologicpy.Vector import Vector
+        from topologicpy.Wire import Wire
+        from topologicpy.Topology import Topology
+
+        eb = face.ExternalBoundary()
+        f_dir = Face.Normal(face)
+        faceVertices = Topology.Vertices(eb)
+        temp_face = Face.ByWire(eb)
+        temp_dir = Face.Normal(temp_face)
+        if Vector.IsAntiParallel(f_dir, temp_dir):
+            faceVertices.reverse()
+        eb = Wire.ByVertices(faceVertices)
+        return eb
     
     @staticmethod
     def FacingToward(face, direction: list = [0,0,-1], asVertex: bool = False, tolerance: float = 0.0001) -> bool:
@@ -2313,7 +2325,7 @@ class Face():
             The external boundary of the input face.
 
         """
-        return face.ExternalBoundary()
+        return Face.ExternalBoundary(face)
     
     @staticmethod
     def Wires(face) -> list:

topologicpy/Topology.py

@@ -3199,6 +3199,146 @@ class Topology():
             return True
         return False
 
+
+    def ExportToDXF(topologies, path,  overwrite=False):
+        """
+        Exports the input topology to a DXF file. See https://en.wikipedia.org/wiki/AutoCAD_DXF.
+        THe DXF version is 'R2010'
+        This is experimental and only geometry is exported.
+
+        Parameters
+        ----------
+        topologies : list or topologic_core.Topology
+            The input list of topologies. This can also be a single topologic_core.Topology.
+        path : str
+            The input file path.
+        overwrite : bool , optional
+            If set to True the ouptut file will overwrite any pre-existing file. Otherwise, it won't. The default is False.
+
+        Returns
+        -------
+        bool
+            True if the export operation is successful. False otherwise.
+
+        """
+        import os
+        import warnings
+
+        try:
+            import ezdxf
+        except:
+            print("Topology.ExportToDXF - Information: Installing required ezdxf library.")
+            try:
+                os.system("pip install ezdxf")
+            except:
+                os.system("pip install ezdxf --user")
+            try:
+                import ezdxf
+                print("Topology.ExportToDXF - Information: ezdxf library installed successfully.")
+            except:
+                warnings.warn("Topology.ExportToDXF - Error: Could not import ezdxf library. Please install it manually. Returning None.")
+                return None
+        
+        from topologicpy.Vertex import Vertex
+        from topologicpy.Edge import Edge
+        from topologicpy.Cluster import Cluster
+        from topologicpy.Topology import Topology
+
+        from os.path import exists
+        if not isinstance(topologies, list):
+            topologies = [topologies]
+        topologies = [topology for topology in topologies if Topology.IsInstance(topology, "Topology")]
+        if len(topologies) < 1:
+            print("Topology.ExportToDXF - Error: The inupt list parameter topologies does not contain any valid topologies. Returning None.")
+        # Make sure the file extension is .brep
+        ext = path[len(path)-4:len(path)]
+        if ext.lower() != ".dxf":
+            path = path+".dxf"
+        if not overwrite and exists(path):
+            print("Topology.ExportToDXF - Error: a file already exists at the specified path and overwrite is set to False. Returning None.")
+            return None
+        
+        def add_vertices(vertices, msp):
+            for v in vertices:
+                if Topology.IsInstance(v, "Vertex"):
+                    msp.add_point((Vertex.X(v), Vertex.Y(v), Vertex.Z(v)))
+        def add_edges(edges, msp):
+            for e in edges:
+                if Topology.IsInstance(e, "Edge"):
+                    sv = Edge.StartVertex(e)
+                    ev = Edge.EndVertex(e)
+                    start = (Vertex.X(sv), Vertex.Y(sv), Vertex.Z(sv))
+                    end = (Vertex.X(ev), Vertex.Y(ev), Vertex.Z(ev))
+                    msp.add_line(start, end)
+        def add_wires(wires, msp):
+            for i, w in enumerate(wires):
+                if Topology.IsInstance(w, "Wire"):
+                    block_name = "Wire_"+str(i+1).zfill(3)
+                    block = doc.blocks.new(name=block_name)
+                    # Add edges to the block
+                    edges = Topology.Edges(w)
+                    for edge in edges:
+                        sv = Edge.StartVertex(edge)
+                        ev = Edge.EndVertex(edge)
+                        start = (Vertex.X(sv), Vertex.Y(sv), Vertex.Z(sv))
+                        end = (Vertex.X(ev), Vertex.Y(ev), Vertex.Z(ev))
+                        block.add_line(start, end)
+                    # Insert the block into the model space
+                    msp.add_blockref(block_name, insert=(0, 0, 0))
+        def add_meshes(meshes, msp):
+            for m in meshes:
+                data = Topology.Geometry(m)
+                vertices = data['vertices']
+                faces = data['faces']
+                mesh = msp.add_mesh()
+                mesh.dxf.subdivision_levels = 0
+                with mesh.edit_data() as mesh_data:
+                    mesh_data.vertices = vertices
+                    mesh_data.faces = faces
+        # Create a new DXF document
+        doc = ezdxf.new(dxfversion='R2010')
+        msp = doc.modelspace()
+        i = 1
+        for topology in topologies:
+            if Topology.IsInstance(topology, "Vertex"):
+                add_vertices([topology], msp)
+            elif Topology.IsInstance(topology, "Edge"):
+                add_edges([topology], msp)
+            elif Topology.IsInstance(topology, "Wire"):
+                add_wires([topology], msp)
+            elif Topology.IsInstance(topology, "Face"):
+                add_meshes([topology], msp)
+            elif Topology.IsInstance(topology, "Shell"):
+                add_meshes([topology], msp)
+            elif Topology.IsInstance(topology, "Cell"):
+                add_meshes([topology], msp)
+            elif Topology.IsInstance(topology, "CellComplex"):
+                add_meshes([topology], msp)
+            elif Topology.IsInstance(topology, "Cluster"):
+                cellComplexes = Topology.CellComplexes(topology)
+                add_meshes(cellComplexes, msp)
+                cells = Cluster.FreeCells(topology)
+                add_meshes(cells, msp)
+                shells = Cluster.FreeShells(topology)
+                add_meshes(shells, msp)
+                faces = Cluster.FreeFaces(topology)
+                add_meshes(faces, msp)
+                wires = Cluster.FreeWires(topology)
+                add_wires(wires, msp)
+                edges = Cluster.FreeEdges(topology)
+                add_edges(edges, msp)
+                vertices = Cluster.FreeVertices(topology)
+                add_vertices(vertices, msp)
+        
+        # Save the DXF document
+        status = False
+        try:
+            doc.saveas(path)
+            status = True
+        except:
+            status = False
+        return status
+
     '''
     @staticmethod
     def ExportToIPFS(topology, url, port, user, password):
@@ -3979,6 +4119,8 @@ class Topology():
 
         """
         from topologicpy.Vertex import Vertex
+        from topologicpy.Face import Face
+        from topologicpy.Vector import Vector
         
         def getSubTopologies(topology, subTopologyClass):
             topologies = []
@@ -4053,6 +4195,7 @@ class Topology():
         elif (Topology.Type(topology) > Topology.TypeID("Face")):
             _ = topology.Faces(None, topFaces)
         for aFace in topFaces:
+            f_dir = Face.Normal(aFace)
             ib = []
             _ = aFace.InternalBoundaries(ib)
             if(len(ib) > 0):
@@ -4060,6 +4203,10 @@ class Topology():
                 for aTriFace in triFaces:
                     wire = aTriFace.ExternalBoundary()
                     faceVertices = getSubTopologies(wire, topologic.Vertex)
+                    temp_face = Face.ByWire(wire)
+                    temp_dir = Face.Normal(temp_face)
+                    if Vector.IsAntiParallel(f_dir, temp_dir):
+                        faceVertices.reverse()
                     f = []
                     for aVertex in faceVertices:
                         try:
@@ -4073,6 +4220,10 @@ class Topology():
                 wire =  aFace.ExternalBoundary()
                 #wire = topologic.WireUtility.RemoveCollinearEdges(wire, 0.1) #This is an angle Tolerance
                 faceVertices = getSubTopologies(wire, topologic.Vertex)
+                temp_face = Face.ByWire(wire)
+                temp_dir = Face.Normal(temp_face)
+                if Vector.IsAntiParallel(f_dir, temp_dir):
+                    faceVertices.reverse()
                 f = []
                 for aVertex in faceVertices:
                     try:

topologicpy/Vector.py

@@ -526,7 +526,6 @@ class Vector(list):
         if np.isclose(dot_product, 1.0) or np.isclose(dot_product, -1.0):
             return True
         else:
-            # Compute bisecting vector
             return False
     
     @staticmethod

topologicpy/Wire.py

@@ -485,7 +485,7 @@ class Wire(Topology):
                     if len(st) > 1:
                         e1 = st[0]
                         e2 = st[1]
-                        if not Edge.IsCollinear(e1, e2, angTolerance=angTolerance, tolerance=tolerance):
+                        if not Edge.IsCollinear(e1, e2, tolerance=tolerance):
                             e1 = Edge.Reverse(e1, tolerance=tolerance)
                             bisector = Edge.ByVertices([vertices[i], newVertices[i]], tolerance=tolerance)
                             nv = Edge.VertexByDistance(bisector, distance=finalMiterThreshold, origin=Edge.StartVertex(bisector), tolerance=0.0001)
@@ -510,7 +510,7 @@ class Wire(Topology):
             for v in vertices:
                 edges = Topology.SuperTopologies(v, c, topologyType="edge")
                 if len(edges) == 2:
-                    if not Edge.IsCollinear(edges[0], edges[1], angTolerance=angTolerance, tolerance=tolerance):
+                    if not Edge.IsCollinear(edges[0], edges[1], tolerance=tolerance):
                         adjacentVertices = Topology.AdjacentTopologies(v, c)
                         total = 0
                         for adjV in adjacentVertices:
@@ -2210,7 +2210,7 @@ class Wire(Topology):
                 edges = []
                 _ = aVertex.Edges(wire, edges)
                 if len(edges) > 1:
-                    if not Edge.IsCollinear(edges[0], edges[1], angTolerance=angTolerance, tolerance=tolerance):
+                    if not Edge.IsCollinear(edges[0], edges[1], tolerance=tolerance):
                         wire_verts.append(aVertex)
                 else:
                     wire_verts.append(aVertex)

topologicpy/version.py

@@ -1 +1 @@
-__version__ = '0.7.0'
+__version__ = '0.7.3'

{topologicpy-0.7.0.dist-info → topologicpy-0.7.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: topologicpy
-Version: 0.7.0
+Version: 0.7.3
 Summary: An Advanced Spatial Modelling and Analysis Software Library for Architecture, Engineering, and Construction.
 Author-email: Wassim Jabi <wassim.jabi@gmail.com>
 License: GNU AFFERO GENERAL PUBLIC LICENSE

{topologicpy-0.7.0.dist-info → topologicpy-0.7.3.dist-info}/RECORD

@@ -6,9 +6,9 @@ topologicpy/Color.py,sha256=UlmRcCSOhqcM_OyMWz4t3Kr75KcgXDhz3uctAJ2n7Ic,18031
 topologicpy/Context.py,sha256=ppApYKngZZCQBFWaxIMi2z2dokY23c935IDCBosxDAE,3055
 topologicpy/DGL.py,sha256=RpkLnAzjg6arY7cEMs2pDFYzRdkerVg1Wbm9hcE3QaM,138991
 topologicpy/Dictionary.py,sha256=pMbfE2RYGCNpVr2x58qiHRc-aBWnp1jLlyzwS9nz6-w,25891
-topologicpy/Edge.py,sha256=1Kg1-xkz1_9SVF9M2UuSf04uMuearTfhXMhM9XD3poA,49477
+topologicpy/Edge.py,sha256=jbJHyPHlLF94RkUKL479B4dPC202K0Xd7bQsPcFwoHc,51296
 topologicpy/EnergyModel.py,sha256=UBLim01lZLikVQmJAHEeja-KvF4tgzTxsSxwNDaezz4,52500
-topologicpy/Face.py,sha256=AgCtQ96QRb5Vej1fpiYp4S_GSDtxu50SKV0XgTVgnqQ,96750
+topologicpy/Face.py,sha256=BOU6hN-4uCaXW_lbvuOlntrz_Q0jEiL5tIHAFuqWw_U,97195
 topologicpy/Graph.py,sha256=GsWySXL-cTd7bl-QWpl3Aw9Hr0KyACzIa4iRenTB63s,368030
 topologicpy/Grid.py,sha256=XM0iQtQbMQoYHf7S0ppSe-dxy93Y9VpI_vUkElWqnYA,17050
 topologicpy/Helper.py,sha256=07V9IFu5ilMpvAdZVhIbdBOjBJSRTtJ0BfR1IoRaRXU,17743
@@ -20,14 +20,14 @@ topologicpy/Polyskel.py,sha256=MYHKFOQBlUNqoUhAdOcKRIHpSk0dWWVrZgXK34NkvFM,15936
 topologicpy/Shell.py,sha256=6hidTQ6MW5q_es-WbTeI_yt7Sd7BMxWU_qfoherVZhE,74343
 topologicpy/Speckle.py,sha256=rUS6PCaxIjEF5_fUruxvMH47FMKg-ohcoU0qAUb-yNM,14267
 topologicpy/Sun.py,sha256=3tYb8kssU882lE1gEWg2mxDvCCY_LAVElkyUT6wa-ZU,36935
-topologicpy/Topology.py,sha256=t8Zc25xRNsJWu2e4n_hQ08WjRFBVRm1ujuSueBYauiU,292372
-topologicpy/Vector.py,sha256=4mvvValXDCXFzLupjC_lpsr7LxKVXLmVZUe6H3Fo7f8,29604
+topologicpy/Topology.py,sha256=limcBROvZgD5xanXLi57WgOvR9e93-MkYZcFtFgctu0,298938
+topologicpy/Vector.py,sha256=FHbrCb9GVLOUV_kqcplh4D88CVxlID6qX_wEQOw4rD0,29565
 topologicpy/Vertex.py,sha256=YgbbCcqABvb97z2-yEytpp5T3yoZPlQplR_vMQkQ9OA,65180
-topologicpy/Wire.py,sha256=Q8jpH797oSkjso-ipAxsCs6ta2BxTHrYuBbSLewmV0w,138505
+topologicpy/Wire.py,sha256=MUEboxo11kMgwnZySSkwiyzBG2wv0wPiinf2cW4UVv8,138424
 topologicpy/__init__.py,sha256=D7ky87CAQMiS2KE6YLvcTLkTgA2PY7rASe6Z23pjp9k,872
-topologicpy/version.py,sha256=09hol5ovL6ArILJIJZ-MiscnOjKa6RjaETmtdgJpW2c,22
-topologicpy-0.7.0.dist-info/LICENSE,sha256=TfPDBt3ar0uv_f9cqCDMZ5rIzW3CY8anRRd4PkL6ejs,34522
-topologicpy-0.7.0.dist-info/METADATA,sha256=zHBR-EYUrl66gK-J54RV5KYBg_h2bxXX5tfQMXJ4El0,46950
-topologicpy-0.7.0.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-topologicpy-0.7.0.dist-info/top_level.txt,sha256=J30bDzW92Ob7hw3zA8V34Jlp-vvsfIkGzkr8sqvb4Uw,12
-topologicpy-0.7.0.dist-info/RECORD,,
+topologicpy/version.py,sha256=i5zhgfscuWC8vabp6pP1S6z_O5355x2aIv0so00_cmo,22
+topologicpy-0.7.3.dist-info/LICENSE,sha256=TfPDBt3ar0uv_f9cqCDMZ5rIzW3CY8anRRd4PkL6ejs,34522
+topologicpy-0.7.3.dist-info/METADATA,sha256=TdJMT-64j4LX5i_DSX8qg-5jNgUSETUla-Mi9AMjkic,46950
+topologicpy-0.7.3.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+topologicpy-0.7.3.dist-info/top_level.txt,sha256=J30bDzW92Ob7hw3zA8V34Jlp-vvsfIkGzkr8sqvb4Uw,12
+topologicpy-0.7.3.dist-info/RECORD,,