topic Call MKL Lapack subroutines and run on HPC machine in Intel® oneAPI Math Kernel Library

Call MKL Lapack subroutines and run on HPC machine

diandian — Thu, 04 Oct 2012 16:13:12 GMT

Hi,

I have problems to call the fortran Lapack subroutines. I am running my fortran 95 code on the HPC machine at the University of Arizona. I need to call SYTRI and SYTRF. The sample codes attached in the machine do not work. Could you let me know what is the correct statement in fortran code to call those subroutines and what is the compiling command?

What I am doing is the following:

!module load intel
ifort -O prgram.f90

!I call the subroutines in my code as (this is the same as the sample codes in the HPC documentation):

USE mkl95_PRECISION, ONLY: WP => SP
USE mkl95_LAPACK, ONLY: SYTRI, SYTRF

Implicit none

! then
CALL sytrf(MATA, 'U', IPIV, INFO)
CALL SYTRI(MATA, IPIV, 'U', INFO)

! The errors are:

matrix_inverse.f90(12): error #6285: There is no matching specific subroutine for this generic subroutine call. [SYTRF]
CALL sytrf(MATA, 'U', IPIV, INFO)
---------^
matrix_inverse.f90(13): error #6285: There is no matching specific subroutine for this generic subroutine call. [SYTRI]
CALL SYTRI(MATA, IPIV, 'U', INFO)
---------^
compilation aborted for matrix_inverse.f90 (code 1)

Thanks for your help

There's an example int the

TimP — Thu, 04 Oct 2012 19:53:28 GMT

There's an example int the mkl examples lapack95 folder which comes with MKL. The message appears to indicate that the data types and array ranks don't match one of the supported combinations (IPIV a rank 1 integer array, INFO a rank 0 integer, MATA a rank 2 real or complex array). This would happen if you omitted any of those declarations, if implicit none didn't complain first.

Thanks for your comments. I

diandian — Thu, 04 Oct 2012 22:22:48 GMT

Hi, Sir: Thanks for your comments. I think I define IPIV(3) , INFO as integer and MATA is a 3X3 symetric real matrix. Now the problem becomes: matrix_inverse.f90:(.text+0xfc): undefined reference to `dsytrf_mkl95_' matrix_inverse.f90:(.text+0x1bf): undefined reference to `dsytri_mkl95_' So, there must be a problem in the US statement: USE mkl95_LAPACK, ONLY: SYTRI, SYTRF But this is the code I get from the documentation within the HPC sever.

You have to specify to the

mecej4 — Fri, 05 Oct 2012 10:26:44 GMT

You have to specify to the linker that you want the Lapack-95 library to be searched. The "USE mkl95_LAPACK" directive instructs the compiler to look up the correct interfaces to the xxx_mkl95_ routines, but the libraries containing the object codes of those routines must be specified by the user.