Scale up the summing of local models

[Rolf Jagerman]

Currently models are summed by returning and adding the β matrices element-wise. This only really works for small β matrices. To scale this up we would want to emit all the elements (j,k) of β, and use a map-reduce style summation.