#12521: evaluate log gamma for complex input
------------------------------------+-----------------------------
Reporter: kcrisman | Owner: AlexGhitza
Type: defect | Status: needs_work
Priority: critical | Milestone: sage-5.13
Component: basic arithmetic | Resolution:
Keywords: lgamma log_gamma | Merged in:
Authors: Eviatar Bach | Reviewers: Burcin Erocal
Report Upstream: N/A | Work issues:
Branch: | Commit:
Dependencies: | Stopgaps:
------------------------------------+-----------------------------
Changes (by jdemeyer):
* status: needs_review => needs_work
Comment:
Also, it seems the algorithm for complex numbers is much faster than for
real numbers, so we should probably use the complex algorithm for all
inputs:
{{{
sage: a=RealField(53)(2.5); timeit('log_gamma(a)')
625 loops, best of 3: 118 µs per loop
sage: a=ComplexField(53)(2.5); timeit('log_gamma(a)')
625 loops, best of 3: 112 µs per loop
}}}
{{{
sage: a=RealField(10**4)(2.5); timeit('log_gamma(a)')
5 loops, best of 3: 506 ms per loop
sage: a=ComplexField(10**4)(2.5); timeit('log_gamma(a)')
625 loops, best of 3: 782 µs per loop
}}}
--
Ticket URL: <http://trac.sagemath.org/ticket/12521#comment:21>
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/groups/opt_out.
