Thanks for the interesting writeup. I wonder if it is because glibc has a less-optimized atan2f (for floats). The double version seems quite involved in glibc

Thanks for the interesting writeup. I wonder if it is because glibc has a less-optimized atan2f (for floats). The double version seems quite involved in glibc anyway: https://sourceware.org/git/?p=glibc.git;a=blob;f=sysdeps/ieee754/dbl-64/e_atan2.c;h=b77ad2dc6b9561e6396a76db11f6db686f170bb1;hb=ae37d06c7d127817ba43850f0f898b793d42aea7 https://sourceware.org/git/?p=glibc.git;a=blob;f=sysdeps/iee... ...where the atan2f version ends up calling atan, which doesn't seem as sophisticated: https://sourceware.org/git/?p=glibc.git;a=blob;f=sysdeps/ieee754/flt-32/s_atanf.c;h=f61b0f5ef64866fb053a87f9c1a968165ba4de2a;hb=ae37d06c7d127817ba43850f0f898b793d42aea7 https://sourceware.org/git/?p=glibc.git;a=blob;f=sysdeps/iee...