I started playing with this. We can implement `Base.sum_seq{T<:BlasFloat}(::Base.Abs2Fun, a::Array{T}, ifirst::Int, ilast::Int)` to call BLAS.dot and then we get pairwise sum with the same code.