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

improve series composition in the linear case

Browse source
This commit is contained in:
Fredrik Johansson 2013-07-19 16:33:14 +02:00
parent 6231061911
commit a6dc4358cc
1 changed files with 24 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

25
fmprb_poly/compose_series.c
View file

@ -29,10 +29,33 @@ void
_fmprb_poly_compose_series(fmprb_ptr res, fmprb_srcptr poly1, long len1,
fmprb_srcptr poly2, long len2, long n, long prec)
{
if (len1 < 6 || n < 6)
if (len2 == 2) /* fast linear case (TODO: all monomials) */
{
long i;
fmprb_t t;
fmprb_init(t);
fmprb_set(t, poly2 + 1);
fmprb_set_round(res, poly1, prec);
for (i = 1; i < n; i++)
{
fmprb_mul(res + i, poly1 + i, t, prec);
if (i + 1 < n)
fmprb_mul(t, t, poly2 + 1, prec);
}
fmprb_clear(t);
}
else if (len1 < 6 || n < 6)
{
_fmprb_poly_compose_series_horner(res, poly1, len1, poly2, len2, n, prec);
}
else
{
_fmprb_poly_compose_series_brent_kung(res, poly1, len1, poly2, len2, n, prec);
}
}
void