def METHOD_NAME(self):
    """
    Test volume integration for splines using numerical integration of the
    Jacobi-Determinant
    """
    # Test 1D
    bezier = c.splinepy.Bezier(
        degrees=[1, 1], control_points=[[0, 0], [2, 0], [0, 1], [2, 1]]
    )
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 2.0))
    # test for other types same spline
    self.assertTrue(c.np.allclose(bezier.bspline.integrate.volume(), 2.0))
    self.assertTrue(
        c.np.allclose(bezier.rationalbezier.integrate.volume(), 2.0)
    )
    self.assertTrue(c.np.allclose(bezier.nurbs.integrate.volume(), 2.0))
    # Check if equal after refinement
    bezier.elevate_degrees([0, 0, 1])
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 2.0))
    # Test knot insertion
    bspline = bezier.bspline
    bspline.insert_knots(0, c.np.random.rand(10))
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 2.0))
    # Move control points along axis does not change volume
    mi = bspline.multi_index
    bspline.cps[mi[3, :]][:, 1] += 10
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 2.0))