[SciPy-user] Possible bug in linalg.lu
Pietro Berkes
p.berkes at biologie.hu-berlin.de
Mon Jul 14 10:23:15 EDT 2003
the function linalg.lu seems to be defined for real matrices only. if you want
to use it for complex ones, the fastest solution is probably to write a
wrapper to the lapack function ZGETRF. it should work like in the real case
(cfr. linalg/src/lu.f).
regards,
pietro.
On Monday 14 July 2003 15:29, Nils Wagner wrote:
> Hi all,
>
> >>> C
>
> array([[ 3.003+0.j, 1. +0.j, 0. +0.j, 0. +0.j],
> [ 3. +0.j, 3.003+0.j, 2. +0.j, 0. +0.j],
> [ 0. +0.j, 2. +0.j, 3.003+0.j, 3. +0.j],
> [ 0. +0.j, 0. +0.j, 1. +0.j, 2. +0.j]])
>
> >>> linalg.lu(C)
>
> Traceback (most recent call last):
> File "<stdin>", line 1, in ?
> File "/usr/local/lib/python2.1/site-packages/scipy/linalg/decomp.py",
> line 238, in lu
> p,l,u,info = flu(a1,permute_l=permute_l,overwrite_a = overwrite_a)
> TypeError: object of type 'None' is not callable
>
>
> For what reason ?
>
> Nils
>
> _______________________________________________
> 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