Error estimation of floating-point summation and dot product
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.
floating-point summation dot product rounding unit in the last place (ulp) unit in the first place (ufp) error analysis error bounds