有大量关于鞍点问题的有效预处理器的文献。在计算物理学中，几乎总是考虑拉格朗日乘数具有弱公式（例如斯托克斯方程）的情况。在这种情况下，可以设计非常快速的求解器。

不幸的是，在结构力学中，拉格朗日乘数经常被引入离散水平。例如考虑自由度之间的多点约束（此功能在工业软件中非常常见）。我知道这不是做事的好方法，但这就是通常的做法。

我知道 St-Georges、Notay、Warzée 的古代作品，或者 Rees 和 Greif、Gould、Hribar 和 Nocedal 的近期作品。但我想知道是否有人经历过其他策略（例如 PETSc 的 FieldSplit 预处理器中可用的策略）用那些可怕的离散拉格朗日乘数对鞍点系统进行预处理？我正在寻找一种在这些鞍点系统上使用高效预处理器（比如多重网格）的方法。

最后，我想知道使用具有拉格朗日乘数的高效预处理器的最佳方法是否不是可以自然处理乘数的域分解方法，例如 FETI 或 BDD。