Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- bug in MKL 11.0 LAPACK function 'gels '

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

Highlighted
##

dan0112

Beginner

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

10-16-2012
10:05 AM

16 Views

bug in MKL 11.0 LAPACK function 'gels '

I am attaching a minimal Fortran program that links to the Intel MKL to make use of the LAPACK function gels, which solves a linear least squares problem.

When I compile and run this program with Intel Fortran Composer XE 2013, the output is totally wrong. Compiling and running with Intel Fortran Composer XE 2011 gives the correct results. The output from both is attached.

The very strange thing is that the results are correct on both compilers if you change line 3 of the code to nt = 513 or any smaller number, whereas they stay incorrect for any nt >= 514. This is not due to weird values - the values being read in are all real, finite numbers.

The output also shows the compilation and linking lines (I used the most recent MKL link line helper). There you see the exact version numbers of the compilers. I should note that the Composer XE 2013 runs on Ubuntu 12.04, whereas the Composer XE 2011 runs on CentOS 6 (and used to run on Ubuntu 10.10, also giving the correct results). So it could be an issue with Ubuntu 12.04.

Let me know if I can be of any help.

EDIT: If you call gelsy instead of gels, then the results are correct. I updated the main.f90, so there is a logical 'use_gelsy_instead', and also it now outputs the info flag of the gels command.

Best,

Daniel

12 Replies

Highlighted
##

Daniel,
the same results I see on Windows 64 bit with the latest 11.0 ( bundled with Composer 2013).
That's not clear to me how did you decided that that result is wrong ( -39507.9652423625 -20829.7946097958 ) ?

Gennady_F_Intel

Moderator

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

10-16-2012
10:38 AM

16 Views

Highlighted
##

Gennady,
in the last 2 years I ran this code on various machines / architectures, so I immediately see it is wrong. Here is an intuitive way to see it: series.txt contains data points that lie (strictly) between 0.0 and 20.0 (if I remember the upper bound correctly). The result can be interpreted as the coefficients of a statistical Least Squares regression - correct me if I am wrong. They cannot be all negative and of this large, absolute value, if both the dependent and independent variables lie in (0.0, 20.0), because then the estimated/forecast y-values would always be negative.

dan0112

Beginner

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

10-16-2012
10:56 AM

16 Views

Highlighted
##

Gennady_F_Intel

Moderator

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

10-16-2012
11:01 AM

16 Views

Ok, thanks for clarification.
We will check the code and let you know the results.
--Gennady

Highlighted
##

Daniel,yes, this the problem.
We confirmed the issue. the issue is escalated. We will let you know if any updates.
--Gennady

Gennady_F_Intel

Moderator

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

10-17-2012
03:46 AM

16 Views

Highlighted
##

I can reproduce the errors on our cluster (running CentOS 6), where I can load Intel Composer XE 2012 or 2013, each with corresponding MKL. So this is not an issue with Ubuntu or with my machine.

dan0112

Beginner

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

10-17-2012
06:56 AM

16 Views

Highlighted
##

dan0112

Beginner

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

10-24-2012
06:26 AM

16 Views

Gennady,
yes, please let me know of updates, then I will try to test as quick as possible.

Highlighted
##

MKL 11.0 update 1 just came out. Is this issue fixed? I have to know before I request to update our cluster... If it is not fixed, could you provide me with a way to track the bug?

dan0112

Beginner

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

10-24-2012
11:45 PM

16 Views

Highlighted
##

no, the fix of the problem is not available in update1. The fix is targeted to the next update ( 11.0 update 2).

Gennady_F_Intel

Moderator

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

10-25-2012
03:12 AM

16 Views

Highlighted
##

As I just found out, the problem doesn't appear when I link to a 'parallel' version. Specifically, in my original program (not the toy program posted here), the following link options do create the error:
[bash]
-lmkl_intel_lp64 -lmkl_sequential -lpthread -lm -lmkl_core -L/usr/local/intel/composer_xe_2013.0.079/mkl/lib/intel64 /usr/local/intel/composer_xe_2013.0.079/mkl/lib/intel64/libmkl_lapack95_lp64.a
[/bash]
The next link line seems to work fine, the error specified above does not occur:
[bash]
-lmkl_intel_lp64 -lmkl_intel_thread -liomp5 -lpthread -lm -lmkl_core -L/usr/local/intel/composer_xe_2013.0.079/mkl/lib/intel64 /usr/local/intel/composer_xe_2013.0.079/mkl/lib/intel64/libmkl_lapack95_lp64.a
[/bash]

dan0112

Beginner

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

11-08-2012
12:38 AM

16 Views

Highlighted
##

that's correct :the issue happens only for sequential code (libmkl_sequential.a lib to use). Threaded version is OK.
--Gennady

Gennady_F_Intel

Moderator

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

11-08-2012
01:10 AM

16 Views

Highlighted
##

update 2 just came out, but I can't find this on the fixes list. Is it fixed? I want to know before I push our support to install it on the cluster. Thanks.

dan0112

Beginner

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

02-04-2013
08:14 AM

16 Views

Highlighted
##

Gennady_F_Intel

Moderator

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

02-04-2013
08:34 AM

16 Views

yes, the problem has been fixed in 11.0 update 2. please check and let us know the result.

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