Le 04/02/2012 06:00, Dima Pasechnik a écrit :
It's well-known how to compute gamma() better than it is implemented in
(e)glibc, the prevalent Linux libc implementation,
which computes exp(lgamma()) rather than gamma() directly.
See e.g.
http://oai.cwi.nl/oai/asset/10080/10080A.pdf
Perhaps I should give it to a student, to re-write in C or Cython :)
1. Do the libc used in the *BSD family of systems do better?
2. If there are better algorightms, they may not be known to the various
libc authors, and they'd probably gladly accept patches.
Snark on #sagemath
--
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
