如何在 Linux 集群上构建 BLAS 和 LAPACK 以在 C++ 中使用？

Question

我正在处理一个很大的计算问题。 为了降低方阵中一组线性方程的计算速度，我使用了lapack和blas 。 为了在我的笔记本电脑 (Ubuntu 2020) 上获取库，我运行了以下命令

sudo apt-get install libblas-dev liblapack-dev

然后我在编译时通过输入以下内容链接代码

g++ main.cpp -llapack -lblas

但是，我正在处理的集群似乎没有安装这两个库。 它在集群上要慢得多，但仍然是一个更好的芯片。 它运行，所以我认为它安装了lapack库，但没有安装blas 。 我想安装两个。

我将如何构建和编译既不能访问 root 也不能访问apt-get的lapack和blas库？

这是一个用于测试的简短脚本。

#include <iostream>
#include <vector>

extern "C" void dgesv_( int *n, int *nrhs, double  *a, int *lda, int *ipiv, double *b, int *lbd, int *info  );

int main() {
    int SIZE = 3;
    int nrhs = 1; // one column in b
    int lda = SIZE;
    int ldb = SIZE;
    std::vector<int> i_piv(SIZE, 0);  // pivot column vector
    int info;
    std::vector<double> A(SIZE*SIZE, 0); // sq mat with 0's
    A = {5, 2, 8, 9, 7, 2, 10, 3, 4};
    std::vector<double> b(SIZE);
    b = {22, 13, 17};

    dgesv_( &SIZE, &nrhs, &*A.begin(), &lda, &*i_piv.begin(), &*b.begin(), &ldb, &info );
    return 0;
}

我想用

g++ main.cpp -L/path/to/lapack -L/path/to/blas -llapack -lblas

其中 b 矩阵被替换为解，解为1.71, 1.29, 0.18 （这有点随意，所以我没有在代码中提供“print_matrix”函数来减少混乱）。

感谢您的时间。

Answer 1

您可以通过 git 将官方 repo 克隆到您的集群：

git clone https://github.com/xianyi/OpenBLAS.git

现在 cd 到克隆的 repo，编译并安装它：

cd OpenBLAS
make 
make install PREFIX=~/blas

现在您可以通过这种方式链接 BLAS：

g++ main.cpp -L~/blas/lib -lblas

更多关于编译 BLAS 的信息可以在官方 repo 中找到。
LAPACK 也一样。

根据 repo 中的README.md ，可以使用 CMAKE 编译 LAPACK：

cd lapack
mkdir build
cd build
cmake -DCMAKE_INSTALL_LIBDIR=/home/user/lapack ..
cmake --build . -j --target install

Answer 2

BLAS

• 下载最新版本的BLAS

• 打开终端并转到保存它的目录

tar -xvf blas-3.8.0.tgz  # unzip the blas source files
cd BLAS-3.8.0/ 
make
mv blas_LINUX.a libblas.a
mv *.a path/to/lib  # move the blas lib to the library you will be including at compile

---

拉帕克

• 下载最新版本的LAPACK

tar -xvf lapack-3.9.0.tar.gz
cd lapack-3.9.0/
cp make.inc.example make.inc  # use example make as make
make
cp *.a path/to/lib

---

现在已经构建了库并存储在path/to/lib ，可以编译问题中的简短示例代码。

g++ main.cpp -L/path/to/lib -llapack -lblas -lgfortran  # compiles the code
./a.out  # runs the code