Beispiel #1
0
 def test(N, rtol=1e-15, atol=1e-14):
     x, w = orth._roots_hermite_asy(N)
     H = hermite_recursion(N+1, x)
     assert_allclose(H[-1,:], np.zeros(N), rtol, atol)
     assert_allclose(sum(w), sqrt(np.pi), rtol, atol)
 def test(N, rtol=1e-15, atol=1e-14):
     x, w = orth._roots_hermite_asy(N)
     H = hermite_recursion(N+1, x)
     assert_allclose(H[-1,:], np.zeros(N), rtol, atol)
     assert_allclose(sum(w), sqrt(np.pi), rtol, atol)