Community
cancel
Showing results for 
Search instead for 
Did you mean: 
68 Views

Problem in using dgeev

Dear all,

I am trying to calculate only the eigenvalues of a general 33 matrix. I don't need the eigenvectors. Please find below my code.

[cpp]#define LAPACK_N 3 #define LAPACK_LDA LAPACK_N #define LAPACK_LDVL 1 #define LAPACK_LDVR 1 // Locals. int n = LAPACK_N; int lda = LAPACK_LDA, ldvl = LAPACK_LDVL, ldvr = LAPACK_LDVR; int info, lwork; double wkopt; double* work; double a[ LAPACK_LDA*LAPACK_N ] = { 10.01, 3.98, 3.30, 0.86, 1.53, 8.26, 4.60, 7.04, 3.89 }; // Local arrays. double wr[LAPACK_N], wi[LAPACK_N]; double vl[LAPACK_LDVL*LAPACK_N], vr[LAPACK_LDVR*LAPACK_N]; // Query and allocate the optimal workspace. lwork = -1; dgeev( "N", "N", &n, a, &lda, wr, wi, vl, &ldvl, vr, &ldvr, &wkopt, &lwork, &info ); lwork = (int) wkopt; work = (double*) malloc( lwork*sizeof(double) ); // Solve eigenproblem. dgeev( "N", "N", &n, a, &lda, wr, wi, vl, &ldvl, vr, &ldvr, work, &lwork, &info );[/cpp]

When the compiler reaches the first dgetev function on line 23, I get the following error:
Unhandled exception at 0x000000013f14a35f in mycode.exe: 0xC0000005: Access violation reading location 0x0000000800169eb0. Where the address of a is 0x0000000000169eb0. I am using Intel MKL 10.3 in Visual Studio 2008 in Windows 7 64bit.

I would be thankful if someone could help me.
0 Kudos
1 Reply
68 Views

I found the solution. It was because of wrong integer size.
Reply