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)