Abstract
In many research areas like structural mechanics, economics, meteorology, and fluid dynamics, problems are mapped to large sparse linear systems via discretization. The resulting matrices are often ill-conditioned with condition numbers of about 1016 and higher. Usually these systems are preconditioned before they are fed to an iterative solver.
Especially for ill-conditioned systems, we show that we have to be careful with these three classical steps — discretization, preconditioning, and (iterative) solving. For Krylov subspace solvers we give some detailed analysis and show possible improvements based on a multiple precision arithmetic. This special arithmetic can be easily implemented using the exact scalar product — a technique for computing scalar products of floating point vectors exactly.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
American National Standards Institute/Institute of Electrical and Electronic Engineers, New York. A Standard for Binary Floating-Point Arithmetic,1985. ANSI/IEEE Std. 754–1985.
L. Collatz. Differentialgleichungen. 6. Aufl. Teubner, Stuttgart, 1981.
G.H. Golub and C.F. van Loan. Matrix Computations. Johns Hopkins, third edition, 1996.
A. Greenbaum. Iterative Methods for Solving Linear Systems. SIAM, Philadelphia, 1997.
R. Gregory and D. Karney. A Collection of Matrices for Testing Computational Algorithms. Wiley Interscience, 1969.
R. Hammer, M. Hocks, U. Kulisch, and D. Ratz. Toolbox for Verified Computing. Springer, Berlin Heidelberg, 1993.
R. Hammer, M. Hocks, U. Kulisch, and D. Ratz. C++ Toolbox for Verified Computing. Springer, Berlin Heidelberg, 1995.
M.R. Hestenes and E. Stiefel. Methods of conjugate gradients for solving linear systems. Journal of Research of the National Bureau of Standards, 49: 409–436, 1952.
U. Kulisch, editor. Wissenschaftliches Rechnen mit Ergebnisverifikation — Eine Einführung. Akademie Verlag, Ost-Berlin, Vieweg, Wiesbaden, 1989.
U. Kulisch and W.L. Miranker. Computer Arithmetic in Theory and Practice. Academic Press, New York, 1981.
U. Kulisch and W.L. Miranker, editors. A New Approach to Scientific Computation. Academic Press, New York, 1983.
C. Lanczos. An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. Journal of Research of the National Bureau of Standards, 45: 255–282, 1950.
R. Lohner. A verified solver for linear systems with band structure. to be included in: Toolbox for Verified Computing II.
The matrix market: a visual web database for numerical matrix data.
C.C. Paige. The Computation of Eigenvalues and Eigenvectors of Very Large Sparse Matrices. PhD thesis, University of London, 1971.
B.N. Parlett. The Symmetric Eigenvalue Problem. Prentice-Hall, Englewood Cliffs, NJ., 1980.
H. Simon. The Lanczos Algorithm for Solving Symmetric Linear Systems. PhD thesis, University of California, Berkley, 1982.
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 1999 Springer Science+Business Media Dordrecht
About this chapter
Cite this chapter
Facius, A. (1999). Influences of Rounding Errors in Solving Large Sparse Linear Systems. In: Csendes, T. (eds) Developments in Reliable Computing. Springer, Dordrecht. https://doi.org/10.1007/978-94-017-1247-7_2
Download citation
DOI: https://doi.org/10.1007/978-94-017-1247-7_2
Publisher Name: Springer, Dordrecht
Print ISBN: 978-90-481-5350-3
Online ISBN: 978-94-017-1247-7
eBook Packages: Springer Book Archive