[Numpy-discussion] lapack_lite dgesv

[R. David]

Hello,

I am trying to use the lapack_lite dgesv routine.

The following sample code :

from numpy import *
[....]
a=zeros((nbrows,nbcols),float,order='C')
[....]
ipiv=zeros((DIM),int,order='C')
[....]
linalg.lapack_lite.dgesv(DIM,1,a,DIM,asarray(ipiv),b,DIM,info)

leads do the followin error message :
lapack_lite.LapackError: Parameter ipiv is not of type PyArray_INT in lapack_lite.dgesv

I don't understand the type problem for ipiv !
Indeed, the type of 'a' is OK, ipiv is created the same way than a, but something goes
wrong.
Do you have a clue for this ?

	Regards,
	Romaric
-- 
--------------------------------------
   R. David - david at icps.u-strasbg.fr
   Tel. : 03 90 24 45 48  (Fax 45 47)
--------------------------------------