我有一个问题，我需要找到一个小的（通常小于 60x60）非对称矩阵的所有正（如特征值为正）特征对。当特征值小于某个阈值时，我可以停止计算。我知道特征值是真实的。关于我可以用来尝试挤出最佳性能的算法的任何建议？我必须做几千次这样的分解，所以速度很重要。

先感谢您。

编辑：我需要在共享内存中的 GPU 上执行此操作。矩阵也不一定大小相同。我目前不知道有任何图书馆这样做。非常适合该问题的算法建议将不胜感激。