Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- Lapack Diagonalization for Large Matrix

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

cd1054

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-18-2020
05:22 PM

153 Views

Hi,

I am working on finding the eigenvalues and eigenvectors of a matrix with dimension 52920. However, when I run the LAPACK routine, my program does not stop running after 16 hours. Do you know what might be wrong? I think there might be an infinite loop that is being entered, but it is difficult for me to track down the problem.

Here is part of my code for reference. Any help would be appreciated! Thank you.

void diagonalize(int D, double* ground_state, double* hamiltonian, double* z, double* w, int* ifail) {

cout << "dimension: " << D << endl;

int matrix_layout = LAPACK_ROW_MAJOR;

int dim = D;

int lda = D;

int ldz = D;

int il;

int iu;

long long int info;

int m;

double abstol = -1;

double vl = 0.0;

double vu = 100.0;

cout << "calling LAPACK..." << endl;

info = LAPACKE_dsyevx(matrix_layout,'V','A','L', dim, hamiltonian, lda, vl, vu, il, iu, abstol, &m, w, z, ldz,ifail);

Link Copied

Accepted Solutions

cd1054

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-19-2020
09:44 AM

85 Views

Thank you!

It turns out that I was just being slightly impatient. The routine completed successfully after ~20 hours.

-Christina

2 Replies

PrasanthD_intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-18-2020
11:28 PM

111 Views

Hi,

There is a dedicated forum for mkl issues(https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/bd-p/oneapi-math-kernel-library), and we are redirecting issue to that forum.

Regards

Prasanth

cd1054

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

11-19-2020
09:44 AM

86 Views

Thank you!

It turns out that I was just being slightly impatient. The routine completed successfully after ~20 hours.

-Christina

For more complete information about compiler optimizations, see our Optimization Notice.