def test_cubature_two_exact(): radius = 1. d = 1 exact = 2 val = ti.cubature_two_exact(d, radius) assert np.allclose([exact], [val]) d = 2 exact = np.pi val = ti.cubature_two_exact(d, radius) assert np.allclose([exact], [val]) d = 3 exact = 4 / 3 * np.pi val = ti.cubature_two_exact(d, radius) assert np.allclose([exact], [val])
def test_cubature_two_exact(): radius = 1. d = 1 exact = 2 val = ti.cubature_two_exact(d, radius) assert np.allclose([exact], [val]) d = 2 exact = np.pi val = ti.cubature_two_exact(d, radius) assert np.allclose([exact], [val]) d = 3 exact = 4/3*np.pi val = ti.cubature_two_exact(d, radius) assert np.allclose([exact], [val])
def test_hcubature_cubature_two(): radius = 0.68244456511919859846 d = 2 xmin = -np.ones((d,)) xmax = np.ones((d,)) exact = ti.cubature_two_exact(d, radius) val, err = cubature(ti.cubature_two, d, 1, xmin, xmax, args=(radius,), abserr=1e-4, relerr=1e-4, maxEval=1000000) true_error = np.abs(val - exact) assert true_error < 1e-4