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

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- LAPACKE_dtpqrt and LAPACKE_dtpmqrt bugs report

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

hussam_a_

Beginner

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

08-23-2019
09:03 AM

78 Views

LAPACKE_dtpqrt and LAPACKE_dtpmqrt bugs report

machine: MAC

mkl version: 2019.4.233

compiler: clang-1001.0.46.4

LAPACKE_dtpqrt:

layout row major

prameters m = n = l = lda = ldb = ldt = 12

if we set nb > 1,

the routines generates nan in the lower part of the triangular matrices stored in t.

This makes LAPACKE_dtpmqrt returns an error indicating invalid t array.

LAPACKE_dtpmqrt:

layout row major, left side multiplication, no transposition

case 1:

prameters m = k = n = l = ldv = ldt = lda = ldb = 12

If nb != k

The routines does not compute the correct product.

case 2:

prameters m = k = l = nb = ldv = ldt = 12,

n = lda = ldb

if n != k

the routines returns an error -14 meaning that lda does not have a valid value,

however, it should be valid

Link Copied

6 Replies

Gennady_F_Intel

Moderator

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

08-26-2019
12:28 AM

78 Views

hussam_a_

Beginner

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

08-26-2019
04:43 AM

78 Views

Here I attach the file.

I compute the QR decomposition of the matrix [a ; b] (vertical block concatenation) by using dtpqrt.

I set a matrix [c ; d] to match the first n columns of the identity and apply Q from the left by using dtpmqrt.

I get the error from dtpmqrt.

I form the concatenated matrices from their original data (i.e., before any computation).

I call dgeqrf and dormqr to compare data and clearly it is wrong

Note: In this example I didn't get nans from dtpqrt in my other code I did have some.

Gennady_F_Intel

Moderator

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

08-26-2019
08:49 PM

78 Views

Gennady_F_Intel

Moderator

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

09-11-2019
07:40 PM

78 Views

Gennady_F_Intel

Moderator

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

09-11-2019
07:41 PM

78 Views

pls look at the mkl reference : https://software.intel.com/en-us/mkl-developer-reference-c-tpmqrt

hussam_a_

Beginner

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

09-17-2019
01:55 AM

78 Views

Following the mkl reference the argument lda is the leading dimension of the array c in the code which is equal to n in that case (number of columns). I do pass this argument.

Excuse me I did not get what is wrong.

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

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