链接到C 的Fortran库(Lapack)

Link to Fortran library (Lapack) from C++

本文关键字:Lapack Fortran 链接      更新时间:2023-10-16

我在我的C 代码中使用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

我知道这是一个名称mangling 问题。如果我更改函数名称末尾添加下划线的代码,则使用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下的英特尔编译器一起编译。在那里,乱七八糟的是不同的,我必须将其更改为lapack::ILAVER,然后编译。

我的代码必须在几个配置(Linux/Mac/Windows)的几个编译器(GCC,Intel,MSVC)下编译。我如何更通用,并确保它将在大型编译器面板下进行编译?

正如您已经发现的那样,诸如添加下划线的hacks具有特定于平台的本质。以"手工"方式支持各种平台和编译器需要大量的不愉快和乏味的工作。

获得可移植性的最简单方法是使用拉帕克(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;
}

更多信息可以在官方文档中找到。

请注意,Lapacke仅处理Lapack,如果您还需要Blas例程,则可以从CBLA中获得。