On Saturday, 4 February 2012 11:53:35 UTC+8, Dr. David Kirkby wrote: > > > > On 2 February 2012 23:05, rjf <> wrote: > >> >> I don't know about arithmetic on ARM specifically, but there is >> something >> wrong with a gamma() function that fails to return an integer (perhaps >> in float format) when it is >> given an integer argument (perhaps in float format), and the answer is >> exactly representable >> as an integer in float format. >> > > Would you say this is true for any function foobar(), or gamma() in > particular? >
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 :) -- 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
