'Log of the multivariate normal probability density function.
Parameters
x : array_like
Quantiles, with the last axis of `x` denoting the components.
%(_doc_default_callparams)s
Notes
%(_doc_callparams_note)s
Returns
pdf : ndarray
Log of the probability density function evaluated at `x`'

def logpdf(self, x, mean, cov):

    (dim, mean, cov) = _process_parameters(None, mean, cov)
    x = _process_quantiles(x, dim)
    (prec_U, log_det_cov) = _psd_pinv_decomposed_log_pdet(cov)
    out = self._logpdf(x, mean, prec_U, log_det_cov)
    return _squeeze_output(out)
