Well, i am starting with MKL, using C++, and I test the "dsyevr " of Lapack, to compute the eigenvalues of a matrix, but i was "shock" when I try a matrix with dimension bigger or equal than 710, after that the program said "segmentation fault", but when I use a matrix of 709 it works, really good... so my question is that if the command "dsyevr" has a limitation for the size of the matrix to calculate the corresponding EV

The part of the code is :

N=709

char jobz = 'N'; char range = 'A'; char uplo = 'U'; int n = N; double H[N][N]; int lda = N; int lwork = 26*N; double vl, vu; int il, iu; double abstol = 1.0; int m; double w[N]; double z[N][N]; int isuppz; int ldz = N; double work[lwork]; int liwork = 10*N; int iwork[liwork]; int info; dsyevr(&jobz, &range, &uplo, &n, *H, &lda, &vl, &vu, &il, &iu, &abstol, &m, w, *z, &ldz, &isuppz, work, &lwork, iwork,&liwork, &info);

Thanks in advances