scipy.special.log1p(inf) returns inf.


On 24 January 2014 13:05, Evgeni Burovski <notifications@github.com> wrote:

> In scipy/stats/_continuous_distns.py:
>
> >
> >      def _logpdf(self, x, c):
> > -        return (-1.0-1.0/c) * np.log1p(c*x)
> > +        return -(c + 1.) * self._log1pcx(x, c)
>
> @pbrod <https://github.com/pbrod> what does scipy.special.log1p(inf)return on your system? --- apparently, there's one from cephes/unity.c
> (I cannot test it myself, since I'm on ubuntu precise at the moment)
>
> —
> Reply to this email directly or view it on GitHub<https://github.com/scipy/scipy/pull/3225/files#r9146444>
> .
>