Comment #1 on issue 1664 by aaron.j.gardner: Matrix([ [1, 2, 3], [4, 5, 6],
[7, 8, 9] ]).eigenvects() fails
http://code.google.com/p/sympy/issues/detail?id=1664
This 3x3 matrix is not the only one producing an eigenvects() error.
I discovered that the problem exists for some 2x2 matrices I've been
working with in
my differential enqs class.
Examples:
=================
>>> M = Matrix([ [1,-1], [1,3] ])
>>> M.eigenvals()
{2: 2}
>>> M.eigenvects()
Traceback (most recent call last):
File "<console>", line 1, in <module>
File "./sympy/matrices/matrices.py", line 1426, in eigenvects
assert len(basis) == k
AssertionError
=================
>>> A = Matrix([ [-2,-3], [-4,1] ])
>>> A.eigenvals()
____ ____
\/ 57 \/ 57
{-1/2 + ------: 1, -1/2 - ------: 1}
2 2
>>> A.eigenvects()
Traceback (most recent call last):
File "<console>", line 1, in <module>
File "./sympy/matrices/matrices.py", line 1426, in eigenvects
assert len(basis) == k
AssertionError
I think the problem lies in the calculations done by nullspace() when
eigenvects()
is called. When I can devote more time, I'll see what I can come up with.
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to sympy-issues@googlegroups.com
To unsubscribe from this group, send email to
sympy-issues+unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en
-~----------~----~----~----~------~----~------~--~---
