嗨,任何人都可以为我提供一个示例 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
)。