Use the Einstein summation convention to avoid the creation of temporary arrays when computing the likelihood-weighted means and standard deviations. This speeds up the evaluation of the physical properties and ultimately pdf_analysis can run more than 25% faster when evaluating a large number of physical properties with more than a million models.