added the option to compute the "exact optimal decomposition" including the full dependence on the number of frequency bins and the number of polarizations used in pareto and pareto_amplitudes. This is not used anywhere else, although the natural place would be in posteriors.log_posterior via a delegation to hPrior. Possibly extending the hPrior object to a specifice pareto_hPrior object?