I think @dengemann just merged a decimation function. Might be worth seeing if his works better. If it doesn't, we should consider putting this one and his in the same place (`decimation.py`?).