我想计算稀疏矩阵A和B的广义奇异值分解(GSVD)。因此,我正在寻找一种能够使用稀疏矩阵的特殊数据结构的实现。我找到的唯一实现(here)是LAPACK包的一部分,它是用Fortran 77编写的。它工作得很好,但不幸的是它无法处理稀疏矩阵。