WIP: Remove fracture tips as mesh size points in 2d domains

Implementation must be cleaned up

Author: keileg

src/porepy/fracs/fracture_network_2d.py

@@ -564,9 +564,11 @@ def _set_mesh_size_fields(
             extra_points = (
                 np.array([d[0] for d in info]).T if len(info) > 0 else np.empty((3, 0))
             )
+            if extra_points.size == 0:
+                continue
 
             points, _, ind_map = pp.array_operations.uniquify_point_set(
-                np.hstack((end_points, extra_points)), tol=tol
+                extra_points, tol=tol
             )
             # Distance to other objects for each point, as computed previously. Assign
             # h_frac or h_bound to the endpoints, depending on whether the line is a
@@ -578,7 +580,7 @@ def _set_mesh_size_fields(
 
             other_object_distances_all = np.hstack(
                 (
-                    np.array([h_end, h_end]),
+                    # np.array([h_end, h_end]),
                     np.array(
                         [
                             d[1] if d[1] > 0 else mesh_size_computer.h_frac()
@@ -595,7 +597,7 @@ def _set_mesh_size_fields(
                 min_dist = np.min(other_object_distances_all[inds])
                 other_object_distances.append(min_dist)
 
-            if points.shape[1] > 0:
+            if points.size > 0:
                 # If there is more than one point in addition to the end points, we can
                 # compute the point-point distances in pairs along this line.
                 point_point_distances = pp.distances.pointset(points, max_diag=True)