test_precompute_expn_asy.py 583 B

123456789101112131415161718192021222324
  1. from numpy.testing import assert_equal
  2. from scipy.special._testutils import check_version, MissingModule
  3. from scipy.special._precompute.expn_asy import generate_A
  4. try:
  5. import sympy
  6. from sympy import Poly
  7. except ImportError:
  8. sympy = MissingModule("sympy")
  9. @check_version(sympy, "1.0")
  10. def test_generate_A():
  11. # Data from DLMF 8.20.5
  12. x = sympy.symbols('x')
  13. Astd = [Poly(1, x),
  14. Poly(1, x),
  15. Poly(1 - 2*x),
  16. Poly(1 - 8*x + 6*x**2)]
  17. Ares = generate_A(len(Astd))
  18. for p, q in zip(Astd, Ares):
  19. assert_equal(p, q)