On Tue, 5 Dec 2017, Paul Clarke wrote: > + unsigned int n = ((unsigned int) (abstheta * inv_PI_4)) + 1.0; > + double x = n / 2.0; I'd expect adding integer 1 and dividing by integer 2, as in the current sinf version, unless there's some reason that's incorrect in this case. -- Joseph S. Myers joseph@codesourcery.com