我正在玩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。