function c = hard(X, compParams)
  scaledpix = mixtureScaledPiX(X,compParams);
  [dummy,c] = max(scaledpix,[],2);
  