Python genz_oscillatory_exactの例

コード例 #1

ファイル: test_cubature.py プロジェクト: antonl/cubature

def test_genz_oscillatory_exact_d4():
    exact = np.sin(29988 - 7*np.pi/34) * \
            np.sin(51)*np.sin(2601/2)*np.sin(6273/2)*np.sin(25500) / \
            784125
    calculated = ti.genz_oscillatory_exact(51, np.array([2, 123, 51, 1000]),
            11/17)
    assert np.allclose([exact], [calculated])

コード例 #2

ファイル: test_cubature.py プロジェクト: B-Rich/cubature-1

def test_genz_oscillatory_exact_d4():
    exact = np.sin(29988 - 7*np.pi/34) * \
            np.sin(51)*np.sin(2601/2)*np.sin(6273/2)*np.sin(25500) / \
            784125
    calculated = ti.genz_oscillatory_exact(51, np.array([2, 123, 51, 1000]),
                                           11 / 17)
    assert np.allclose([exact], [calculated])

コード例 #3

ファイル: test_cubature.py プロジェクト: antonl/cubature

def test_hcubature_genz_oscillatory_d2():
    u = 2*np.pi*15/609
    a = np.array([15.51, 2], dtype=float)
    n = 3

    xmin = np.zeros((2,))
    xmax = np.ones((2,)) * n

    exact = ti.genz_oscillatory_exact(n, a, u)

    # check that integrand is callable
    ti.genz_oscillatory(np.ones((2,), dtype=float), a, u)
    
    val, err = cubature(ti.genz_oscillatory, 2, 1, xmin, xmax, args=(a, u),
            adaptive='h')
    assert np.allclose(exact, val) 

コード例 #4

def test_hcubature_genz_oscillatory_d2():
    u = 2*np.pi*15/609
    a = np.array([15.51, 2], dtype=float)
    n = 3

    xmin = np.zeros((2,))
    xmax = np.ones((2,)) * n

    exact = ti.genz_oscillatory_exact(n, a, u)

    # check that integrand is callable
    ti.genz_oscillatory(np.ones((2,), dtype=float), a, u)

    val, err = cubature(ti.genz_oscillatory, 2, 1, xmin, xmax, args=(a, u),
            adaptive='h')
    assert np.allclose(exact, val)

コード例 #5

ファイル: test_cubature.py プロジェクト: antonl/cubature

def test_genz_oscillatory_exact_d3():
    exact = 4 * np.cos(1056 + np.pi/7)*np.sin(12)*np.sin(306)*np.sin(738)/6273
    calculated = ti.genz_oscillatory_exact(12, np.array([2, 123, 51]), 1/14)
    assert np.allclose([exact], [calculated])

コード例 #6

ファイル: test_cubature.py プロジェクト: B-Rich/cubature-1

def test_genz_oscillatory_exact_d3():
    exact = 4 * np.cos(1056 + np.pi /
                       7) * np.sin(12) * np.sin(306) * np.sin(738) / 6273
    calculated = ti.genz_oscillatory_exact(12, np.array([2, 123, 51]), 1 / 14)
    assert np.allclose([exact], [calculated])