Issue No. 08 - August (1994 vol. 43)
DOI Bookmark: http://doi.ieeecomputersociety.org/10.1109/12.295860
<p>An attractive paradigm for building fast numerical algorithms is the following: 1) try a fast but occasionally unstable algorithm, 2) test the accuracy of the computed answer, and 3) recompute the answer slowly and accurately in the unlikely event it is necessary. This is especially attractive on parallel machines where the fastest algorithms may be less stable than the best serial algorithms. Since unstable algorithms can overflow or cause other exceptions, exception handling is needed to implement this paradigm safely. To implement it efficiently, exception handling cannot be too slow. We illustrate this paradigm with numerical linear algebra algorithms from the LAPACK library.</p>
linear algebra; parallel algorithms; digital arithmetic; exception handling; convergence of numerical methods; eigenvalues and eigenfunctions; exception handling; fast numerical algorithms; parallel machines; unstable algorithms; numerical linear algebra; LAPACK library; IEEE floating point arithmetic.
J.W. Demmel, Xiaoye Li, "Faster Numerical Algorithms Via Exception Handling", IEEE Transactions on Computers, vol. 43, no. , pp. 983-992, August 1994, doi:10.1109/12.295860