[SciPy-user] Field report of iterative solvers in scipy and a possibly bug in mmio.py

[Nils Wagner]

Hi all,

I am wondering if someone is able to reproduce the frequency response 
function
of the test set available through

http://www.imtek.uni-freiburg.de/simulation/benchmark/?li=30&nr=38886

with the currently implemented iterative solvers in scipy ?

Any field report would be appreciated.

Thanks in advance

                                       Nils
from scipy import *
print 'Reading stiffness matrix'
K=io.mmread('windscreen.K.mtx')
print 'Reading mass matrix'
M=io.mmread('windscreen.M.mtx')
print 'Reading load vector'
f=io.mmread('windscreen.B.mtx')
omega = 0.5*arange(1,401)
for i in omega:
  x,info = linalg.bicgstab(K-i**2*M,f)

There seems to be a bug in mmio.py

Traceback (most recent call last):
  File "frf.py", line 7, in ?
    f=io.mmread('windscreen.B.mtx')
  File "/usr/lib/python2.3/site-packages/scipy/io/mmio.py", line 103, in 
mmread
    rows,cols,entries,rep,field,symm = mminfo(source)
  File "/usr/lib/python2.3/site-packages/scipy/io/mmio.py", line 65, in 
mminfo
    assert len(line)==2,`line`
AssertionError: ['22692', '1', '22692']