def dist_mahalanobis(self, X):
    Q = dotABT(self.cholL(), X-self.m)
    Q *= Q
    return np.sum(Q, axis=0)
    
  ######################################################### Global updates
  #########################################################   M-step
Given NxD matrix X, compute  Nx1 vector Dist
dataset/ETHPy150Open daeilkim/refinery/refinery/bnpy/bnpy-dev/bnpy/distr/GaussDistr.py/GaussDistr.dist_mahalanobis