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

Question

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

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

Answer 1

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

Answer 2

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