#12595: eigenmatrix of complex floating-point matrix is wrong
------------------------------------+---------------------------------------
Reporter: dkrenn | Owner: jason, was
Type: defect | Status: new
Priority: major | Milestone: sage-5.1
Component: linear algebra | Resolution:
Keywords: eigenmatrix, CDF | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
------------------------------------+---------------------------------------
Comment (by dsm):
Somewhat smaller example, and evidence that something's wrong in RDF too:
{{{
sage: M = [[1,-1],[1, 1]]
sage: rings = SR, ZZ, QQ, RDF, CDF
sage: for ring in rings:
....: m = Matrix(ring, M)
....: d, p = m.left_eigenmatrix()
....: print m- (~p)*d*p
....:
[0 0]
[0 0]
[0 0]
[0 0]
[0 0]
[0 0]
[ 0.0 -2.0]
[ 2.0 0.0]
[2.22044604925e-16 + 1.96068977308e-16*I -2.0 -
2.22044604925e-16*I]
[ 2.0 + 8.14867788484e-17*I
-2.22044604925e-16]
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/12595#comment:2>
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.
