Abstract:Multilevel methods represent a powerful approach in numerical solution of partial differential equations. The multilevel structure can also be used to construct estimates for total and algebraic errors of computed approximations. This paper deals with residual-based error estimates that are based on properties of quasi-interpolation operators, stable-splittings, or frames. We focus on the settings where the system matrix on the coarsest level is still large and the associated terms in the estimates can only be approximated. We show that the way in which the error term associated with the coarsest level is approximated is substantial. It can significantly affect both the efficiency (accuracy) of the overall error estimates and their robustness with respect to the size of the coarsest problem. The newly proposed approximation of the coarsest-level term is based on using the conjugate gradient method with an appropriate stopping criterion. We prove that the resulting estimates are efficient and robust with respect to the size of the coarsest-level problem. Numerical experiments illustrate the theoretical findings.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to effectively estimate and approximate the total approximate error and algebraic error in multi - level decomposition methods when the problem size at the coarsest level is still very large. Specifically, the article focuses on posteriori error estimates based on residuals, which depend on the properties of quasi - interpolation operators, stable splittings or frames. When the system matrix at the coarsest level is so large that the relevant estimate terms can only be approximated, the effectiveness and robustness of error estimates are significantly affected in this case.
### Main problems of the paper
1. **Impact of large - scale coarsest - level problems**: In multi - level methods, the problem at the coarsest level is usually small and can be directly solved. However, when dealing with complex domains or large - scale parallel computing, the problem at the coarsest level may still be large and cannot be directly solved. Therefore, effective methods need to be found to approximate the solutions of these large - scale problems and ensure the accuracy and robustness of error estimates.
2. **Effectiveness and robustness of error estimates**: Traditional error estimation methods assume that the problem at the coarsest level can be solved exactly, but in practical applications, this is often impossible. Therefore, new methods need to be developed to approximate the error terms at the coarsest level and prove the effectiveness and robustness of these approximation methods.
3. **Application of the conjugate gradient method**: The paper proposes a new method based on the conjugate gradient method (CG) to approximate the error terms at the coarsest level and ensure its efficiency and robustness through appropriate stopping criteria. This method is not only suitable for large - scale problems but also can guarantee the quality of error estimates.
### Formula representation
The formulas involved in the article include:
- Definition of the residual \( r \):
\[
\langle r, w \rangle=\int_{\Omega} f w-\nabla v_J\cdot\nabla w\quad\forall w\in H^1_0(\Omega)
\]
- Residual equation:
\[
\|\nabla e\|^2 = \langle r, e \rangle
\]
- Riesz representation at the coarsest level:
\[
\langle r, w_0 \rangle=\int_{\Omega}\nabla r_0\cdot\nabla w_0\quad\forall w_0\in V_0
\]
- Multi - level error estimate:
\[
\|\nabla e\|\leq C_{cls}\eta_J + C_{I,2lvl}\sum_{j = 1}^J\|h^{-1}_j r_j\|+C_{IV0,4}\|\nabla r_0\|
\]
### Summary
The core problem of this paper is to solve the error estimation problem when the problem size at the coarsest level in multi - level decomposition methods is large. By introducing a new approximation method based on the conjugate gradient method, the author has successfully improved the effectiveness and robustness of error estimates, which is of great significance for complex domains and large - scale parallel computing.