pub fn sindf(x: f64) f32
pub fn sindf(x: f64) f32 { // |sin(x)/x - s(x)| < 2**-37.5 (~[-4.89e-12, 4.824e-12]). const S1 = -0x15555554cbac77.0p-55; // -0.166666666416265235595 const S2 = 0x111110896efbb2.0p-59; // 0.0083333293858894631756 const S3 = -0x1a00f9e2cae774.0p-65; // -0.000198393348360966317347 const S4 = 0x16cd878c3b46a7.0p-71; // 0.0000027183114939898219064 // Try to optimize for parallel evaluation as in __tandf.c. const z = x * x; const w = z * z; const r = S3 + z * S4; const s = z * x; return @floatCast((x + s * (S1 + z * S2)) + s * w * r); }