Re: [SciPy-Dev] [SciPy-User] polylogarithm?
![](https://secure.gravatar.com/avatar/a2cef18ab0a86a4e92d82e6bfea198eb.jpg?s=120&d=mm&r=g)
For the record, I created a first pull request, with limited scope, in order to make the C call to cephes described below possible. Further details in https://github.com/scipy/scipy/pull/84 I will switch to the dev list from now on. best, Johann On 09/23/2011 06:18 PM, Pauli Virtanen wrote:
zetac(z) knows how to eat negative arguments, but not zeta(z,1).... Is there a reason why special.zeta does not default to 1+special.zetac for s=1? This would make the behavior of the 2 functions more identical. Probably no reason, except that it wasn't implemented. mpmath is impressive, and in several ways ahead of scipy.special --- or at least in
Fri, 23 Sep 2011 15:21:12 +0200, Johann Cohen-Tanugi wrote: [clip] the parts where the problems overlap, as you can do tricks with arbitrary precision that are not really feasible.
Note that if you need to call the zeta function from the Cython extension, call the C library directly:
cdef extern from "cephes.h": double zeta(double x, double q) double zetac(double x)
and link the extension with the "sc_cephes" library.
***
There's a formula (look in the mpmath sources ;) for the transform from x< 0 to x> 0 for zeta(x, a) for general a. But that needs polylog.
The implementation for zetac(x) for x< 0 seems also a bit incomplete, as it goes only down to -30.8148. It seems this is due to a silly reason, it should use gammaln instead of Gamma to avoid the overflow.
Pauli
_______________________________________________ SciPy-User mailing list SciPy-User@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-user
participants (1)
-
Johann Cohen-Tanugi