ENH: new function for decimating surfces

This allows to decimate surfaces, e.g., to reduce computation time when mapping (many) digitization points to the high-resolution scalp.
Not yet ready, sometimes crashes for some reason.