Abstract:Solutions to differential equations, which are used to model physical systems, are computed numerically by solving a set of discretized equations. This set of discretized equations is reduced to a large linear system, whose solution is typically found using an iterative solver. We start with an initial guess, $x_0$, and iterate the algorithm to obtain a sequence of solution vectors, $x_k$, which are approximations to the exact solution of the linear system, $x$. The iterative algorithm is said to converge to $x$, in the field of reals, if and only if $x_k$ converges to $x$ in the limit of $k \to \infty$. In this paper, we formally prove the asymptotic convergence of a particular class of iterative methods called the stationary iterative methods, in the Coq theorem prover. We formalize the necessary and sufficient conditions required for the iterative convergence, and extend this result to two classical iterative methods: the Gauss--Seidel method and the Jacobi method. For the Gauss--Seidel method, we also formalize a set of easily testable conditions for iterative convergence, called the Reich theorem, for a particular matrix structure, and apply this on a model problem of the one-dimensional heat equation. We also apply the main theorem of iterative convergence to prove convergence of the Jacobi method on the model problem.

On the Convergence of General Stationary Iterative Methods for Range‐hermitian Singular Linear Systems

Convergence of General Nonstationary Iterative Methods for Solving Singular Linear Equations

Convergence and Quotient Convergence of Iterative Methods for Solving Singular Linear Equations with Index One

On the Convergence of the Accelerated Riccati Iteration Method.

Convergence Criterion of Newton's Method for Singular Systems with Constant Rank Derivatives

The Stationary Iterations Revisited

Newton-type Methods for Inverse Singular Value Problems with Multiple Singular Values

Extended Newton-type Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

On the Hermitian positive definite solution and Newton's method for a nonlinear matrix equation

The Convergence of Least-Squares Progressive Iterative Approximation for Singular Least-Squares Fitting System.

Quotient Convergence and Multi-Splitting Methods for Solving Singular Linear Equations

Convergence Behavior of Gauss-Newton's Method and Extensions of the Smale Point Estimate Theory.

Formalization of Asymptotic Convergence for Stationary Iterative Methods

On the Local Convergence of an Iterative Approach for Inverse Singular Value Problems

Computing the Hermitian Positive Definite Solutions of a Nonlinear Matrix Equation

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Convergence of the Newton-type methods for the square inverse singular value problems with multiple and zero singular values

Stationary splitting iterative methods for the matrix equation A X B = C

On Semi-Convergence of the Generalized Shift-Splitting Iteration Method for Singular Nonsymmetric Saddle Point Problems

Convergence Properties of an Iterative Method for Solving Symmetric Non-linear Equations