It would be interesting to take that C version and hammer on it a bit for speed.

...and then add SIMD.