#5378: Sage 3.3: numerical noise in rings/polynomial/complex_roots.py on cicero
&
fulvia
----------------------+-----------------------------------------------------
Reporter: mabshoff | Owner: mabshoff
Type: defect | Status: new
Priority: blocker | Milestone: sage-3.4.1
Component: doctest | Keywords:
----------------------+-----------------------------------------------------
Comment(by mabshoff):
The following illustrates the problem and a potential solution:
{{{
----------------------------------------------------------------------
| Sage Version 3.4.1.rc3, Release Date: 2009-04-16 |
| Type notebook() for the GUI, and license() for information. |
----------------------------------------------------------------------
Loading Sage library. Current Mercurial branch is: review
sage: from sage.rings.polynomial.complex_roots import complex_roots
sage: x = polygen(ZZ)
sage: complex_roots(x^2 + 27*x + 181)
[(-12.3819660112501?, 1), (-14.61803398874990? + 0.?e-27*I, 1)]
sage: v=complex_roots(x^2 + 27*x + 181)
sage: sorted((v[0][0].real(),v[1][0].real()))
[-14.61803398874990?, -12.3819660112501?]
}}}
On another machine we get:
{{{
sage: from sage.rings.polynomial.complex_roots import complex_roots
sage: x = polygen(ZZ)
sage: complex_roots(x^2 + 27*x + 181)
[(-14.61803398874990? + 0.?e-27*I, 1), (-12.38196601125010? + 0.?e-27*I,
1)]
sage: v=complex_roots(x^2 + 27*x + 181)
sage: sorted((v[0][0].real(),v[1][0].real()))
[-14.61803398874990?, -12.38196601125010?]
}}}
Patch coming up.
Cheers,
Michael
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/5378#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
-~----------~----~----~----~------~----~------~--~---
