Abstract
We improve the well-known Wilkinson-type estimates for the error of standard floating-point recursive summation and dot product by up to a factor 2. The bounds are valid when computed in rounding to nearest, no higher order terms are necessary, and they are best possible. For summation there is no restriction on the number of summands. The proofs are short by using a new tool for the estimation of errors in floating-point computations which cures drawbacks of the “unit in the last place (ulp)”. The presented estimates are nice and simple, and closer to what one may expect.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have