Algorithms for data assimilation try to predict the most likely state of a dynamical system by combining information from observations and prior models. Variational approaches, such as the weak-constraint four-dimensional variational data assimilation formulation considered in this problem, can ultimately be interpreted as a minimization problem. One of the main challenges of such a formulation is the solution of large linear systems of equations which arise within the inner linear step of the adopted nonlinear solver. Depending on the adopted approach, these linear algebraic problems amount to either a saddle point linear system or a symmetric positive definite (SPD) one. Both formulations can be solved by means of a Krylov method, like GMRES or CG, that needs to be preconditioned to ensure fast convergence in terms of the number of iterations. In this paper we illustrate novel, efficient preconditioning operators which involve the solution of certain Stein matrix equations. In addition to achieving better computational performance, the latter machinery allows us to derive tighter bounds for the eigenvalue distribution of the preconditioned linear system for certain problem settings. A panel of diverse numerical results displays the effectiveness of the proposed methodology compared to current state-of-the-art approaches.
翻译:数据同化的解算法试图通过综合来自观测和先前模型的信息来预测动态系统的最可能状态。不同方法,例如这一问题中考虑的弱的限制四维变异数据同化配方,最终可以被解释为一个最小化问题。这种配方的主要挑战之一是在采用的非线性求解器的内线步骤内产生的大型线性方程系统的解决办法。根据采用的方法,这些线性代数问题相当于一个支撑点线性线性系统或一个对称正数性确定(SPD)1。两种配方都可以通过Krylov方法(如GMRES或CG)来解决,这需要有一个先决条件,以确保迭代数的快速趋同。在本文中,我们举例说明了涉及某些Stein矩阵方程式解决办法的新型、高效的前提条件操作者。除了实现更好的计算性能外,后一种机制还使我们能够为某些问题设置的前提线性线性系统(SPD)的乙基值分布得出更紧的界限。一个不同数字方法小组比较了当前方法的有效性。