C/C++ interface

[csv610]

Hello,

I am now using Intel MKL instead of public domain Clapack.

Can someone tell me why I get the following error when I run the program: This error doesn't come with

CLAPACK.

** On entry to DGESV , parameter number 4 had an illegal value

typedef MKL_INT integer;

integer nrhs = 1;
integer lda = m;
integer ldb = n;
integer info;
integer *ipiv;

ipiv = ( integer *) malloc( m*sizeof(integer) );

// A = gen_gen_matrix(m,n);
// b = gen_gen_matrix(m, nrhs);

integer M = m;
dgesv_(&M, &nrhs, A, &lda, ipiv, b, &ldb, &info);

[Michael_C_Intel4]

Hello.

If you storethe matrix Arow-wise you need to provide lda=n. Probaly, it would help.

Michael.