dstedc and workspace array

Open discussion regarding features, bugs, issues, vendors, etc.

dstedc and workspace array

Postby treq11 » Mon Mar 30, 2009 3:15 pm


I'm running into some errors using the routine dstedc. After querying for the workspace sizes and executing dstedc my program crashes. The c++ code I'm running is basically the following (for several different values of n)

Code: Select all
char type = 'V';
int info;
int n = 30;
double* work = new double[lwork];  // lwork = 3091
double* iwork = new int[liwork];     // liwork = 936
// eVals   - double[n]
// subDiag - double[n-1]
// eVecs   - double[n*n]

dstedc( &type, &n, eVals, subDiag, eVecs, &n, work, &lwork, iwork, &liwork, &info );

If I increase lwork and liwork the routine works properly. This happens with MKL (10.1) as well as the netlib version (3.1). ACML (4.2.0), which doesn't require inputting the workspace sizes, gives a similar error.

Thanks in advance for any help!
Posts: 3
Joined: Mon Mar 30, 2009 1:34 pm

Re: dstedc and workspace array

Postby treq11 » Wed May 06, 2009 12:01 pm

Any idea why the problem above happens?
Posts: 3
Joined: Mon Mar 30, 2009 1:34 pm

Re: dstedc and workspace array

Postby Julien Langou » Mon May 11, 2009 3:18 pm

I did not have time to check your bug report. I took the time though to add it to the list of Errata of LAPACK-3.2.
See: http://www.netlib.org/lapack/Errata/
Hopefully someone will have time to investigate.
Thanks for reporting this.
Best wishes,
Julien Langou
Posts: 821
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA

Re: dstedc and workspace array

Postby CyLith » Wed Aug 04, 2010 2:56 pm

I cannot reproduce this. The code above uses a double (instead of int) pointer for iwork. If that is actually the case in the real code you will get a buffer overflow.
Posts: 40
Joined: Sun Feb 08, 2009 7:23 am
Location: Stanford, CA

Re: dstedc and workspace array

Postby admin » Wed Jun 22, 2011 7:35 am

We may have found your problem!!!
The fix was committed into the LAPACK repository.
See http://www.netlib.org/lapack/#_svn_access to access it
Site Admin
Posts: 608
Joined: Wed Dec 08, 2004 7:07 pm

Re: dstedc and workspace array

Postby rmjarvis » Sun Feb 05, 2012 3:17 pm

I ran across the same problem in MKL with icpc 12.0. I also have the workaround. The recommended size of the lwork is listed as 1+3n+2n*lgn+3n^2. However, this is too low by n^2. The complex varieties get it right, since they have the same recommendation for lrwork, plus an additional n^2 in lwork. For real varieties, the right value is the sum of these two: 1+3n+2*lgn+4n^2. So the bug is that a workspace query returns the wrong value in work[0]. And also that a non-query doesn't check for the right value. Presumably that is what Julie is reporting as fixed, but until the fix makes its way into things like MKL and ACML, hopefully this post will help people who may be searching for a solution.
Posts: 1
Joined: Sun Feb 05, 2012 3:10 pm

Re: dstedc and workspace array

Postby zotkevich » Mon Feb 06, 2012 11:49 pm

Intel® MKL documentation contains the fix already. See http://software.intel.com/sites/product ... /index.htm for details.
Posts: 13
Joined: Tue Dec 21, 2010 2:41 am
Location: Intel Co.

Return to User Discussion

Who is online

Users browsing this forum: Google [Bot], Yahoo [Bot] and 4 guests