从 C++ 链接到 Fortran 库 (Lapack)

Question

我在我的 C++ 代码中使用Lapack。我很困惑如何正确链接到图书馆。这是一个小例子，对应于我从 Lapack 调用函数的代码：

#include <iostream>

namespace lapack { extern "C" {
  void ilaver(int* major, int* minor, int* patch); } }

int main()
{
    int major = 0;
    int minor = 0;
    int patch = 0;
    lapack::ilaver(&major, &minor, &patch);
    std::cout << major << "." << minor << "." << patch << std::endl;
    return 0;
}

如果我尝试使用 GCC 4.8.5 (Linux openSUSE) 编译它，我会收到以下错误：

> g++ ilaver.cpp -o ilaver -L /softs/lapack/3.7.1/64/gcc/4.8.5/lib64 -l lapack
/tmp/ccHvDCAh.o: In function `main':
ilaver.cpp:(.text+0x33): undefined reference to `ilaver'
collect2: error: ld returned 1 exit status

我知道这是一个名称修改问题。如果我更改我的代码，在函数名称的末尾添加一个下划线，它可以使用 GCC 正确编译：

#include <iostream>

namespace lapack { extern "C" {
  void ilaver_(int* major, int* minor, int* patch); } }

int main()
{
    int major = 0;
    int minor = 0;
    int patch = 0;
    lapack::ilaver_(&major, &minor, &patch);
    std::cout << major << "." << minor << "." << patch << std::endl;
    return 0;
}

但它不能在 Windows 下使用 Intel 的编译器进行编译。在那里，mangling 不同，我必须将其更改为lapack::ILAVER，然后才能编译。

我的代码必须使用多个编译器（GCC、Intel、MSVC）在多个配置（Linux/Mac/Windows）下编译。我怎样才能更通用并确保它可以在大量编译器下编译？

Answer 1

正如您已经发现的那样，添加下划线之类的技巧本质上是特定于平台的。以这种方式“手动”支持大量平台和编译器需要大量令人不快和乏味的工作。

获得可移植性的最简单方法是使用LAPACKE，这是LAPACK 的官方C 接口。这还有一个额外的好处，就是让您不必重新声明您需要的所有功能。

这是一个简单的例子：

#include <iostream>
#include <lapacke.h>    

int main()
{
    // By using lapack_int, we also support LAPACK-ILP64
    lapack_int major = 0;
    lapack_int minor = 0;
    lapack_int patch = 0;
    LAPACKE_ilaver(&major, &minor, &patch);
    std::cout << major << "." << minor << "." << patch << std::endl;
    return 0;
}

更多信息可以在official documentation找到。

请注意，LAPACKE 只处理 LAPACK，如果您还需要 BLAS 例程，可以从 CBLAS 获取。