#7099: serious incomplete gamma function precision bugs
-------------------------------------+-------------------------------------
Reporter: was | Owner: jkantor
Type: defect | Status: needs_review
Priority: critical | Milestone: sage-6.3
Component: numerical | Resolution:
Keywords: | Merged in:
Authors: Peter Bruin | Reviewers:
Report Upstream: N/A | Work issues:
Branch: | Commit:
u/pbruin/7099-incomplete_gamma | 78ba9638700d8cd9d4b18d83323af8fe9507a7ba
Dependencies: | Stopgaps:
-------------------------------------+-------------------------------------
Comment (by rws):
With #16697 the output is as expected:
{{{
sage: C = ComplexField(1000)
sage: gamma_inc(C(2+I),C(3+I))
0.121515644664508695525971545977439666159749344176962379708992904126499444842886620664991650378432544392118359044438541514683402245033018771644222346410367471459456844674335147722343580581945662693850674590491020834632434082710800093315646442975240326569517738365018117780134100101636704042869033248174
+
0.101533909079826033296475736021224621546966200987295663190553587086145836461236284668967411665020429964946098113930918849948956425984499549094441904693395768367238320065064071027383069839637218088862214571990869510941211277488169032567679631037683814516738122300220474252081775895835843619616213883517*I
}}}
This is because I changed the default to mpmath for the gamma functions.
The ticket is still relevant for `C(2+I).gamma_inc(C(3+I))` however.
--
Ticket URL: <http://trac.sagemath.org/ticket/7099#comment:18>
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.
