def METHOD_NAME(tree: BoundingBoxTree, midpoint_tree: BoundingBoxTree, mesh: Mesh,
                           points: npt.NDArray[np.floating]) -> npt.NDArray[np.int32]:
    """Compute closest mesh entity to a point.
    Args:
        tree: bounding box tree for the entities.
        midpoint_tree: A bounding box tree with the midpoints of all
            the mesh entities. This is used to accelerate the search.
        mesh: The mesh.
        points: The points to check for collision, ``shape=(num_points,3)``.
    Returns:
        Mesh entity index for each point in ``points``. Returns -1 for a
        point if the bounding box tree is empty.
    """
    return _cpp.geometry.METHOD_NAME(tree._cpp_object, midpoint_tree._cpp_object, mesh._cpp_object, points)