plot all 91 models
added MPI capability to measure_features
rename compute_constraints --> likelihood
test with power spectrum only