vectorize operations where possible (allowing tensors, etc...)

[kratsg]

is it possible to somehow vectorize combine/reduce the distributions as part of graph optimization or is it unnecessary -- such as how pyhf has a GaussianCombined and PoissonCombined to handle all of this. Follow up from #110. See also #18.