I recommended Juan to post to sage-devel because there are people here
who are knowledgeable about zeta and L-function computations. The code
and results look good to me, and I would like to add it to mpmath, but
maybe some people here have comments and can provide a more in-depth
review of the algorithm. It would also be interesting to know how this
implementation compares to zeta implementations already present in
Sage (Pari, lcalc?), with which I'm not too familiar.
Most time is spent in zetasum which just computes \sum_{k=1}^n k^{-s}.
It ought to be next to trivial to add an MPFR-based version of this
function to sage.libs.mpmath, which should give a large speedup. Later
on it could be sped up further by using fastfunlib (when it reaches a
working state).
On Thu, Dec 3, 2009 at 7:55 PM, j. arias-de-reyna <ar...@us.es> wrote:
> [...]
>
> I pretended to send with this message my program and a pdf file with
> some benchmarks. But I do not know how to do it. I will send it to you
> if you are interested.
Juan, if nothing else works, please send the updated files to me and I
will upload them somewhere.
Fredrik
--
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to
sage-devel-unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
