*
* return the remainder of x rem pi/2 in *y
* use double precision for everything except passing x
- * use __rem_pio2_slow() for large x
+ * use __rem_pio2_large() for large x
*/
#include "libm.h"
e0 = (ix>>23) - 150; /* e0 = ilogb(|x|)-23; */
SET_FLOAT_WORD(z, ix - ((int32_t)(e0<<23)));
tx[0] = z;
- n = __rem_pio2_slow(tx,ty,e0,1,0);
+ n = __rem_pio2_large(tx,ty,e0,1,0);
if (hx < 0) {
*y = -ty[0];
return -n;