code: def _log_multivariate_normal_density_diag(x, means, covars):
      """Compute Gaussian log-density at x for a diagonal model."""
      _, n_dim = x.shape
      lpr = -0.5 * (
          n_dim * np.log(2 * np.pi)
          + np.sum(np.log(covars), 1)
          + np.sum((means**2) / covars, 1)
          - 2 * np.dot(x, (means / covars).T)
          + np.dot(x**2, (1.0 / covars).T)
      )
      return lpr
  
file_path: madmom/madmom/ml/gmm.py
imports: from __future__ import absolute_import, division, print_function
  import numpy
  from scipy import linalg
label_desc: Use image node_id 4 for calling a global function with example usage: _log_multivariate_normal_density_diag(x, means, covars) and returns: lpr