code: def _log_multivariate_normal_density_tied(X, means, covars):
      
      cv = np.tile(covars, (means.shape[0], 1, 1))
      return _log_multivariate_normal_density_full(X, means, cv)
docs: Compute Gaussian log-density at X for a tied model.