Elemental的分布式SVD在什么情况下不起作用?

计算科学 线性代数 C++ mpi 拉帕克 svd
2021-12-23 14:17:58

我正在玩Elemental 的分布式奇异值分解,并遇到了两个特殊问题。

tests/lapack_like/SVD.cpp构建测试,并运行

mpirun -np 4 ./SVD --height 3 --width 3

我明白了

SVD residual was unacceptably large

对于分布式Complex<float>测试(但顺序是可以的)。

此外,当我运行时--height 50 --width 52,我遇到了一些问题

H and A must have the same height

SVD 的 Elemental 文档已经过时,我不知道对矩阵的维度有什么要求才能采用其 SVD。

我用拉取的 OpenBLAS、Scalapack、Metis 和 Parmetis 构建了 Elemental。即使使用参考 BLAS/LAPACK 实现,我也无法连接 libFLAME。

0个回答
没有发现任何回复~