more sophisticated algorithm reduces the bias in the density matching estimation and ensures that all points are used, by including an extra sensitivity parameter:
# Increase each centroid's sensitivity by a small amount
# Pick a sample point at random
# For each quantization vector centroid , let denote the distance of and
# Find the centroid for which is the smallest
# Move towards by a small fraction of the distance
# Set to zero
# Repeat

It is desirable to use a cooling schedule to produce