从C++打电话给斯卡拉帕克

Calling ScaLAPACK from C++

本文关键字：帕克 C++ 打电话      更新时间：2023-10-16

嗨，任何人都可以为我提供一个示例 Makefile 来从C++调用 ScaLAPACK，我遇到了困难。

我正确编译了最新版本，并通过了所有测试。我在Fedora上用GCC和OpenMPI编译了它，我也尝试使用存储库中预先构建的二进制文件，但没有运气。

我对ScaLAPACK过程使用了以下声明：

extern "C" void pdgesvd_(char *jobu, char *jobvt, int *M, int *N, double *A, int *ia, int *ja, int *desca, double *s, double *u, int *iu, int *ju,
            int *descu, double *vt, int *ivt, int *jvt, int *descvt, double *work, int *lwork, int *info);
extern "C" void pdgemv_(char *t, int *M, int *N, double *alpha, double *A, int *ia, int *ja, int *desca, double *X, int *ix, int *jx,int *descx,
            int *incx, double *beta, double *Y, int *iy, int *jy, int *descy, int *incy);

extern "C" void pdgemm_(char *transa, char *transb, int *M, int *N, int *K, double *alpha, double *A, int *ia, int *ja, int *desca, double *B,
            int *ib, int *jb, int *descb, double *beta, double *C, int *ic, int *jc, int *descc);

BLACS 程序以相同的方式声明：

extern "C" void blacs_get_(/*in*/int *icontxt, /*in*/int *what, /*out*/int *val);

我将"_"字符附加到过程名称中，因为 BLAS 和 ScaLAPACK 是使用 -DAdd_ 编译器选项构建的。对于您的系统，这可能有所不同。如果是这样，请尝试objdump -t libscalapack.a和其他静态库文件以确定实际的过程名称。

带有 ScaLAPACK 调用的程序可以使用如下内容进行编译：

mpicxx -o your_binary_name -O2 your_source.cpp -lscalapack -lblacs -lcblacs -lblacs -llapack -lblas -lgfortran

请注意，g++ 调用行 （ -lblacs -lcblacs -lblacs -llapack -lblas -lgfortran） 中静态库的顺序很重要！还要检查库目录中是否存在所有这些静态库文件（libblas.a、libcblacs.a等）。不要忘记Fortran软件包（-lgfortran）。

相关文章：