1

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

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

4

1 回答 1

3

我对 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 包 (for -lgfortran)。

于 2012-09-26T17:00:33.417 回答