def test_scipy_poisson_limit():
'''Test that the lower-level routine gives the snae number.
Test numbers are from table1 1, 3 in
Kraft, Burrows and Nousek in
`ApJ 374, 344 (1991) <https://ui.adsabs.harvard.edu/abs/1991ApJ...374..344K>`_
'''
assert_allclose(funcs._scipy_kraft_burrows_nousek(5, 2.5, .99),
(0, 10.67), rtol=1e-3)
assert_allclose(funcs._scipy_kraft_burrows_nousek(np.int32(5.), 2.5, .99),
(0, 10.67), rtol=1e-3)
assert_allclose(funcs._scipy_kraft_burrows_nousek(np.int64(5.), 2.5, .99),
(0, 10.67), rtol=1e-3)
assert_allclose(funcs._scipy_kraft_burrows_nousek(5, np.float32(2.5), .99),
(0, 10.67), rtol=1e-3)
assert_allclose(funcs._scipy_kraft_burrows_nousek(5, np.float64(2.5), .99),
(0, 10.67), rtol=1e-3)
assert_allclose(funcs._scipy_kraft_burrows_nousek(5, 2.5, np.float32(.99)),
(0, 10.67), rtol=1e-3)
assert_allclose(funcs._scipy_kraft_burrows_nousek(5, 2.5, np.float64(.99)),
(0, 10.67), rtol=1e-3)
conf = funcs.poisson_conf_interval([5, 6], 'kraft-burrows-nousek',
background=[2.5, 2.],
confidence_level=[.99, .9])
assert_allclose(conf[:, 0], (0, 10.67), rtol=1e-3)
assert_allclose(conf[:, 1], (0.81, 8.99), rtol=5e-3)
def test_scipy_poisson_limit():
'''Test that the lower-level routine gives the snae number.
Test numbers are from table1 1, 3 in
Kraft, Burrows and Nousek in
`ApJ 374, 344 (1991) <http://adsabs.harvard.edu/abs/1991ApJ...374..344K>`_
'''
assert_allclose(funcs._scipy_kraft_burrows_nousek(5., 2.5, .99),
(0, 10.67), rtol=1e-3)
conf = funcs.poisson_conf_interval([5., 6.], 'kraft-burrows-nousek',
background=[2.5, 2.],
conflevel=[.99, .9])
assert_allclose(conf[:, 0], (0, 10.67), rtol=1e-3)
assert_allclose(conf[:, 1], (0.81, 8.99), rtol=5e-3)