Math 广义奇异值分解&；稀疏矩阵

我想计算稀疏矩阵的广义奇异值分解（GSVD）

我发现的唯一实现（）是用Fortran 77编写的LAPACK包的一部分。

---

它工作得很好，但不幸的是它不能处理稀疏矩阵。

MATLAB的接受。我相信Octave（免费提供）也支持gsvd。

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

为什么不发布此帖子？”包含四种数值（迭代）方法，用于使用双精度ANSI Fortran-77计算大型稀疏矩阵的奇异值分解（SVD）。“可能会有所帮助。我没有这个软件包的经验。谢谢，但我正在寻找广义奇异值分解。可能不是你想要的答案，但是你考虑过写你自己的吗？这是我第一次尝试做的。但我仍然不理解算法，我的主要问题。接受稀疏矩阵只是意味着将它们转换为完全存储并计算所有奇异值。这里的问题是使用与

svds

相同语义的广义SV分解，它只返回几个奇异值，而不是全部。