Adds a constant to each double-precision floating-point strided array element and computes the sum using an improved Kahan–Babuška algorithm.