pub(crate) fn rem_pio2f(x: f32) -> (i32, f64)
Return the remainder of x rem pi/2 in *y
use double precision for everything except passing x use __rem_pio2_large() for large x