Adds a constant to each double-precision floating-point strided array element and computes the sum using pairwise summation.