function [p, c] = evalmixture(X, compParams) ;
  pix = mixturePiX(X, compParams) ;
  p = sum(pix,2) ;
  [dummy, c] = max(pix,[],2);