[SciPy-dev] Scipy cannot handle singular eigenvalue problems

[Nils Wagner]

Hi all,

Scipy cannot handle singular eigenvalue problems. Note, that both
matrices
are singular.

from scipy import *
A = array(([1.,-1.],[-1.,1.]))
B = array(([1.0,0.0],[0.0,0.0]))
w,vr = linalg.eig(A,B)
print w
print vr

Traceback (most recent call last):
  File "sevp.py", line 4, in ?
    w,vr = linalg.eig(A,B)
  File "/usr/lib/python2.3/site-packages/scipy/linalg/decomp.py", line
109, in eig
    return _geneig(a1,b,left,right,overwrite_a,overwrite_b)
  File "/usr/lib/python2.3/site-packages/scipy/linalg/decomp.py", line
68, in _geneig
    vr = _make_complex_eigvecs(w, vr, t)
  File "/usr/lib/python2.3/site-packages/scipy/linalg/decomp.py", line
29, in _make_complex_eigvecs
    vnew.real = scipy_base.take(vin,ind[::2],1)
ValueError: matrices are not aligned for copy

Any suggestion ?

The output of MATLAB is enclosed as well (out.ps)

Nils
-------------- next part --------------
A non-text attachment was scrubbed...
Name: out.ps
Type: application/postscript
Size: 4412 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20031210/59143b01/attachment.ps>