topic Maybe I am not making myself in Intel® oneAPI Math Kernel Library

MKL 2017 Update 1, MKL_DIRECT_CALL and zgetri issue

AndrewC — Tue, 21 Feb 2017 18:18:23 GMT

I have found that when MKL_DIRECT_CALL is on, there seems to be an issue with one of the zgetrf/zgetri and 5x5 matrices. I am getting quite different results compared to MKL 2016 ( which did not have LAPACK direct call functionality) and I have tracked it down to the DIRECT_CALL implementation of theses LAPACK functions. Is this a known issue?

It would be really nice if the DIRECT _CALL allowed , via the use of #DEFINEs to be able to turn off the functionality in a more fine grained manner rather than all or nothing.

Just as a hint. I am pretty

AndrewC — Tue, 21 Feb 2017 19:07:08 GMT

Just as a hint. I am pretty sure the problem is not in 4x4 or 6x6 matrices as I would have had 100's of QA failures. It **seems** to be related to 5x5 which is an unusual size in my QA cases so only a couple failed. Maybe someone could look at this and see something by inspection.

Here is a reproducer. Just

AndrewC — Tue, 21 Feb 2017 20:10:00 GMT

Here is a reproducer. Just comment out MKL_DIRECT_CALL and the results change significantly

Tested on WIndows, Visual Studio 2015, MKL 2017 Update 1, 64-bit

Test project MUST be using the Intel Compiler!

Dear customer,

Zhen_Z_Intel — Wed, 22 Feb 2017 03:31:00 GMT

Dear customer,

I could not reproduce your problem in same platform(VS 2015& MKL2017u1&Intel64), could you please provide more info:

1. which version of MKL you compared? you said MKL2016, could you please provide specific version number.
2. what's the result you get? please paste here.
3. how's your compiling & linking commands? please specify your linked libraries.

I am looking forward your reply. Thanks.

Best regards,
Fiona

Maybe I am not making myself

AndrewC — Wed, 22 Feb 2017 06:03:32 GMT

Maybe I am not making myself clear.

I am using Intel C++ Composer XE 2017 Update 1 with MKL 2017 Update 1
Comment the line #define MKL_DIRECT_CALL and build/run the test, saving the output
Un-comment the line , so that MKL_DIRECT_CALL is defined and build/run, saving the output.
The output is different. The correct result is when MKL_DIRECT_CALL is NOT #defined/used

.I have attached the complete Visual Studio Project.

Dear customer,

Zhen_Z_Intel — Wed, 22 Feb 2017 07:17:00 GMT

Dear customer,

I found the problem, thank you for posting the issue. We will investigate on it.

Best regards,
Fiona

I have opened a premier

AndrewC — Wed, 22 Feb 2017 20:32:05 GMT

I have opened a premier support ticket #02634419 on this.

Hi vasci_, Regards,

Eugene_C_Intel1 — Thu, 23 Feb 2017 00:13:52 GMT

vasci_,

We investigated the issue. It is caused by a bug in implementation of LU (?GETRF routine) for MKL_DIRECT_CALL for size 5x5. It will be fixed in the nearest MKL release. But as a workaround you can fix it in your own mkl_direct_lapack.h header. Please find attached a patch with the fix.

Thank you for your report!

Regards,
Eugene

OK, thanks! Hopefully this

AndrewC — Thu, 23 Feb 2017 01:15:41 GMT

OK, thanks! Hopefully this will make it into Update 2.

Hi Andrew, We have added the

Gennady_F_Intel — Fri, 12 May 2017 04:46:05 GMT

Hi Andrew, We have added the fix in the latest MKL version ( MKL v.2017 update 3) which has been released yesterday. Please update, take a look and make sure that we have indeed fixed the problem. Thank you!

OK, I am going through the

AndrewC — Mon, 15 May 2017 18:16:46 GMT

OK, I am going through the process of updating to "Composer Update 4"

The issue is fixed!

AndrewC — Mon, 15 May 2017 18:43:29 GMT

The issue is fixed!