math - 广义奇异值分解和稀疏矩阵

Question

我想计算稀疏矩阵 A 和 B 的广义奇异值分解 (GSVD)。因此，我正在寻找一种能够为稀疏矩阵使用特殊数据结构的实现。

我发现的唯一实现（这里）是用 Fortran 77 编写的 LAPACK 包的一部分。它工作得很好，但不幸的是它不能处理稀疏矩阵。

Answer 1

MATLABgsvd接受稀疏矩阵。我相信 Octave（免费提供）也支持gsvd。

Answer 2

我在 Scicomp 上问了同样的问题并得到了很好的答案。帖子可以在这里找到。