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]