Score: 0

Stabilization of the Gradient Method for Solving Linear Algebraic Systems -- A Method Related to the Normal Equation

Published: May 31, 2025 | arXiv ID: 2506.00702v1

By: Ibrahima Dione

Potential Business Impact:

Makes computer math problems solve much faster.

Business Areas:
A/B Testing Data and Analytics

Although it is relatively easy to apply, the gradient method often displays a disappointingly slow rate of convergence. Its convergence is specially based on the structure of the matrix of the algebraic linear system, and on the choice of the stepsize defining the new iteration. We propose here a simple and robust stabilization of the gradient method, which no longer assumes a structure on the matrix (neither symmetric, nor positive definite) to converge, and which no longer requires an approach on the choice of the stepsize. We establish the global convergence of the proposed stabilized algorithm under the only assumption of nonsingular matrix. Several numerical examples illustrating its performances are presented, where we have tested small and large scale linear systems, with and not structured matrices, and with well and ill conditioned matrices.

Country of Origin
🇨🇦 Canada

Page Count
25 pages

Category
Mathematics:
Numerical Analysis (Math)