def METHOD_NAME(self):
    # Fracture plane defined by x + y + z = 1
    f_1 = pp.PlaneFracture(
        np.array([[0, 1, 1, 0], [0, 0, 1, 1], [1, 0, -1, 0]]), check_convexity=False
    )
    c_known = np.array([0.5, 0.5, 0]).reshape((3, 1))
    self.assertTrue(np.allclose(c_known, f_1.center))