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

def logpdf(self, x, mean=None, cov=1, allow_singular=False):

    (dim, mean, cov) = self._process_parameters(None, mean, cov)
    x = self._process_quantiles(x, dim)
    psd = _PSD(cov, allow_singular=allow_singular)
    out = self._logpdf(x, mean, psd.U, psd.log_pdet, psd.rank)
    return _squeeze_output(out)
