Example #1
0
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])
Example #2
0
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])
Example #3
0
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
Example #4
0
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