#8321: numerical integration with arbitrary precision
-------------------------+--------------------------------------------------
Reporter: burcin | Owner: burcin
Type: defect | Status: needs_work
Priority: major | Milestone: sage-4.7.2
Component: symbolics | Keywords: numerics,integration
Work_issues: | Upstream: N/A
Reviewer: | Author: Stefan Reiterer
Merged: | Dependencies:
-------------------------+--------------------------------------------------
Comment(by kcrisman):
{{{
# f = e^(-x^2)*log(x) on [17, 42]
sage: num_int_test(e^(-x^2)*log(x), 17, 42)
Exact = integrate(e^(-x^2)*log(x), x, 17, 42)
Exact .n() = 2.59225286296247e-127
GSL = 2.5657285007e-127
mpmath = 2.59225286296247e-127
}}}
Gee, I don't like this at all. So do you think this is a bug? (And if
so, in which program?) GSL is pretty stable, though...
{{{
try:
precision = parent.prec()
mpmath.mp.prec = precision
if precision == RDF.precision():
mp_f = fast_callable(f, vars=[x], domain=RDF)
}}}
What if we called GSL for precisely this use case, not
{{{fast_callable}}}? Also, what does Maxima ({{{.nintegrate()}}}) do?
Good analysis!
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/8321#comment:32>
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 post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
