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

Question

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

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

Answer 1

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