| Karim Belabas on Mon, 05 Apr 2004 21:01:03 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| Re: patch to speed up binomial |
* Bill Allombert [2004-03-20 14:54]:
> Hello PARI developers,
>
> Here a patch to speed up binomial() the same way as mpfact:
>
> Basically it use divide_conquer_prod. That should be around
> 60% faster on my tests.
>
> This patch can probably be made cleaner.
> I don't think it affects the accuracy of the result in any way
> (when n is a real).
Applied in CVS after further small improvements ( use diviiexact(),
use seq_umul() for small inputs. )
Thanks !
Karim.
--
Karim Belabas Tel: (+33) (0)1 69 15 57 48
Dep. de Mathematiques, Bat. 425 Fax: (+33) (0)1 69 15 60 19
Universite Paris-Sud http://www.math.u-psud.fr/~belabas/
F-91405 Orsay (France) http://pari.math.u-bordeaux.fr/ [PARI/GP]