Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Highlighted
Beginner
102 Views

Potential bug in complex versions of mkl_?omatadd

Hi all,

I think came upon a bug in the mkl_?omatadd routine, especially the complex versions. If both matrices A and B are conjugate transposed, the matrix A is only conjugated. I configured the zomatadd.c example file that intel provided in order to show that this is the case:

Screenshot 2020-10-20 at 13.26.30.png

Screenshot 2020-10-20 at 13.23.00.png

As you can see the results are the same, which is incorrect. The second computation gives the wrong answer.

There are actually two more cases I found, where the result is wrong:

  1. A transposed and B conjugated: This combination conjugate transposes A instead of just transposing it.
  2. A conjugated and B transposed: This combination conjugate transposes B instead of just transposing it.

The example file is also attached for you to test it.

 

Best regards,

Simon

0 Kudos
4 Replies
Highlighted
Moderator
94 Views

what version of mkl do you use?


0 Kudos
Highlighted
Beginner
87 Views

I'm using the mkl that comes with 2020 update 2 parallel studio. I'm on MacOS. 

0 Kudos
Highlighted
Moderator
66 Views

Yes, you are right, that's the issue. You may escalate the case to Intel Online Service Center against Intel MKL.


0 Kudos
Highlighted
Beginner
57 Views

Thanks for confirming the issue. There was already a ticket created by Khoros, is that enough? 

0 Kudos