Tuning of Fracture projection to 2d plane tolerance.

src/porepy/fracs/fractures.py

@@ -397,7 +397,7 @@ def intersects(self, other, tol, check_point_contact=True):
         if int_points.shape[1] > 1:
             int_points, _, _ \
                 = setmembership.unique_columns_tol(int_points, tol=tol)
-                
+
         # There should be at most two of these points.
         # In some cases, likely involving extrusion, several segments may lay
         # essentially in the fracture plane, producing more than two segments.
@@ -1613,8 +1613,8 @@ def _points_2_plane(self, p_loc, edges_loc, p_ind_loc):
         p_2d = rot.dot(p_loc)
 
         extent = p_2d.max(axis=1) - p_2d.min(axis=1)
-
-        assert extent[2] < np.max(extent[:2]) * self.tol * 5
+        lateral_extent = np.maximum(np.max(extent[2]), 1)
+        assert extent[2] < lateral_extent * self.tol * 10
         # Dump third coordinate
         p_2d = p_2d[:2]