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
arb/acb_dirichlet/pairing_conrey.c

30 lines
861 B
C
Raw Normal View History

update licence
2016-09-06 15:32:23 +02:00
/*
huge refactor - split files, one for each function - use acb_dirichlet prefix everywhere - distinguish pairing / character value + fix bug on exponent
2016-03-22 11:29:21 +01:00
Copyright (C) 2016 Pascal Molin
[dirichlet] add conrey type to handle logs + char type - try to handle even and odd components the same way in the dirichlet group - switch from phi_q_odd to smaller expo = exponent of the group all character orders divide this number, and a character of that order exists - use conrey logarithm to reuse log and to loop efficiently over the group (see the diff on l.c, only 1 log in computed instead of 2 * q) - NOT TESTED, for the moment it just compiles, I know some errors (e.g. the FIXME in group_init.c : the generators have to be lifted mod q) this commit is just a proof of concept.
2016-02-24 11:06:49 +01:00
update licence
2016-09-06 15:32:23 +02:00
This file is part of Arb.
Arb is free software: you can redistribute it and/or modify it under
the terms of the GNU Lesser General Public License (LGPL) as published
by the Free Software Foundation; either version 2.1 of the License, or
(at your option) any later version. See <http://www.gnu.org/licenses/>.
*/
[dirichlet] add conrey type to handle logs + char type - try to handle even and odd components the same way in the dirichlet group - switch from phi_q_odd to smaller expo = exponent of the group all character orders divide this number, and a character of that order exists - use conrey logarithm to reuse log and to loop efficiently over the group (see the diff on l.c, only 1 log in computed instead of 2 * q) - NOT TESTED, for the moment it just compiles, I know some errors (e.g. the FIXME in group_init.c : the generators have to be lifted mod q) this commit is just a proof of concept.
2016-02-24 11:06:49 +01:00
#include "acb_dirichlet.h"
change interface to switch acb / ui output
2016-04-05 12:18:39 +02:00
void
update doc and change chi_char to pairing_char
2016-10-08 22:43:29 +02:00
acb_dirichlet_pairing_char(acb_t res, const dirichlet_group_t G, const dirichlet_char_t a, const dirichlet_char_t b, slong prec)
[dirichlet] add conrey type to handle logs + char type - try to handle even and odd components the same way in the dirichlet group - switch from phi_q_odd to smaller expo = exponent of the group all character orders divide this number, and a character of that order exists - use conrey logarithm to reuse log and to loop efficiently over the group (see the diff on l.c, only 1 log in computed instead of 2 * q) - NOT TESTED, for the moment it just compiles, I know some errors (e.g. the FIXME in group_init.c : the generators have to be lifted mod q) this commit is just a proof of concept.
2016-02-24 11:06:49 +01:00
{
change interface to switch acb / ui output
2016-04-05 12:18:39 +02:00
ulong expo;
update doc and change chi_char to pairing_char
2016-10-08 22:43:29 +02:00
expo = dirichlet_pairing_char(G, a, b);
separate dirichlet/acb_dirichlet
2016-10-06 14:19:39 +02:00
if (expo == DIRICHLET_CHI_NULL)
change interface to switch acb / ui output
2016-04-05 12:18:39 +02:00
acb_zero(res);
else
{
fmpq_t t;
fmpq_init(t);
fmpq_set_si(t, 2 * expo, G->expo);
arb_sin_cos_pi_fmpq(acb_imagref(res), acb_realref(res), t, prec);
fmpq_clear(t);
}
[dirichlet] add conrey type to handle logs + char type - try to handle even and odd components the same way in the dirichlet group - switch from phi_q_odd to smaller expo = exponent of the group all character orders divide this number, and a character of that order exists - use conrey logarithm to reuse log and to loop efficiently over the group (see the diff on l.c, only 1 log in computed instead of 2 * q) - NOT TESTED, for the moment it just compiles, I know some errors (e.g. the FIXME in group_init.c : the generators have to be lifted mod q) this commit is just a proof of concept.
2016-02-24 11:06:49 +01:00
}
Reference in a new issue Copy permalink