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

length estimate fot odd theta series

Browse source
This commit is contained in:
Pascal 2016-06-15 15:24:49 +02:00
parent c37c09522f
commit ed86f818c8
3 changed files with 8 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
acb_dirichlet/char_next.c
View file

@ -42,6 +42,8 @@ acb_dirichlet_char_next(acb_dirichlet_char_t chi, const acb_dirichlet_group_t G)
chi->expo[k] = 0;
}
/* should do it with log, but not kept yet */
chi->conductor = acb_dirichlet_conductor_ui(G, chi->n);
acb_dirichlet_char_normalize(chi, G);
/* return last index modified */

1
acb_dirichlet/parity_ui.c
View file

@ -36,6 +36,7 @@ acb_dirichlet_parity_ui(const acb_dirichlet_group_t G, ulong a)
if (G->neven && a % 4 == 3)
par++;
/* could replace by jacobi(a, core(q_odd)) */
for (k = G->neven; k < G->num; k++)
{
if (n_jacobi_unsigned(a, G->primes[k]) == -1)

7
acb_dirichlet/theta_length.c
View file

@ -31,8 +31,11 @@
ulong
acb_dirichlet_theta_length_d(ulong q, double x, slong prec)
{
double a = PI / (double)q * x * x;
return ceil(sqrt(((double)prec * LOG2 - log(2 * a)) / a));
double a, la;
a = PI / (double)q * x * x;
la = (a>.3) ? -log(2*a*(1-a)) : .8;
la = ((double)prec * LOG2 + la) / a;
return ceil(sqrt(la)+.5);
}
ulong