Options
Linear systems in matlab with zero residual
Publikationstyp
Journal Article
Date Issued
2018
Sprache
English
Author(s)
Institut
TORE-URI
Journal
Volume
26
Start Page
56
End Page
66
Citation
Reliable Computing (26): 56-66 (2018)
Scopus ID
Several examples of linear systems are given for which the residual of the computed approximation (computed in working precision) is entirely zero. Nevertheless, the computed approximation by Matlab’s backslash operator is highly inaccurate, and often no warning is issued. Examples for matrix inversion in that spirit are given as well. Moreover, it is demonstrated that rounding the exact real result to the nearest floating-point numbers may be worse than using an approximation computed by some floating-point algorithm. Concluding, some advice is given to users and developers.