Getting MPI output/combined output

[cmbant]

We don't have a very clear way to get combined output, and currently some samplers produce the same result on all processes (minimize), and some produce per-process result on each process (post/mcmc). Either add combined_products() to sampler base, or "combined" argument to products() or similar?