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, 1],
        control_points=[
            [0, 0, 0],
            [1, 0, 0],
            [0, 1, 0],
            [1, 1, 0],
            [0, 0, 3],
            [1, 0, 3],
            [0, 1, 3],
            [1, 1, 3],
        ],
    )
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 3.0))
    # Check if equal after refinement
    bezier.elevate_degrees([0, 0, 1, 1, 2, 2])
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 3.0))
    # Test knot insertion
    bspline = bezier.bspline
    bspline.insert_knots(0, c.np.random.rand(10))
    self.assertTrue(c.np.allclose(bezier.integrate.volume(), 3.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(), 3.0))