LAPACK与gcc创业指南

Question

我需要帮助建立LAPACK在Linux gcc。我是LAPACK的新手，没有使用Fortran的知识。

我下载了lapack 3.4.0，并使库获得liblapack.a和librefblas.a。

之后，我将这些库链接到我的程序：-llapack -lrefblas

我想使用LAPACK函数，如dpotrf、dgetrf、dgetri等，如何将头文件包括在内，以便编译器找到函数？我有必要使用LAPACK的C接口吗？

Answer 1

在c++中使用LAPACK有几种方法。下面是我要做的事情(假设您使用的是某个*nix系统)。

1. 确保您拥有正确的库，并且知道正确的编译器/链接器选项集。由于这些代码是用Fortran编写的，所以我从Fortran代码开始。就像this one。确保您可以使用gfortran编译它。可能的链接器选项可以是(取决于您的系统)：-llapack、-lblas或其中的一些组合。然后，
2. 使用C++接口。同样，有多种方法可供选择。我个人认为使用clapack接口是最简单的，您在其中声明的LAPACk函数如下：

extern "C" void dsyev_( char *jobz, char *uplo, int &n, double *a, int &lda, double *w, double *work, int &lwork, int &info);

正确的链接器选项集再次取决于您的系统，并且可以是这样的：-llapack -lf77blas -latlas (这个集合在我的Ubuntu上工作，其中LAPACK来自atlas包)。