Did you test it for the negative input as well? Unfortunately, I cannot compile anything after 7.2.beta6 on my machine, so I'm stuck with reading the code---but from what I read, numeric evaluation from beta in pynac calls numeric evaluation of lgamma in pynac, which in turn (correct me if I'm wrong) calls py_lgamma in pynac.pyx, which then just calls x.log_gamma(), without conversion to a complex argument. Then, the documented behavior (I did miss that, sorry) of log_gamma comes into play.
But of course, if everything works in 7.3.beta4, then there is no need to discuss or fix this any further. ;-) Benjamin Am Dienstag, 14. Juni 2016 10:28:04 UTC+2 schrieb Ralf Stephan: > > On Tuesday, June 14, 2016 at 10:18:01 AM UTC+2, Benjamin Hackl wrote: >> >> This seems to come from this (7.2.beta6): >> >> sage: (-0.9).log_gamma() >> NaN >> > > No, this is independent and even documented: > To get evaluation of input for which gamma > is negative and the ceiling is even, we must > explicitly make the input complex. This is > a known issue, see :trac:`12521`:: > > sage: log_gamma(-2.1) > NaN > sage: log_gamma(CC(-2.1)) > 1.53171380819509 + 3.14159265358979*I > > The beta() issue works in 7.3beta4 so where is the problem exactly? > -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To post to this group, send email to sage-devel@googlegroups.com. Visit this group at https://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/d/optout.
