Structured Backward Error Analysis with Preserving Sparsity for Three-by-Three Block Saddle Point Problems

Sk. Safique Ahmad,Pinki Khatun
2024-08-26
Abstract:Backward error (\textit{BE}) analysis emerges as a powerful tool for assessing the backward stability and strong backward stability of iterative algorithms. In this paper, we explore structured \textit{BEs} for a class of three-by-three block saddle point problems (\textit{TBSPPs}), aiming to assess the strong backward stability of iterative algorithms devised to find their solution. Our investigations preserve the inherent matrix structure and sparsity pattern in the corresponding perturbation matrices and derive explicit formulae for the structure \textit{BEs}. Moreover, we provide formulae for the structure-preserving minimal perturbation matrices for which the structured \textit{BE} is attained. Furthermore, using the obtained structured \textit{BEs}, we propose three new termination criteria for the iterative algorithms to solve the \textit{TBSPP.} Numerical experiments are performed to test the strong backward stability of various iterative algorithms, as well as the reliability and effectiveness of the proposed termination criteria.
Numerical Analysis
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the structured backward error analysis problem of the three - by - three block saddle point problems (TBSPP) in order to evaluate the strong backward stability of iterative algorithms. Specifically, the authors focus on: 1. **Structured Backward Error (SBE)**: - Research how to calculate the structured backward error of TBSPP while maintaining the inherent structure and sparse pattern of the matrix. - Derive specific formulas to calculate these structured backward errors. 2. **Minimal Perturbation Matrices**: - Provide explicit formulas for the minimal perturbation matrices that maintain the original matrix structure and sparse pattern, making the approximate solution the actual solution of the slightly perturbed TBSPP. 3. **New Termination Criteria**: - Utilize the obtained structured backward error to propose three new termination criteria for iterative algorithms, in order to improve the effectiveness and reliability of solving TBSPP. 4. **Numerical Experiments**: - Verify the proposed theoretical results through numerical experiments, test the backward stability and strong backward stability of different iterative algorithms, and evaluate the effectiveness of the newly proposed termination criteria. ### Background and Motivation The three - by - three block saddle point problem (TBSPP) is an important class of linear systems and has wide applications in the field of scientific computing, such as solving PDE - constrained optimization problems, liquid crystal orientation models, Stokes - Darcy equations, etc. In order to ensure the accuracy and reliability of iterative algorithms, backward error analysis (BEA) is a powerful tool. However, most of the existing research focuses on two - by - two block saddle point problems and does not fully consider the sparsity and structure preservation of matrices. ### Main Contributions 1. **Structured Backward Error Analysis**: - Under the premise of maintaining the structure and sparse pattern of the coefficient matrix A, the structured backward error in three different situations is studied. 2. **Minimal Perturbation Matrix Formulas**: - Provide explicit formulas for the minimal perturbation matrices that maintain the original matrix structure and sparse pattern. 3. **New Termination Criteria**: - Propose three new termination criteria for iterative algorithms, which improve the effectiveness and reliability of solving TBSPP. 4. **Numerical Experiments**: - Verify the theoretical results through numerical experiments and show the effectiveness and superiority of the new method. ### Conclusion By introducing structured backward error analysis and maintaining the sparsity and structure of matrices, this paper provides a more accurate and reliable method for solving TBSPP, and also provides new ideas and tools for research in related fields.