Agner's CPU blog

I'm working on patching an old piece of code used to reduce the angles of a 3D float vector to the range ±Pi. The original code used loops to implement a horribly inaccurate version of IEEE remainder. Replacing the loops with FPREM1 2Pi has been working well so far, but I'd really like to use SSE in...