def METHOD_NAME(self):
    p = np.array([[1, 2], [5, 1], [2, 7], [7, 5]]).T
    sp, sort_ind = sort_points.sort_point_pairs(p)
    # Use numpy arrays to ease comparison of points
    known_lines = np.array([[1, 2], [2, 7], [7, 5], [5, 1]]).T
    known_sort_ind = np.array([0, 2, 3, 1])
    self.assertTrue(np.allclose(known_lines, sp))
    self.assertTrue(np.allclose(known_sort_ind, sort_ind))