Fri Sep 16 12:24:41 CDT 2011
On Fri, Sep 16, 2011 at 7:11 PM, Ralf Gommers
> On Thu, Sep 15, 2011 at 8:09 AM, Johann Cohen-Tanugi
> <email@example.com> wrote:
>> hi there, any chance for a polylog implementation in scipy.special? I
>> know it is there in mpmath, but I thought I would ask anyway.
> If someone (you?) contributes a patch, that would be a great addition to
> scipy.special imho. mpmath is nice, but it doesn't understand ndarrays and
> is way too slow when you want to use the polylog for something like fitting
> Bose-Einstein or Fermi-Dirac distributions.
> It looks like the implementation in mpmath is quite clean and could provide
> a starting point for a Cython/C version.
Maybe it's still too slow, but mpmath.fp.polylog is about 100 times
faster than the multiprecision version (and usually gives nearly full
double-precision accuracy anyway).
It should indeed be quite simple to translate to Cython/C.
More information about the SciPy-User