Merge branch 'reduction_kernels' into pure_scalar