Cover image for Accuracy and stability of numerical algorithms
Title:
Accuracy and stability of numerical algorithms
Personal Author:
Edition:
2nd ed.
Publication Information:
Philadelphia, PA : Society for Industrial and Applied Mathematics, 2002
ISBN:
9780898715217

Available:*

Library
Item Barcode
Call Number
Material Type
Item Category 1
Status
Searching...
30000010134104 QA297 H53 2002 Open Access Book Book
Searching...

On Order

Summary

Summary

This book gives a thorough, up-to-date treatment of the behaviour of numerical algorithms in finite precision arithmetic. It combines algorithmic derivations, perturbation theory, and rounding error analysis, all enlivened by historical perspective and informative quotations. The coverage of the first edition has been expanded and updated, involving numerous improvements. Two new chapters treat symmetric indefinite systems and skew-symmetric systems, and nonlinear systems and Newton's method. Twelve new sections include coverage of additional error bounds for Gaussian elimination, rank revealing LU factorizations, weighted and constrained least squares problems, and the fused multiply-add operation found on some modern computer architectures. This new edition is a suitable reference for an advanced course and can also be used at all levels as a supplementary text from which to draw examples, historical perspective, statements of results, and exercises. In addition the thorough indexes and extensive, up-to-date bibliography are in a readily accessible form.


Author Notes

Nicholas J. Higham is Richardson Professor of Applied Mathematics at the University of Manchester, England. He is the author of more than 80 publications and is a member of the editorial boards of Foundations of Computational Mathematics , the IMA Journal of Numerical Analysis, Linear Algebra and Its Applications , and the SIAM Journal on Matrix Analysis and Applications . His book Handbook of Writing for the Mathematical Sciences (second edition) was published by SIAM in 1998, and his book MATLAB Guide , co-authored with Desmond J. Higham, was published by SIAM in 2000.


Reviews 1

Choice Review

Higham provides a comprehensive exposition of the behavior of numerical algorithms in finite precision arithmetic. The book is highly technical and requires adequate working knowledge of numerical linear algebra and mathematical analysis. Derivation of algorithms, perturbation theory, and rounding error analysis are presented in a unified form, and the use of software tools, especially LAPACK and MATLAB, are emphasized throughout. Comparable to James H. Wilkinson's Rounding Errors in Algebraic Processes (CH, Mar'65), it is unique in the sense that motivation to problems and algorithmic derivations are given in a succinct form: although rounding errors, for example, can be difficult to understand, Higham attempts to treat this topic clearly and concisely. Extensive citations and historical perspectives make it a suitable resource for an advanced course in numerical linear algebra. Recommended. Graduate; faculty; professional. D. E. Bentil University of Vermont


Table of Contents

List of Figuresp. xvii
List of Tablesp. xix
Preface to Second Editionp. xxi
Preface to First Editionp. xxv
About the Dedicationp. xxix
1 Principles of Finite Precision Computationp. 1
1.1 Notation and Backgroundp. 2
1.2 Relative Error and Significant Digitsp. 3
1.3 Sources of Errorsp. 5
1.4 Precision Versus Accuracyp. 6
1.5 Backward and Forward Errorsp. 6
1.6 Conditioningp. 8
1.7 Cancellationp. 9
1.8 Solving a Quadratic Equationp. 10
1.9 Computing the Sample Variancep. 11
1.10 Solving Linear Equationsp. 12
1.11 Accumulation of Rounding Errorsp. 14
1.12 Instability Without Cancellationp. 14
1.13 Increasing the Precisionp. 17
1.14 Cancellation of Rounding Errorsp. 19
1.15 Rounding Errors Can Be Beneficialp. 22
1.16 Stability of an Algorithm Depends on the Problemp. 24
1.17 Rounding Errors Are Not Randomp. 25
1.18 Designing Stable Algorithmsp. 26
1.19 Misconceptionsp. 28
1.20 Rounding Errors in Numerical Analysisp. 28
1.21 Notes and Referencesp. 28
Problemsp. 31
2 Floating Point Arithmeticp. 35
2.1 Floating Point Number Systemp. 36
2.2 Model of Arithmeticp. 40
2.3 IEEE Arithmeticp. 41
2.4 Aberrant Arithmeticsp. 43
2.5 Exact Subtractionp. 45
2.6 Fused Multiply-Add Operationp. 46
2.7 Choice of Base and Distribution of Numbersp. 47
2.8 Statistical Distribution of Rounding Errorsp. 48
2.9 Alternative Number Systemsp. 49
2.10 Elementary Functionsp. 50
2.11 Accuracy Testsp. 51
2.12 Notes and Referencesp. 52
Problemsp. 57
3 Basicsp. 61
3.1 Inner and Outer Productsp. 62
3.2 The Purpose of Rounding Error Analysisp. 65
3.3 Running Error Analysisp. 65
3.4 Notation for Error Analysisp. 67
3.5 Matrix Multiplicationp. 69
3.6 Complex Arithmeticp. 71
3.7 Miscellanyp. 73
3.8 Error Analysis Demystifiedp. 74
3.9 Other Approachesp. 76
3.10 Notes and Referencesp. 76
Problemsp. 77
4 Summationp. 79
4.1 Summation Methodsp. 80
4.2 Error Analysisp. 81
4.3 Compensated Summationp. 83
4.4 Other Summation Methodsp. 88
4.5 Statistical Estimates of Accuracyp. 88
4.6 Choice of Methodp. 89
4.7 Notes and Referencesp. 90
Problemsp. 91
5 Polynomialsp. 93
5.1 Horner's Methodp. 94
5.2 Evaluating Derivativesp. 96
5.3 The Newton Form and Polynomial Interpolationp. 99
5.4 Matrix Polynomialsp. 102
5.5 Notes and Referencesp. 102
Problemsp. 104
6 Normsp. 105
6.1 Vector Normsp. 106
6.2 Matrix Normsp. 107
6.3 The Matrix p-Normp. 112
6.4 Singular Value Decompositionp. 114
6.5 Notes and Referencesp. 114
Problemsp. 115
7 Perturbation Theory for Linear Systemsp. 119
7.1 Normwise Analysisp. 120
7.2 Componentwise Analysisp. 122
7.3 Scaling to Minimize the Condition Numberp. 125
7.4 The Matrix Inversep. 127
7.5 Extensionsp. 128
7.6 Numerical Stabilityp. 129
7.7 Practical Error Boundsp. 130
7.8 Perturbation Theory by Calculusp. 132
7.9 Notes and Referencesp. 132
Problemsp. 134
8 Triangular Systemsp. 139
8.1 Backward Error Analysisp. 140
8.2 Forward Error Analysisp. 142
8.3 Bounds for the Inversep. 147
8.4 A Parallel Fan-In Algorithmp. 149
8.5 Notes and Referencesp. 151
Problemsp. 153
9 LU Factorization and Linear Equationsp. 157
9.1 Gaussian Elimination and Pivoting Strategiesp. 158
9.2 LU Factorizationp. 160
9.3 Error Analysisp. 163
9.4 The Growth Factorp. 166
9.5 Diagonally Dominant and Banded Matricesp. 170
9.6 Tridiagonal Matricesp. 174
9.7 More Error Boundsp. 176
9.8 Scaling and Choice of Pivoting Strategyp. 177
9.9 Variants of Gaussian Eliminationp. 179
9.10 A Posteriori Stability Testsp. 180
9.11 Sensitivity of the LU Factorizationp. 181
9.12 Rank-Revealing LU Factorizationsp. 182
9.13 Historical Perspectivep. 183
9.14 Notes and Referencesp. 187
Problemsp. 192
10 Cholesky Factorizationp. 195
10.1 Symmetric Positive Definite Matricesp. 196
10.2 Sensitivity of the Cholesky Factorizationp. 201
10.3 Positive Semidefinite Matricesp. 201
10.4 Matrices with Positive Definite Symmetric Partp. 208
10.5 Notes and Referencesp. 209
Problemsp. 211
11 Symmetric Indefinite and Skew-Symmetric Systemsp. 213
11.1 Block LDL[superscript T] Factorization for Symmetric Matricesp. 214
11.2 Aasen's Methodp. 222
11.3 Block LDL[superscript T] Factorization for Skew-Symmetric Matricesp. 225
11.4 Notes and Referencesp. 226
Problemsp. 228
12 Iterative Refinementp. 231
12.1 Behaviour of the Forward Errorp. 232
12.2 Iterative Refinement Implies Stabilityp. 235
12.3 Notes and Referencesp. 240
Problemsp. 242
13 Block LU Factorizationp. 245
13.1 Block Versus Partitioned LU Factorizationp. 246
13.2 Error Analysis of Partitioned LU Factorizationp. 249
13.3 Error Analysis of Block LU Factorizationp. 250
13.4 Notes and Referencesp. 256
Problemsp. 257
14 Matrix Inversionp. 259
14.1 Use and Abuse of the Matrix Inversep. 260
14.2 Inverting a Triangular Matrixp. 262
14.3 Inverting a Full Matrix by LU Factorizationp. 267
14.4 Gauss-Jordan Eliminationp. 273
14.5 Parallel Inversion Methodsp. 278
14.6 The Determinantp. 279
14.7 Notes and Referencesp. 281
Problemsp. 283
15 Condition Number Estimationp. 287
15.1 How to Estimate Componentwise Condition Numbersp. 288
15.2 The p-Norm Power Methodp. 289
15.3 LAPACK 1-Norm Estimatorp. 292
15.4 Block 1-Norm Estimatorp. 294
15.5 Other Condition Estimatorsp. 295
15.6 Condition Numbers of Tridiagonal Matricesp. 299
15.7 Notes and Referencesp. 301
Problemsp. 303
16 The Sylvester Equationp. 305
16.1 Solving the Sylvester Equationp. 307
16.2 Backward Errorp. 308
16.3 Perturbation Resultp. 313
16.4 Practical Error Boundsp. 315
16.5 Extensionsp. 316
16.6 Notes and Referencesp. 317
Problemsp. 318
17 Stationary Iterative Methodsp. 321
17.1 Survey of Error Analysisp. 323
17.2 Forward Error Analysisp. 325
17.3 Backward Error Analysisp. 330
17.4 Singular Systemsp. 331
17.5 Stopping an Iterative Methodp. 335
17.6 Notes and Referencesp. 337
Problemsp. 337
18 Matrix Powersp. 339
18.1 Matrix Powers in Exact Arithmeticp. 340
18.2 Bounds for Finite Precision Arithmeticp. 346
18.3 Application to Stationary Iterationp. 351
18.4 Notes and Referencesp. 351
Problemsp. 352
19 QR Factorizationp. 353
19.1 Householder Transformationsp. 354
19.2 QR Factorizationp. 355
19.3 Error Analysis of Householder Computationsp. 357
19.4 Pivoting and Row-Wise Stabilityp. 362
19.5 Aggregated Householder Transformationsp. 363
19.6 Givens Rotationsp. 365
19.7 Iterative Refinementp. 368
19.8 Gram-Schmidt Orthogonalizationp. 369
19.9 Sensitivity of the QR Factorizationp. 373
19.10 Notes and Referencesp. 374
Problemsp. 378
20 The Least Squares Problemp. 381
20.1 Perturbation Theoryp. 382
20.2 Solution by QR Factorizationp. 384
20.3 Solution by the Modified Gram-Schmidt Methodp. 386
20.4 The Normal Equationsp. 386
20.5 Iterative Refinementp. 388
20.6 The Seminormal Equationsp. 391
20.7 Backward Errorp. 392
20.8 Weighted Least Squares Problemsp. 395
20.9 The Equality Constrained Least Squares Problemp. 396
20.10 Proof of Wedin's Theoremp. 400
20.11 Notes and Referencesp. 402
Problemsp. 405
21 Underdetermined Systemsp. 407
21.1 Solution Methodsp. 408
21.2 Perturbation Theory and Backward Errorp. 409
21.3 Error Analysisp. 411
21.4 Notes and Referencesp. 413
Problemsp. 414
22 Vandermonde Systemsp. 415
22.1 Matrix Inversionp. 416
22.2 Primal and Dual Systemsp. 418
22.3 Stabilityp. 423
22.4 Notes and Referencesp. 428
Problemsp. 430
23 Fast Matrix Multiplicationp. 433
23.1 Methodsp. 434
23.2 Error Analysisp. 438
23.3 Notes and Referencesp. 446
Problemsp. 448
24 The Fast Fourier Transform and Applicationsp. 451
24.1 The Fast Fourier Transformp. 452
24.2 Circulant Linear Systemsp. 454
24.3 Notes and Referencesp. 456
Problemsp. 457
25 Nonlinear Systems and Newton's Methodp. 459
25.1 Newton's Methodp. 460
25.2 Error Analysisp. 461
25.3 Special Cases and Experimentsp. 462
25.4 Conditioningp. 464
25.5 Stopping an Iterative Methodp. 467
25.6 Notes and Referencesp. 468
Problemsp. 469
26 Automatic Error Analysisp. 471
26.1 Exploiting Direct Search Optimizationp. 472
26.2 Direct Search Methodsp. 474
26.3 Examples of Direct Searchp. 477
26.4 Interval Analysisp. 481
26.5 Other Workp. 484
26.6 Notes and Referencesp. 486
Problemsp. 487
27 Software Issues in Floating Point Arithmeticp. 489
27.1 Exploiting IEEE Arithmeticp. 490
27.2 Subtleties of Floating Point Arithmeticp. 493
27.3 Cray Peculiaritiesp. 493
27.4 Compilersp. 494
27.5 Determining Properties of Floating Point Arithmeticp. 494
27.6 Testing a Floating Point Arithmeticp. 495
27.7 Portabilityp. 496
27.8 Avoiding Underflow and Overflowp. 499
27.9 Multiple Precision Arithmeticp. 501
27.10 Extended and Mixed Precision BLASp. 503
27.11 Patriot Missile Software Problemp. 503
27.12 Notes and Referencesp. 504
Problemsp. 505
28 A Gallery of Test Matricesp. 511
28.1 The Hilbert and Cauchy Matricesp. 512
28.2 Random Matricesp. 515
28.3 "Randsvd" Matricesp. 517
28.4 The Pascal Matrixp. 518
28.5 Tridiagonal Toeplitz Matricesp. 521
28.6 Companion Matricesp. 522
28.7 Notes and Referencesp. 523
Problemsp. 525
A Solutions to Problemsp. 527
B Acquiring Softwarep. 573
B.1 Internetp. 574
B.2 Netlibp. 574
B.3 MATLABp. 575
B.4 NAG Library and NAGWare F95 Compilerp. 575
C Program Librariesp. 577
C.1 Basic Linear Algebra Subprogramsp. 578
C.2 EISPACKp. 579
C.3 LINPACKp. 579
C.4 LAPACKp. 579
D The Matrix Computation Toolboxp. 583
Bibliographyp. 587