hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-06 01:41:39 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

optimize Euler product

Browse source
This commit is contained in:
Fredrik Johansson 2013-03-15 15:07:07 +01:00
parent 8b3689d43f
commit f9612fe95d
1 changed files with 6 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
zeta/zeta_ui_euler_product.c
View file

@ -61,7 +61,7 @@ fmprb_zeta_inv_ui_euler_product(fmprb_t z, ulong s, long prec)
{
long wp, powprec;
double powmag;
fmprb_t t;
fmprb_t t, u;
ulong M;
mp_limb_t p;
@ -75,6 +75,7 @@ fmprb_zeta_inv_ui_euler_product(fmprb_t z, ulong s, long prec)
wp = prec + FLINT_BIT_COUNT(prec) + (prec/s) + 4;
fmprb_init(t);
fmprb_init(u);
/* z = 1 - 2^(-s) */
fmprb_set_ui(z, 1UL);
@ -95,8 +96,10 @@ fmprb_zeta_inv_ui_euler_product(fmprb_t z, ulong s, long prec)
break;
M = p;
fmprb_ui_pow_ui(t, p, s, powprec);
fmprb_div(t, z, t, powprec);
fmprb_set_round(u, z, powprec);
fmprb_div(t, u, t, powprec);
fmprb_sub(z, z, t, wp);
p = n_nextprime(p, 0);
@ -104,6 +107,7 @@ fmprb_zeta_inv_ui_euler_product(fmprb_t z, ulong s, long prec)
add_error(z, M, s);
fmprb_clear(t);
fmprb_clear(u);
}
void