mirror of
https://github.com/vale981/arb
synced 2025-03-06 09:51:39 -05:00
44 lines
1.1 KiB
C
44 lines
1.1 KiB
C
/*
|
|
Copyright (C) 2016 Pascal Molin
|
|
|
|
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/>.
|
|
*/
|
|
|
|
#include "acb_poly.h"
|
|
#include "acb_dirichlet.h"
|
|
|
|
/* all roots are already computed */
|
|
void
|
|
_acb_dirichlet_dft_pol(acb_ptr w, acb_srcptr v, acb_srcptr z, slong len, slong prec)
|
|
{
|
|
/* FIXME: huge accuracy loss */
|
|
#if 0
|
|
_acb_poly_evaluate_vec_fast(w, v, len, z, len, prec);
|
|
#elif 0
|
|
_acb_poly_evaluate_vec_iter(w, v, len, z, len, prec);
|
|
#else
|
|
slong i, j;
|
|
|
|
for (i = 0; i < len; i++)
|
|
{
|
|
acb_zero(w + i);
|
|
for (j = 0; j < len; j++)
|
|
acb_addmul(w + i, v + j, z + (i * j % len), prec);
|
|
}
|
|
#endif
|
|
}
|
|
|
|
void
|
|
acb_dirichlet_dft_pol(acb_ptr w, acb_srcptr v, slong len, slong prec)
|
|
{
|
|
acb_ptr z;
|
|
z = _acb_vec_init(len);
|
|
acb_dirichlet_vec_nth_roots(z, len, prec);
|
|
_acb_dirichlet_dft_pol(w, v, z, len, prec);
|
|
_acb_vec_clear(z, len);
|
|
}
|