#4942: find_root() is broken when interval borders cannot be evaluated
-----------------------+----------------------------------------------------
Reporter: mabshoff | Owner: mhansen
Type: defect | Status: assigned
Priority: blocker | Milestone: sage-3.3
Component: numerical | Resolution:
Keywords: |
-----------------------+----------------------------------------------------
Comment (by mhansen):
It seems this is a problem with Scipy:
{{{
In [16]: def f(x):
....: return 1.0/(x-1.0)+1.0
....:
In [17]: import scipy.optimize
In [18]: scipy.optimize.brentq(f, 0, 2)
Out[18]: 0.0
In [19]: f(0.001)
Out[19]: -0.0010010010010010895
In [20]: f(2)
Out[20]: 2.0
In [21]: scipy.optimize.brentq(f, 0.001, 2)
Out[21]: 1.0000000000007283
In [22]: f(1.0000000000007283)
Out[22]: 1373048666882.2488
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4942#comment:4>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
