Abstract
In this paper, a new iterative refinement for ill-conditioned linear systems is derived based on discrete gradient methods for gradient systems. It is proved that the new method is convergent for any initial values irrespective of the choice of the stepsize h. Some properties of the new iterative refinement are presented. It is shown that the condition number of the coefficient matrix in the linear system to be solved in every step can be improved significantly compared with Wilkinson’s iterative refinement. The numerical experiments illustrate that the new method is more effective and efficient than Wilkinson’s iterative refinement when dealing with ill-conditioned linear systems.
Similar content being viewed by others
References
Golub, G.H., Van Loan, C.F.: Matrix Computations. The John Hapkins University Press, Baltimore (1996)
Volokh, K.Y., Vilnay, O.: Pin-pointing solution of ill-conditioned square systems of linear equations. Appl. Math. Lett. 13, 119–124 (2000)
Varah, J.M.: On the numerical solution of ill-conditioned linear system with applications to ill-posed problems. SIAM J. Numer. Anal. 10(2), 257–267 (1973)
Hansen, P.: Regularization tools: a MATLAB package for analysis and solution of discrete ill-posed problems. Numer. Algorithms 6, 1–35 (1994)
Tikhonov, A.N., Arsenin, V.Y.: Solutions of Ill-posed Problems. Winston and Sons, Washignton D.C (1977)
Liu, C.-S., Atluri, S.N.: An iterative method using an optimal desent vector, for solving an ill-conditioned system \(Bx=b\), better and faster than the conjugate gradient method. Comput. Model. Eng. Sci. 80(4), 275–298 (2011)
Buttari, A., Dongarra, J., Langou, J., Langou, J.L., Luszczek, P.: Mixed precision iterative refinement techniques for the solution of dense linear systems. Int. J. High Perform. C 21(4), 457–466 (2007)
Carson, E., Higham, N.J.: Accelerating the solution of linear systems by iterative refinement in three precisions. SIAM J. Sci. Comput. 40(2), A817–A847 (2018)
Carson, E., Higham, N.J.: A new analysis of iterative refinement and its application to accurate solution of ill-conditioned sparse linear systems. SIAM J. Sci. Comput. 39(6), A2834–A2856 (2017)
Chu, M.T.: On the continuous realization of iterative processes. SIAM Rev. 30, 375–387 (1988)
Chu, M.T.: Linear algebra algorithms as dynamical systems. Acta Numer. 19, 1–86 (2008)
Ramm, A.G.: Dynamical systems method for solving operator equations. Commun. Nonlinear Sci. Numer. Simul. 9, 383–402 (2004)
Wu, X., Shao, R., Zhu, Y.: New iterative improvement of a solution for an ill-conditioned system of linear equations based on a linear dynamic system. Comput. Math. Appl. 44, 1109–1116 (2002)
Wu, X., Wang, Z.: A new iterative refinement with roundoff error analysis. Numer. Linear Algebr. 18, 275–282 (2011)
Wu, X., Fang, Y.: Wilkinson’s iterative refinement of solution with automatic step-size control for linear system of equations. Appl. Math. Comput. 193, 506–513 (2007)
Hairer, E., Lubich, C., Wanner, G.: Geometric numerical integration: structure-preserving algorithms for ordinary differential equations. Springer Science and Business Media 31(2) (2006)
McLachlan, R.I., Quispel, G.R.W.: Splitting methods. Acta Numer. 11, 341–434 (2002)
Miyatake, Y., Sogabe, T., Zhang, S.: On the equivalence between SOR-type methods for linear systems and the discrete gradient methods for gradient systems. J. Comput. Appl. Math. 342, 58–69 (2018)
McLachlan, R.I., Quispel, G.R.W., Robidoux, N.: Geometric integration using discrete gradient. R. Soc. Lond. Philos. Trans. Ser. A Math. Phys. Eng. Sci. 357, 1021–1045 (1999)
Quispel, G.R.W., McLaren, D.I.: A new class of energy-preserving numerical integration methods. J. Phys. A 41, 045206 (2008)
Quispel, G.R.W., Turner, G.S.: Discrete gradient methods for solving ODEs numerically while preserving a first integral. J. Phys. A 29, L341–L349 (1996)
Gonzalez, O.: Time integration and discrete Hamilton systems. J. Nonlinear Sci. 6, 449–467 (1996)
Itoh, T., Abe, K.: Hamiltonian-conserving discrete canonical equations based on variational difference quotients. J. Comput. Phys. 76, 85–102 (1988)
Harten, A., Lax, P.D., van Leer, B.: On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev. 25, 35–61 (1983)
Acknowledgements
The authors sincerely thank the editors and referees for their kind and valuable comments of revision which improved the presentation of the paper.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The research was supported in part by the Natural Science Foundation of China under Grant 11701271 and by the Natural Science Foundation of the Jiangsu Higher Education Institutions under Grant 16KJB110010.
About this article
Cite this article
Liu, K., Yang, J. & Liu, C. A new iterative refinement for ill-conditioned linear systems based on discrete gradient. Japan J. Indust. Appl. Math. 37, 803–818 (2020). https://doi.org/10.1007/s13160-020-00417-z
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13160-020-00417-z
Keywords
- Wilkinson’s iterative refinement
- Ill-conditioned system of linear equations
- Discrete gradient method
- Gradient system