MKL_DCOOMM giving incorrect answers

MKL_DCOOMM giving incorrect answers


composer_xe_2011_sp1.7.256  and MKL_DCOOMM

In a test program attached as a file MKL_DCOOMM and DGEMM give different answers. The answer given by DGEMM is correct. I would be greatly obliged if anybody can guide me in finding the reason behind this inconsistency.



Downloadapplication/octet-stream mkl-dcoomm.f901.47 KB
4 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

There is a problem with passing INTEGER*8 arrays to library routines that expect INTEGER*4 arrays. No such problem occurs with scalar integer arguments as long as the variables do not need more than 32 bits to represent them.

Change INTEGER*4 to INTEGER in your program, and things will work correctly. If you really need to use INTEGER*8 arrays (for some reason that you have not so far explained), you will have to use a version of MKL that supports them (e.g., for X64 targets using the ILP64 libraries).

Either change INTEGER*8 to INTEGER in your code, or turn on the ILP64 support in MKL. Your program will work fine.

The easiest way to do (in this particular case) is to pass '-i8' option to ifort. For example, ifort -i8 -mkl mkl-dcoomm.f90

Many Thanks


Leave a Comment

Please sign in to add a comment. Not a member? Join today