c++ - 从 C++ 调用 ScaLAPACK

Question

您好，谁能给我提供一个从 C++ 调用 ScaLAPACK 的 Makefile 示例，我遇到了困难。

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

Answer 1

我对 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_ 编译器选项构建的。这可能与您的系统不同。如果是这样，请在 libscalapack.a 和其他静态库文件上尝试 objdump -t 以确定实际的过程名称。

您的带有 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)。