Tweaks to sum across dimensions

- Implement `sum(fn, A, region)`
- Implement `Base.sumabs` across dimensions
- Use pairwise summation/BLAS for sums across first dimension