cutThruAll fix for non-planar faces (#604)

* reworked cutThruAll to support non-planar faces

* Added test

Author: adam-urbanczyk

cadquery/cq.py

@@ -3084,20 +3084,9 @@ def cutThruAll(self, clean: bool = True, taper: float = 0) -> "Workplane":
         self.ctx.pendingWires = []
 
         solidRef = self.findSolid()
-        faceRef = self.findFace()
-
-        # if no faces on the stack take the nearest face parallel to the plane zDir
-        if not faceRef:
-            # first select all with faces with good orietation
-            sel1 = PerpendicularDirSelector(self.plane.zDir)
-            faces = sel1.filter(solidRef.Faces())
-            # then select the closest
-            sel2 = NearestToPointSelector(self.plane.origin.toTuple())
-            faceRef = sel2.filter(faces)[0]
-
         rv = []
         for solid in solidRef.Solids():
-            s = solid.dprism(faceRef, wires, thruAll=True, additive=False, taper=-taper)
+            s = solid.dprism(None, wires, thruAll=True, additive=False, taper=-taper)
 
             if clean:
                 s = s.clean()

cadquery/occ_impl/shapes.py

@@ -2637,7 +2637,7 @@ def sweep_multi(
 
     def dprism(
         self,
-        basis: Face,
+        basis: Optional[Face],
         profiles: List[Wire],
         depth: Optional[float] = None,
         taper: float = 0,
@@ -2659,7 +2659,12 @@ def dprism(
         for p in sorted_profiles:
             face = Face.makeFromWires(p[0], p[1:])
             feat = BRepFeat_MakeDPrism(
-                shape, face.wrapped, basis.wrapped, taper * DEG2RAD, additive, False
+                shape,
+                face.wrapped,
+                basis.wrapped if basis else TopoDS_Face(),
+                taper * DEG2RAD,
+                additive,
+                False,
             )
 
             if thruAll or depth is None:

tests/test_cadquery.py

@@ -1409,6 +1409,28 @@ def testCutThroughAll(self):
         self.saveModel(t)
         self.assertEqual(13, t.faces().size())
 
+        # no planar faces
+        sphere_r = 10.0
+
+        r = (
+            Workplane()
+            .sphere(sphere_r)
+            .workplane()
+            .circle(sphere_r / 2.0)
+            .cutThruAll()
+            .workplane()
+            .transformed(rotate=(90, 0, 0))
+            .circle(sphere_r / 2.0)
+            .cutThruAll()
+            .workplane()
+            .transformed(rotate=(0, 90, 0))
+            .circle(sphere_r / 2.0)
+            .cutThruAll()
+        )
+
+        self.assertTrue(r.val().isValid())
+        self.assertEqual(r.faces().size(), 7)
+
     def testCutToFaceOffsetNOTIMPLEMENTEDYET(self):
         """
             Tests cutting up to a given face, or an offset from a face