[SciPy-user] generalized eig algorithm did not converge

Nils Wagner nwagner at mecha.uni-stuttgart.de
Mon Mar 20 09:20:36 CST 2006


Hi all,

What is the reason for this message ?

Traceback (most recent call last):
  File "gh.py", line 4, in ?
    w = linalg.eigvals(G,H)
  File "/usr/lib64/python2.4/site-packages/scipy/linalg/decomp.py", line
244, in eigvals
    return eig(a,b=b,left=0,right=0,overwrite_a=overwrite_a)
  File "/usr/lib64/python2.4/site-packages/scipy/linalg/decomp.py", line
121, in eig
    return _geneig(a1,b,left,right,overwrite_a,overwrite_b)
  File "/usr/lib64/python2.4/site-packages/scipy/linalg/decomp.py", line
72, in _geneig
    if info>0: raise LinAlgError,"generalized eig algorithm did not
converge"
scipy.linalg.basic.LinAlgError: generalized eig algorithm did not converge

Matlab is able to solve this generalized eigenvalue problem without any
warning.

Afaik, Matlab and scipy's eig is based on LAPACK. So what is the reason
for the difference ?

Any idea ?

Nils

-------------- next part --------------
A non-text attachment was scrubbed...
Name: mat.tar.gz
Type: application/x-tar
Size: 20590 bytes
Desc: not available
Url : http://www.scipy.net/pipermail/scipy-user/attachments/20060320/f43be186/mat.tar-0001.tar
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gh.m
Type: application/m-file
Size: 51 bytes
Desc: not available
Url : http://www.scipy.net/pipermail/scipy-user/attachments/20060320/f43be186/gh-0001.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gh.py
Type: text/x-python
Size: 90 bytes
Desc: not available
Url : http://www.scipy.net/pipermail/scipy-user/attachments/20060320/f43be186/gh-0001.py


More information about the SciPy-user mailing list