我正在使用带有 Trilinos 的 Albany 来解决热膨胀不匹配的弹性问题。我正在使用带有 MueLu 预处理的块 GMRES。它适用于数百万自由度的问题大小，但是当问题大小变大时，即使预处理残差收敛，真实残差也不会收敛。我为 MueLu 尝试了不同的参数（sym Gauss-Seidel，Jacobi，...），但它不起作用。关于收敛问题的任何建议？

对于使用 MueLu 与 Block GMRES 收敛的问题大小，我尝试了使用 ILU 预处理器的 Block GMRES，并且真正的残差范数保持为 1。这是否意味着我有一个带有 ILU 的奇异预处理器？为什么会这样？

谢谢。