The main focus of this paper is the study of efficient multigrid methods for large linear systems with a particular saddle-point structure. Indeed, when the system matrix is symmetric, but indefinite, the variational convergence theory that is usually used to prove multigrid convergence cannot be directly applied. However, different algebraic approaches analyze properly preconditioned saddle-point problems, proving convergence of the Two-Grid method. In particular, this is efficient when the blocks of the coefficient matrix possess a Toeplitz or circulant structure. Indeed, it is possible to derive sufficient conditions for convergence and provide optimal parameters for the preconditioning of the saddle-point problem in terms of the associated generating symbols. In this paper, we propose a symbol-based convergence analysis for problems that have a hidden block Toeplitz structure. Then, they can be investigated focusing on the properties of the associated generating function f, which consequently is a matrix-valued function with dimension depending on the block size of the problem. As numerical tests we focus on the matrix sequence stemming from the finite element approximation of the Stokes problem. We show the efficiency of the methods studying the hidden 9-by-9 block multilevel structure of the obtained matrix sequence. Moreover, we propose an efficient algebraic multigrid method with convergence rate independent of the matrix size. Finally, we present several numerical tests comparing the results with state-of-the-art strategies.
翻译:本文主要研究具有特定鞍点结构的大型线性系统的高效多重网格方法。当系统矩阵对称但不定时,通常用于证明多重网格收敛性的变分收敛理论无法直接应用。然而,不同的代数方法通过恰当预处理鞍点问题,证明了二重网格方法的收敛性。特别地,当系数矩阵的块具有Toeplitz或循环结构时,该方法是高效的。实际上,可以推导出收敛的充分条件,并根据关联生成符号提供鞍点问题预处理的最优参数。本文提出了一种基于符号的收敛性分析方法,适用于具有隐式块Toeplitz结构的问题。然后,可通过研究关联生成函数f的性质来分析这些问题,该函数是依赖于问题块大小的矩阵值函数。在数值测试中,我们重点研究了由Stokes问题有限元近似生成的矩阵序列。通过分析所得矩阵序列的隐式9×9块多层结构,展示了该方法的效率。此外,我们提出了一种收敛速率与矩阵规模无关的高效代数多重网格方法。最后,我们通过多项数值测试,将结果与现有先进策略进行了比较。