#17328: incomplete gamma function bugs for certain arguments
-------------------------------------+-------------------------------------
Reporter: kcrisman | Owner:
Type: defect | Status: needs_review
Priority: major | Milestone: sage-6.4
Component: calculus | Resolution:
Keywords: incomplete gamma | Merged in:
function | Reviewers:
Authors: Peter Bruin | Work issues:
Report Upstream: N/A | Commit:
Branch: | 3cae7c8469d6c2f5dd3d15a872965f528f637c5c
u/pbruin/17328-incgam | Stopgaps:
Dependencies: |
-------------------------------------+-------------------------------------
Comment (by pbruin):
Replying to [comment:8 kcrisman]:
> Seems to give correct answers in various precisions.
> {{{
> plot(lambda t: gamma(-1, -log(t)).real() -
numerical_integral(1/ln(x)^2,2,t)[0],2,3,color='red')
> }}}
> is a very nice straight line ''just'' below zero now. And the
imaginary part seems to be consistently pi for t>1. And the fix is a good
one.
OK, thanks!
> > We should also add a doctest for the original bug.
> Did you want to add
> {{{
> N(integral(1/log(x)^2,(x,2,3)))
> }}}
> as well then, or is the example you have with the actual antideriv and
floats sufficient?
I don't think it is necessary to add another doctest for the integral;
this would just test in addition that Maxima computes the correct
antiderivative. Of course the above command was how the bug was
originally found, but I don't really see why we should doctest exactly
that command. What I meant to say is "add a doctest to show that the PARI
bug has been fixed".
--
Ticket URL: <http://trac.sagemath.org/ticket/17328#comment:10>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
