例如,不定矩阵系统出现在通过混合有限元对鞍点问题进行离散化中。然后可以将系统矩阵置于以下形式
其中是负(半)定,是正(半)定,是任意的。当然,根据惯例,您可以使用确定性条件,但这几乎就是这些矩阵的结构。
对于这些方法,可以采用Uzawa的方法,这实际上只是将系统转化为可以通过共轭梯度、梯度下降等方法求解的等价半定系统的“技巧”。
我面临一个没有这种块结构的不确定系统。在这种情况下,宇泽类型的方法不适用。我知道 Paige & Saunders 引入的最小残差方法 (MINRES),它只是一个三项递归,似乎很容易实现。
问题: MINRES 通常是一个不错的选择,例如,用于原型设计吗?有实际意义吗?目前,预处理不是核心问题。