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.

The Stationary Iterations Revisited

Convergence of General Nonstationary Iterative Methods for Solving Singular Linear Equations

Convergence of Newton'S Method and Inverse Function Theorem in Banach Space

Formalization of Asymptotic Convergence for Stationary Iterative Methods

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

On the Convergence of the Accelerated Riccati Iteration Method.

Strong Approximation for A Class of Stationary Processes

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

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

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

Some Stationary Weak Solutions to Inhomogeneous Landau-Lifshitz Equations in Three Dimensions

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

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

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

Convergence of Newton'S Method for Systems of Equations with Constant Rank Derivatives

A New Linear Convergence Result for the Iterative Soft Thresholding Algorithm

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

New strong convergence analysis for variational inequalities and fixed point problems

Convergence Properties Of Krylov Subspace Methods For Singular Linear Systems With Arbitrary Index

On the convergence of splittings for semidefinite linear systems

A Weak Condition for Secant Method to Solve Systems of Nonlinear Equations