使用g++链接intel mkl数学库失败

Question

我尝试使用g++编译Intel MKL 11.1。

g++ -m32 test.c -lmkl_intel -lmkl_intel_thread -lmkl_core -liomp5 -lpthread -lm

错误信息如下：

该错误提示：

/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `logf'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `atan2'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `sin'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `fabs'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `exp'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `cos'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `sqrt'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_intel_thread.so: undefined reference to `log'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `pow'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `log10'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `ceil'
/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/ia32/libmkl_core.so: undefined reference to `expf'

但是，如果我使用以下方式编译：

g++ -m32 test.c -lmkl_intel -lmkl_intel_thread -lmkl_core -lm -liomp5 -lpthread -lm

不会出现错误。顺便说一下，我不明白为什么在-lmkl_core后面添加-lm会有帮助。

gcc -m32 test.c -lmkl_intel -lmkl_intel_thread -lmkl_core -liomp5 -lpthread -lm

原来没问题，我也不明白这里有什么区别。
以下是我的测试代码。

/* C source code is found in dgemm_example.c */

#define min(x,y) (((x) < (y)) ? (x) : (y))

#include <stdio.h>
#include <stdlib.h>
#include "mkl.h"

int main()
{
    double *A, *B, *C;
    int m, n, k, i, j;
    double alpha, beta;

    printf ("\n This example computes real matrix C=alpha*A*B+beta*C using \n"
            " Intel(R) MKL function dgemm, where A, B, and  C are matrices and \n"
            " alpha and beta are double precision scalars\n\n");

    m = 2000, k = 200, n = 1000;
    printf (" Initializing data for matrix multiplication C=A*B for matrix \n"
            " A(%ix%i) and matrix B(%ix%i)\n\n", m, k, k, n);
    alpha = 1.0; beta = 0.0;

    printf (" Allocating memory for matrices aligned on 64-byte boundary for better \n"
            " performance \n\n");
    A = (double *)mkl_malloc( m*k*sizeof( double ), 64 );
    B = (double *)mkl_malloc( k*n*sizeof( double ), 64 );
    C = (double *)mkl_malloc( m*n*sizeof( double ), 64 );
    if (A == NULL || B == NULL || C == NULL) {
        printf( "\n ERROR: Can't allocate memory for matrices. Aborting... \n\n");
        mkl_free(A);
        mkl_free(B);
        mkl_free(C);
        return 1;
    }

    printf (" Intializing matrix data \n\n");
    for (i = 0; i < (m*k); i++) {
        A[i] = (double)(i+1);
    }

    for (i = 0; i < (k*n); i++) {
        B[i] = (double)(-i-1);
    }

    for (i = 0; i < (m*n); i++) {
        C[i] = 0.0;
    }

    printf (" Computing matrix product using Intel(R) MKL dgemm function via CBLAS interface \n\n");
    cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, 
            m, n, k, alpha, A, k, B, n, beta, C, n);
    printf ("\n Computations completed.\n\n");

    printf (" Top left corner of matrix A: \n");
    for (i=0; i<min(m,6); i++) {
        for (j=0; j<min(k,6); j++) {
            printf ("%12.0f", A[j+i*k]);
        }
        printf ("\n");
    }

    printf ("\n Top left corner of matrix B: \n");
    for (i=0; i<min(k,6); i++) {
        for (j=0; j<min(n,6); j++) {
            printf ("%12.0f", B[j+i*n]);
        }
        printf ("\n");
    }

    printf ("\n Top left corner of matrix C: \n");
    for (i=0; i<min(m,6); i++) {
        for (j=0; j<min(n,6); j++) {
            printf ("%12.5G", C[j+i*n]);
        }
        printf ("\n");
    }

    printf ("\n Deallocating memory \n\n");
    mkl_free(A);
    mkl_free(B);
    mkl_free(C);

    printf (" Example completed. \n\n");
    return 0;
}

- cloudygoose

---

谢谢，将-lm重复两次放置让我使用MKL构建了MAGMA。 - SkyWalker

我也遇到了同样的问题。你解决了吗？ - étale-cohomology

有关此事有任何更新吗？我也遇到了同样的问题。 - Faser

Answer 1

我正在使用Debian，按照以下步骤成功运行和编译mkl中的一个示例程序：

• 我使用script.sh自动安装必要的软件包（需要使用sudo）

• 按照以下方式运行其中一个示例程序（例如this one）: (ref)

  1. $ export MKLROOT=/opt/intel/compilers_and_libraries_2018.2.199/linux/mkl
  2. $ cc -fopenmp -m64 -I${MKLROOT}/include test.c -Wl,--no-as-needed -L${MKLROOT}/lib/intel64 -lmkl_intel_lp64 -lmkl_core -lmkl_gnu_thread -lpthread -lm -ldl
  3. $ ./a.out

- xxks-kkk