Problems with using MKL sgetrf

[taniab]

I am compiling a big code base for contour mapping that uses sgetrf and sgetri. The code compiles fine, but I get two error messages as given below, followed by a segmentation fault:
Parameter 1 to routine SGETRF was incorrect
Parameter 1 to routine SGETRI was incorrect

To recreate the error on a smaller test, I have created the attached file test.c.

1) I compiled it using:

g++ -g test.c -lmkl_intel_ilp64 -lmkl_sequential -lmkl_core -lmkl_blacs_intelmpi_ilp64 -lpthread -lm -ldl -lmpi

The original code uses nvcc.

2) MKLROOT is set by setvars.sh to be:
/opt/intel/oneapi/mkl/latest

ls /opt/intel/oneapi/mkl/
2021.2.0/ latest/

Thank you for any help with this.

[mecej4 · Accepted solution]

If the integer arguments being passed to MKL Lapack and BLAS routines are 32-bit integers, you should not be using the ILP64 libraries. You should use the LP64 libraries, instead. Use the MKL Link Line Advisor to obtain the correct command line.

[taniab]

Thank you so much. Using LP64 libraries fixed the problem.

[MRajesh_intel]

 Hi,

Thanks for the confirmation!

As this issue has been resolved, we will no longer respond to this thread. If you require any additional assistance from Intel, please start a new thread. Any further interaction in this thread will be considered community only.

Have a Good day.

Regards

Rajesh