1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
|
/* This file is public domain. Author: Pascal Molin. */
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <flint/profiler.h>
#include <flint/acb.h>
#include <flint/acb_dirichlet.h>
static int usage(char *argv[])
{
printf("Computes central values (s = 0.5) of Dirichlet L-functions.\n\n");
printf("usage: %s [--quiet] [--check] [--prec <bits>] qmin qmax\n", argv[0]);
return 1;
}
int main(int argc, char *argv[])
{
int i, check = 0, out = 1;
slong prec = 100, digits = 30;
ulong qmin, qmax, q;
acb_t s;
acb_dirichlet_hurwitz_precomp_t pre;
if (argc < 3)
return usage(argv);
for (i = 1; i < argc - 2; i++)
{
if (!strcmp(argv[i], "--quiet"))
out = 0;
else if (!strcmp(argv[i], "--check"))
check = 1;
else if (!strcmp(argv[i], "--prec"))
{
i++;
prec = atol(argv[i]);
digits = floor(prec * 0.3);
}
}
if (argc < i + 2)
return usage(argv);
qmin = atol(argv[i]);
qmax = atol(argv[i+1]);
fflush(stdout);
acb_init(s);
TIMEIT_ONCE_START;
acb_one(s);
acb_div_si(s, s, 2, prec);
acb_dirichlet_hurwitz_precomp_init_num(pre, s, 0,
(qmax - qmin + 1) * 0.5 * qmax, prec);
for (q = qmin; q <= qmax; q++)
{
ulong k;
dirichlet_group_t G;
dirichlet_char_t x;
acb_ptr z;
if (q % 4 == 2)
continue;
dirichlet_group_init(G, q);
dirichlet_char_init(x, G);
z = _acb_vec_init(G->phi_q);
acb_dirichlet_l_vec_hurwitz(z, s, pre, G, prec);
if (out || check)
{
k = 0;
dirichlet_char_one(x, G);
while (dirichlet_char_next(x, G) >= 0)
{
k++;
if (dirichlet_conductor_char(G,x) < q)
continue;
if (acb_contains_zero(z + k))
{
flint_printf("Value could be zero!\n");
flint_printf("%wu,%wu: ", q, x->n);
acb_printn(z + k, digits, 0);
flint_printf("\n");
flint_abort();
}
if (out)
{
flint_printf("%wu,%wu: ", q, x->n);
acb_printn(z + k, digits, 0);
flint_printf("\n");
}
}
}
_acb_vec_clear(z, G->phi_q);
dirichlet_char_clear(x);
dirichlet_group_clear(G);
}
acb_dirichlet_hurwitz_precomp_clear(pre);
acb_clear(s);
TIMEIT_ONCE_STOP;
print_memory_usage();
flint_cleanup();
return 0;
}
|
