广义奇异值分解＆amp;稀疏矩阵

Question

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

我找到的唯一实现（here）是用Fortran 77编写的LAPACK包的一部分。 它工作得很好，但不幸的是它无法处理稀疏矩阵。

Answer 1

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

Answer 2

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