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]