def METHOD_NAME(self):
    """
    When more than 1000 points passed they are sub-sampled in order to
    compute a concave hull
    """
    iris = Table.from_file("iris")
    data = np.repeat(iris.X[:, 2:4], 10, axis=0)  # more than 1000 points
    clusters = np.array([0] * 50 * 10 + [1] * 50 * 10 + [2] * 50 * 10)
    hulls = compute_hulls(data, clusters)
    self.assertEqual(3, len(hulls))
    self.assertEqual(2, hulls[0].shape[1])  # hull have x and y
    self.assertEqual(2, hulls[1].shape[1])  # hull have x and y
    self.assertEqual(2, hulls[2].shape[1])  # hull have x and y