我有一个对称的不定矩阵，$H$. 我还有一个可以计算对称不定矩阵的代数最小特征值的例程。我想使用我现有的代码计算最小幅度的特征值。

我可以做到这一点的一种方法是使用我的例程来计算代数上最小的特征对$H^2$. 然后我可以使用瑞利商来计算最小幅度特征值$H$. 这很好用，但速度很慢，因为$cond(H^2) = cond(H)^2$.

有没有更好的办法？