No kidding. It occurred to me that what I should do is the same thing I did for the `_accumulate_normals`. That is, I should have an actual understandable equivalent bit of code commented to explain what's happening. Basically I started from Matti's C code, and hammered away at it until I got it to be efficient in Python, and this was the best way I could find.