[Trilinos-Users] Trilinos and MKL

[Ramsey, James J CIV (US)]

________________________________________________________

From: trilinos-users-bounces at software.sandia.gov [trilinos-users-bounces at software.sandia.gov] on behalf of Villa Andrea (RSE) [Andrea.Villa at rse-web.it]
Sent: Tuesday, August 20, 2013 4:26 AM
To: trilinos-users at software.sandia.gov
Subject: [Trilinos-Users] Trilinos and MKL

Dear Trilinos developers, I'm experiencing some problems using MKL with trilinos. I've used the following commands in the configuration script:


-D BLAS_INCLUDE_DIRS:PATH=/opt/intel/Compiler/cc/composer_xe_2011_sp1.9.293/mkl/include/ \
-D BLAS_LIBRARY_DIRS:PATH=/opt/intel/Compiler/cc/composer_xe_2011_sp1.9.293/mkl/lib/intel64/ \
-D BLAS_LIBRARY_NAMES:STRING="mkl_intel_ilp64; mkl_sequential; mkl_core" \
-D LAPACK_INCLUDE_DIRS:PATH=/opt/intel/Compiler/cc/composer_xe_2011_sp1.9.293/mkl/include/ \
-D LAPACK_LIBRARY_DIRS:PATH=/opt/intel/Compiler/cc/composer_xe_2011_sp1.9.293/mkl/lib/intel64/ \
-D LAPACK_LIBRARY_NAMES:STRING="mkl_intel_ilp64" \

However when I run some simple test codes involving epetra I get the error message attached below. Does anyone has any suggestion to fix the problem? The compiler is gnu.

Thank you very much.

Andrea.


MKL ERROR: Parameter 4 was incorrect on entry to DGETRF
MKL ERROR: Parameter 3 was incorrect on entry to DGETRI
________________________________________________________

I noticed that in the library names for BLAS and LAPACK, you have "mkl_intel_ilp64" rather than "mkl_intel_lp64". See what happens if you use "mkl_intel_lp64" instead. FYI, some documentation on the difference between ILP64 and LP64 is here: http://software.intel.com/sites/products/documentation/hpc/mkl/mkl_userguide_lnx/GUID-87821148-338B-4022-8C90-F24C866F2878.htm

My guess is that the error messages from MKL reflect an issue with integer sizes.