一种可能性是使用 ARPACK 和 ViennaCL 的组合：

• ARPACK 是一个特征求解器。它与回调接口一起使用（您提供一个计算$Ax$对于给定的$x$它通过多次调用此函数来计算特征值）。如果要解决广义特征值问题（$Ax = \lambda Bx$) 你需要为计算提供回调$x \rightarrow Bx$也是。
• 这些回调函数可以在 GPU 中实现（它将与 ARPACK 一起使用），现在实现它的方式将取决于矩阵的结构$A$和$B$. 如果您的矩阵是满的、对角的或稀疏的，您将使用不同的实现。ViennaCL 为不同的矩阵结构提出了几个函数。

http://www.caam.rice.edu/software/ARPACK/

http://viennacl.sourceforge.net/

除了 ARPACK，同样的方法可以应用于任何具有反向通信 API（即，具有计算矩阵向量积的回调）的特征求解器，例如 SLEPc（参见下面 Daniel Shapero 的评论）。