[SciPy-dev] isn't it a bug in scipy<->LAPACK connection?
dmitrey
dmitrey.kroshko at scipy.org
Thu Oct 18 13:06:56 EDT 2007
Hi all,
isn't it a bug in scipy<->LAPACK connection?
from scipy.linalg.flapack import dgelss
#print dgelss.__doc__
from scipy import rand
N = 100
M = 10 #!!with M >= N all works ok!!
A = rand(M,N)
b = 15+rand(M)
v,x,s,rank,info = dgelss(A,b)
So it yields "failed in converting 2nd argument `b' of flapack.dgelss to
C/Fortran array".
dgelss routine allows M to be less, equal or greater to N, you can
easily see it from documentation, for example:
https://web.kudpc.kyoto-u.ac.jp/doc/HPC-WG/Manual/lapack/dgelss.html
Regards, Dmitrey
More information about the SciPy-Dev
mailing list