This is an automated email from the git hooks/post-receive script. ppm-guest pushed a commit to annotated tag v0.37 in repository libmath-prime-util-perl.
commit b61f82e78aef3a6c51e9787eca610d94a3ccc550 Author: Dana Jacobsen <d...@acm.org> Date: Wed Jan 22 23:10:37 2014 -0800 Add some more terms to integer zeta table. Speeds up prime_count_approx by 1.5x --- util.c | 22 +++++++++++++++++++--- 1 file changed, 19 insertions(+), 3 deletions(-) diff --git a/util.c b/util.c index 796a8c6..54466bc 100644 --- a/util.c +++ b/util.c @@ -1460,6 +1460,22 @@ static const long double riemann_zeta_table[] = { 0.0000000000036379795473786511902372363L, 0.0000000000018189896503070659475848321L, 0.0000000000009094947840263889282533118L, + 0.0000000000004547473783042154026799112L, + 0.0000000000002273736845824652515226821L, + 0.0000000000001136868407680227849349105L, + 0.0000000000000568434198762758560927718L, + 0.0000000000000284217097688930185545507L, + 0.0000000000000142108548280316067698343L, + 0.00000000000000710542739521085271287735L, + 0.00000000000000355271369133711367329847L, + 0.00000000000000177635684357912032747335L, + 0.000000000000000888178421093081590309609L, + 0.000000000000000444089210314381336419777L, + 0.000000000000000222044605079804198399932L, + 0.000000000000000111022302514106613372055L, + 0.0000000000000000555111512484548124372374L, + 0.0000000000000000277555756213612417258163L, + 0.0000000000000000138777878097252327628391L, }; #define NPRECALC_ZETA (sizeof(riemann_zeta_table)/sizeof(riemann_zeta_table[0])) @@ -1469,9 +1485,9 @@ static const long double riemann_zeta_table[] = { * The Cephes zeta function uses a series (2k)!/B_2k which converges rapidly * and has a very wide range of values. We use it here for some values. * - * Note: Calculations here are done on long doubles and we try to generate ~17 - * digits of accuracy. When these are returned to Perl they get put in - * a standard 64-bit double, so don't expect more than 15 digits. + * Note: Calculations here are done on long doubles and we try to generate as + * much accuracy as possible. They will get returned to Perl as an NV, + * which is typically a 64-bit double with 15 digits. * * For values 0.5 to 5, this code uses the rational Chebyshev approximation * from Cody and Thacher. This method is extraordinarily fast and very -- Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/pkg-perl/packages/libmath-prime-util-perl.git _______________________________________________ Pkg-perl-cvs-commits mailing list Pkg-perl-cvs-commits@lists.alioth.debian.org http://lists.alioth.debian.org/cgi-bin/mailman/listinfo/pkg-perl-cvs-commits