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

MKL_DCOOMM giving incorrect answers

Hi

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.

 

 

0 Kudos
3 Replies
Highlighted
Black Belt
12 Views

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).

0 Kudos
Highlighted
Employee
12 Views

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

0 Kudos
Highlighted
Beginner
12 Views

Many Thanks

 

0 Kudos