def test_hcubature_cubature_zero(): xmin = np.zeros((4,)) xmax = np.array([12, 4, 0.25, 1], dtype=float) # check that it is possible to run ti.cubature_zero(xmax) exact = ti.cubature_zero_exact(xmax) val, err = cubature(ti.cubature_zero, xmax.shape[0], 1, xmin, xmax) assert np.allclose([exact], [val])
def test_hcubature_cubature_zero(): xmin = np.zeros((4, )) xmax = np.array([12, 4, 0.25, 1], dtype=float) # check that it is possible to run ti.cubature_zero(xmax) exact = ti.cubature_zero_exact(xmax) val, err = cubature(ti.cubature_zero, xmax.shape[0], 1, xmin, xmax) assert np.allclose([exact], [val])
def test_cubature_zero_exact(): d = 1 xmin = np.zeros((d,), dtype=float) xmax = np.ones((d,), dtype=float)*np.pi/2 assert np.allclose(0., ti.cubature_zero_exact(xmin)) assert np.allclose(1., ti.cubature_zero_exact(xmax)) d = 3 xmin = np.zeros((d,), dtype=float) assert np.allclose(0., ti.cubature_zero_exact(xmin)) assert np.allclose(1., ti.cubature_zero_exact(xmax)) d = 4 xmin = np.zeros((d,), dtype=float) assert np.allclose(0., ti.cubature_zero_exact(xmin)) assert np.allclose(1., ti.cubature_zero_exact(xmax))
def test_cubature_zero_exact(): d = 1 xmin = np.zeros((d, ), dtype=float) xmax = np.ones((d, ), dtype=float) * np.pi / 2 assert np.allclose(0., ti.cubature_zero_exact(xmin)) assert np.allclose(1., ti.cubature_zero_exact(xmax)) d = 3 xmin = np.zeros((d, ), dtype=float) assert np.allclose(0., ti.cubature_zero_exact(xmin)) assert np.allclose(1., ti.cubature_zero_exact(xmax)) d = 4 xmin = np.zeros((d, ), dtype=float) assert np.allclose(0., ti.cubature_zero_exact(xmin)) assert np.allclose(1., ti.cubature_zero_exact(xmax))