def test_hcubature_cubature_one():
d = 3
xmax = np.ones((d,), dtype=float)
xmin = np.zeros_like(xmax)
ti.cubature_one(xmin)
exact = ti.cubature_one_exact(xmax)
val, err = cubature(ti.cubature_one, d, 1, xmin, xmax)
assert np.allclose([exact], [val])
def test_hcubature_cubature_one():
d = 3
xmax = np.ones((d, ), dtype=float)
xmin = np.zeros_like(xmax)
ti.cubature_one(xmin)
exact = ti.cubature_one_exact(xmax)
val, err = cubature(ti.cubature_one, d, 1, xmin, xmax)
assert np.allclose([exact], [val])
def test_cubature_zero():
d = 5
k = 2/np.sqrt(np.pi)
x = np.ones((d,), dtype=float)
expected = k**d
assert np.allclose(ti.cubature_one(x), expected)
def test_cubature_zero():
d = 5
k = 2 / np.sqrt(np.pi)
x = np.ones((d, ), dtype=float)
expected = k**d
assert np.allclose(ti.cubature_one(x), expected)
