On Wed, May 02, 2007 at 12:38:12AM +0200, Guillem Borrell i Nogueras wrote:
lapack.argtypes=[c_int,c_int, ndpointer(dtype=float64, ndim=2, flags='FORTRAN'), c_int,c_int, ndpointer(dtype=float64, ndim=1, flags='FORTRAN'), c_int,c_int]
This also isn't correct, according to the dgesv documentation. It should be lapack.dgesv_.argtypes=[POINTER(c_int),POINTER(c_int), ndpointer(dtype=np.float64, ndim=2, flags='FORTRAN'), POINTER(c_int), POINTER(c_int), ndpointer(dtype=np.float64, ndim=2, flags='FORTRAN'), POINTER(c_int),POINTER(c_int)] I attach a working version of your script. Cheers Stéfan