[SciPy-Dev] [SciPy-User] ANN: scipy 0.8.0 release candidate 3

Ralf Gommers ralf.gommers@googlemail....
Sun Jul 18 09:22:43 CDT 2010


On Sun, Jul 18, 2010 at 6:30 AM, Derek Homeier <
derek@astro.physik.uni-goettingen.de> wrote:

> >
> > These should have been fixed by r6520:
> > -        self.check_cephes_vs_amos(iv, iv, rtol=1e-12, atol=1e-305)
> > +        self.check_cephes_vs_amos(iv, iv, rtol=5e-9, atol=1e-305)
> >
> >          data(gammaincinv, 'gamma_inv_big_data_ipp-
> > gamma_inv_big_data',
> > -             (0,1), 2, rtol=5e-12),
> > +             (0,1), 2, rtol=1e-11),
> >
> > Can you check if you have these changes, and with what accuracy the
> > tests pass?
> >
>
> I have (rc3 tarball), but it's a different couple of tests.
> 'gamma_inv_big_data_ipp-gamma_inv_big_data' actually still passes with
> rtol=2.1e-12.
>
> These are the minimal accuracy changes needed:
>
> --- /sw/lib/python2.6/site-packages/scipy/special/tests/test_data.py
> 2010-07-11 17:25:24.000000000 +0200
> +++ test_data.py        2010-07-17 21:13:07.000000000 +0200
> @@ -85,7 +85,7 @@
>          data(gamma, 'test_gamma_data_ipp-near_1', 0, 1),
>          data(gamma, 'test_gamma_data_ipp-near_2', 0, 1),
>          data(gamma, 'test_gamma_data_ipp-near_m10', 0, 1),
> -        data(gamma, 'test_gamma_data_ipp-near_m55', 0, 1),
> +        data(gamma, 'test_gamma_data_ipp-near_m55', 0, 1, rtol=7e-12),
>          data(gamma, 'test_gamma_data_ipp-near_0', 0j, 1, rtol=2e-9),
>          data(gamma, 'test_gamma_data_ipp-near_1', 0j, 1, rtol=2e-9),
>          data(gamma, 'test_gamma_data_ipp-near_2', 0j, 1, rtol=2e-9),
> --- /sw/lib/python2.6/site-packages/scipy/special/tests/test_basic.py
> 2010-07-11 17:25:24.000000000 +0200
> +++ test_basic.py       2010-07-17 21:28:07.000000000 +0200
> @@ -1618,7 +1618,7 @@
>
>          # Most error apparently comes from AMOS and not our
> implementation;
>          # there are some problems near integer orders there
> -        assert dc[k] < 1e-9, (v[k], x[k], iv(v[k], x[k]), iv(v[k],
> x[k]+0j))
> +        assert dc[k] < 1.9e-7, (v[k], x[k], iv(v[k], x[k]), iv(v[k],
> x[k]+0j))
>
>      def test_kv_cephes_vs_amos(self):
>          #self.check_cephes_vs_amos(kv, kn, rtol=1e-9, atol=1e-305)
>
> Of course nearly 2e-7 seems a quite high tolerance - don't know if
> that's acceptable.
>
> That does seem a little high. I guess Pauli has to decide whether this
qualifies as a bug or not.

Cheers,
Ralf
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-dev/attachments/20100718/f790bef7/attachment.html 


More information about the SciPy-Dev mailing list