[SciPy-dev] isn't it a bug in scipy<->LAPACK connection?

dmitrey dmitrey.kroshko@scipy....
Thu Oct 18 12:06:56 CDT 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