[SciPy-user] Problems with complex matrices

Travis Oliphant oliphant at ee.byu.edu
Fri Oct 22 15:20:02 CDT 2004


Nils Wagner wrote:

> Hi all,
>
> Another bug is revealed
>
> Traceback (most recent call last):
>  File "complex_it.py", line 6, in ?
>    res=dot(A,x)-r
>  File "/usr/lib/python2.3/site-packages/Numeric/dotblas/__init__.py", 
> line 55, in dot
>    if multiarray.array(a).shape == () or multiarray.array(b).shape == ():
> TypeError: a float is required
>
> Nils
>
>------------------------------------------------------------------------
>
>from scipy import *
>from RandomArray import *
>A = rand(10,10)+1j*rand(10,10)
>r = rand(10)+1j*rand(10)
>x = linalg.gmres(A,r)
>  
>
should be

x, info = linalg.gmres(A,r)

>res=dot(A,x)-r
>  
>
>------------------------------------------------------------------------
>
>_______________________________________________
>SciPy-user mailing list
>SciPy-user at scipy.net
>http://www.scipy.net/mailman/listinfo/scipy-user
>  
>



More information about the SciPy-user mailing list