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

update conrey profile

Browse source
This commit is contained in:
Pascal 2016-09-05 16:24:44 +02:00
parent 403c2efb6d
commit d73fdb7f10
1 changed files with 132 additions and 34 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

166
acb_dirichlet/profile/p-conrey.c
View file

@ -23,32 +23,34 @@
******************************************************************************/
#include <string.h>
#include "acb_dirichlet.h"
#include "profiler.h"
int main()
#define LOG 0
#define CSV 1
#define JSON 2
typedef ulong (*do_f) (ulong q1, ulong q2);
static ulong
do_gcd(ulong q1, ulong q2)
{
slong iter, k, n, nref;
ulong q, maxq;
ulong n, q, k;
maxq = 4000;
flint_printf("loop over all (Z/q)* for q<=%wu....\n", maxq);
fflush(stdout);
flint_printf("gcd................ ");
TIMEIT_ONCE_START
for (n = 0, q = 2; q <= maxq; q++)
for (n = 0, q = q1; q <= q2; q++)
for (k = 1; k < q; k++)
if (n_gcd(k, q) > 1)
continue;
else
n++;
TIMEIT_ONCE_STOP
nref = n;
n += (n_gcd(k, q) == 1);
flint_printf("conrey.... ");
TIMEIT_ONCE_START
for (n = 0, q = 2; q <= maxq; q++)
return n;
}
static ulong
do_conrey(ulong q1, ulong q2)
{
ulong n, q;
for (n = 0, q = q1; q <= q2; q++)
{
acb_dirichlet_group_t G;
acb_dirichlet_conrey_t x;
@ -63,16 +65,16 @@ int main()
acb_dirichlet_conrey_clear(x);
acb_dirichlet_group_clear(G);
}
TIMEIT_ONCE_STOP
if (n != nref)
{
flint_printf("FAIL: wrong number of elements %wu != %wu\n\n",n, nref);
abort();
}
flint_printf("chars.... ");
TIMEIT_ONCE_START
for (n = 0, q = 2; q <= maxq; q++)
return n;
}
static ulong
do_chars(ulong q1, ulong q2)
{
ulong n, q;
for (n = 0, q = q1; q <= q2; q++)
{
acb_dirichlet_group_t G;
acb_dirichlet_char_t chi;
@ -82,20 +84,116 @@ int main()
acb_dirichlet_char_one(chi, G);
n++;
for (; acb_dirichlet_char_next(chi, G) >= 0; n++);
acb_dirichlet_char_clear(chi);
acb_dirichlet_group_clear(G);
}
TIMEIT_ONCE_STOP
if (n != nref)
return n;
}
static ulong
do_gcdpluscond(ulong q1, ulong q2)
{
ulong n, q, k;
for (n = 0, q = q1; q <= q2; q++)
{
flint_printf("FAIL: wrong number of elements %wu != %wu\n\n",n, nref);
abort();
acb_dirichlet_group_t G;
acb_dirichlet_group_init(G, q);
for (k = 1; k < q; k++)
{
/* known factors -> faster gcd */
slong i;
if (G->q_even > 1 && k % 2 == 0)
continue;
for (i = G->neven; i < G->num; i++)
if (k % G->P[i].p == 0)
break;
if (i == G->num)
acb_dirichlet_ui_conductor(G, k);
}
n += G->phi_q;
acb_dirichlet_group_clear(G);
}
return n;
}
int main(int argc, char *argv[])
{
int out;
ulong n, nref, maxq = 5000;
int l, nf = 4;
do_f func[4] = { do_gcd, do_conrey, do_chars, do_gcdpluscond };
char * name[4] = { "gcd", "conrey", "chars", "gcd+cond" };
int i, ni = 5;
ulong qmin[5] = { 2, 1000, 10000, 100000, 1000000 };
ulong qmax[5] = { 500, 3000, 11000, 100100, 1000010 };
if (argc < 2)
out = LOG;
else if (!strcmp(argv[1], "json"))
out = JSON;
else if (!strcmp(argv[1], "csv"))
out = CSV;
else if (!strcmp(argv[1], "log"))
out = LOG;
else
{
printf("usage: %s [log|csv|json]\n", argv[0]);
abort();
}
if (out == CSV)
flint_printf("# %-6s, %7s, %7s, %7s\n","name", "qmin", "qmax", "time");
for (i = 0; i < ni; i++)
{
if (out == LOG)
{
flint_printf("loop over all (Z/q)* for %wu<=q<=%wu\n", qmin[i], qmax[i]);
}
for (l = 0; l < nf; l++)
{
if (out == LOG)
flint_printf("%-8s ... ",name[l]);
else if (out == CSV)
flint_printf("%-8s, %7d, %7d, ",name[l],qmin[i],qmax[i]);
else if (out == JSON)
flint_printf("{ \"name\": \"%s\", \"qmin\": %d, \"qmax\": %d, \"time\": ",
name[l],qmin[i],qmax[i]);
TIMEIT_ONCE_START
(func[l])(qmin[i], qmax[i]);
TIMEIT_ONCE_STOP
if (l == 0)
nref = n;
else if (n != nref)
{
flint_printf("FAIL: wrong number of elements %wu != %wu\n\n",n, nref);
abort();
}
if (out == JSON)
flint_printf("}\n");
else
flint_printf("\n");
}
}
flint_cleanup();
flint_printf("PASS\n");
return EXIT_SUCCESS;
}