function [p, c] = evalmixtureX(X, compParams) ;
  [scaledpix, maxlog] = mixtureScaledPiX(X, compParams) ;
  p = sum(scaledpix,2) .* exp(maxlog) ;
  [dummy, c] = max(scaledpix,[],2);