[SciPy-user] Iterative solvers and complex arithmetic

[Travis Oliphant]

Nils Wagner wrote:

> Hi all,
>
> Who can resolve this problem ?
>
> python2.3 -i failure.py
>
> Traceback (most recent call last):
>  File "failure.py", line 9, in ?
>    x,info = linalg.gmres(A,r)
>  File "/usr/lib/python2.3/site-packages/scipy/linalg/iterative.py", 
> line 499, in gmres
>    work[slice2] += sclr1*matvec(work[slice1])
> TypeError: return array has incorrect type
>
> Nils

I will try to fix this, for now make sure A and r have the same type.

By the way,  you don't need from RandomArray import *

rand is available from scipy

-Travis