#20439: eigenmatrix_right gives the conjugate of what it should
------------------------------+----------------------------
Reporter: tmonteil | Owner:
Type: defect | Status: new
Priority: major | Milestone: sage-7.2
Component: linear algebra | Keywords:
Merged in: | Authors:
Reviewers: | Report Upstream: N/A
Work issues: | Branch:
Commit: | Dependencies:
Stopgaps: |
------------------------------+----------------------------
Here is the example reported at http://ask.sagemath.org/question/33084
/bug-in-eigenmatrix-command/
{{{
sage: A = matrix(CDF, [[-2.53634347567, 2.04801738686, -0.0,
-62.166145304], [ 0.7, -0.6, 0.0, 0.0], [0.547271128842, 0.0, -0.3015,
-21.7532081652], [0.0, 0.0, 0.3, -0.4]])
sage: D, P = A.eigenmatrix_right()
}}}
According to the documentation, `D,P` should satisfy `A*P == P*D`,
however:
{{{
sage: (A*P - P*D).norm()
7.454841195108627
}}}
The conjugate of `P` seems to be the correct answer though:
{{{
sage: P = P.conjugate()
sage: (A*P - P*D).norm()
6.716506829378007e-15
}}}
--
Ticket URL: <http://trac.sagemath.org/ticket/20439>
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 https://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
