Efficient Calculations of Faithfully Rounded l 2 -Norms of n -Vectors
Stef Graillat, Christoph Lauter, Ping Tak Peter Tang, Naoya Yamanaka, Shin’ichi Oishi · ACM Transactions on Mathematical Software · 2015
In this article, we present an efficient algorithm to compute the faithful rounding of the l 2 -norm of a floating-point vector. This means that the result is accurate to within 1 bit of the underlying floating-point type. This algorithm does not generate overflows or underflows spuriously, but does so when the final result calls for such a numerical exception to be raised. Moreover, the algorithm is well suited for parallel implementation and vectorization. The implementation runs up to 3 times faster than the netlib version on current processors.