See comment above.  For multi-modal distributions this very sensible scheme may miss out large chunks of probability.  I don't think this is a shortcoming that is easily avoided, so maybe we should recommend something about this in the documentation?