Floating-Point Computation of Functions with Maximum Accuracy
Bohlender · IEEE Transactions on Computers · 1977
Algorithms are given that compute multiple sums and products and arbitrary roots of floating-point numbers with maximum accuracy, The summation algorithm can be applied to compute scalar products, matrix products, etc. For all these functions, simple error formulas and the smallest floating-point intervals containing the exact result can be obtained.