cost_likelihood_pmodel function

Cost function computing a log-likelihood for calibration of P-model parameters