Checks to actually follow up on:
[DONE] check/check_dirichlet_eta.cc -- s < 0. Dive into polylog. Could use reflection.
[DONE] check/check_jacobi_cn.cc -- Completely ruined except for |k| == 1. Use Boost (WTF doesn't do k<0).
[DONE] check/check_jacobi_dn.cc -- Completely ruined except for |k| == 1. Use Boost (WTF doesn't do k<0).
[DONE] check/check_jacobi_sn.cc -- Completely ruined except for |k| == 1. Use Boost (WTF doesn't do k<0).
[DONE] check/check_psi.cc -- x < -2.0
[] check/check_sph_bessel.cc -- n == 1, large distances, large orders.
[DONE] check/check_lfalling_factorial -- Logs of negative numbers. Need decision. lgamma_t?
[] check/check_cylhankel[12] -- negative fractional order
Checks to actually follow up on: [DONE] check/check_dirichlet_eta.cc -- s < 0. Dive into polylog. Could use reflection. [DONE] check/check_jacobi_cn.cc -- Completely ruined except for |k| == 1. Use Boost (WTF doesn't do k<0). [DONE] check/check_jacobi_dn.cc -- Completely ruined except for |k| == 1. Use Boost (WTF doesn't do k<0). [DONE] check/check_jacobi_sn.cc -- Completely ruined except for |k| == 1. Use Boost (WTF doesn't do k<0). [DONE] check/check_psi.cc -- x < -2.0 [] check/check_sph_bessel.cc -- n == 1, large distances, large orders. [DONE] check/check_lfalling_factorial -- Logs of negative numbers. Need decision. lgamma_t? [] check/check_cylhankel[12] -- negative fractional order