Status: New
Owner: [email protected]
Labels: Type-Defect Priority-Medium

New issue 2193 by [email protected]: Matrix.eigenvects() return empty vectors in some cases
http://code.google.com/p/sympy/issues/detail?id=2193

In some cases (where matrix is symbolical or with complex) "eigenvects" method return no vectors,

m = Matrix(2,2,[1, 0, 0, I])
m.eigenvects()
[(1/2 + I/2 - (-1)**(3/4)*2**(1/2)/2, 1, []), (1/2 + I/2 + (-1)**(3/4)*2**(1/2)/2, 1, [])]

expected: [(1, 1, [1, 0]), (I, 1, [0, 1])]

The same is with
a = var("a")
m = Matrix(2,2,[1, 0, 0, a])
m.eigenvects()

Probably it is because of the very complex answer ("1/2 + I/2 - (-1)**(3/4)*2**(1/2)/2") given by roots(x**2 + (1 + I)*x + I, x), while must be symply 1 and I.

Also, in any case, if no vectors are found for eigenvalue (by unrecognized reasons) but they must be, then exception must be raised in this method at the right moment.

--
You received this message because you are subscribed to the Google Groups 
"sympy-issues" 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/sympy-issues?hl=en.

Reply via email to