fn sinh64(x: f64) f64
fn sinh64(x: f64) f64 { const u = @as(u64, @bitCast(x)); const w = @as(u32, @intCast(u >> 32)) & (maxInt(u32) >> 1); const ax = @as(f64, @bitCast(u & (maxInt(u64) >> 1))); if (x == 0.0 or math.isNan(x)) { return x; } var h: f32 = 0.5; if (u >> 63 != 0) { h = -h; } // |x| < log(FLT_MAX) if (w < 0x40862E42) { const t = math.expm1(ax); if (w < 0x3FF00000) { if (w < 0x3FF00000 - (26 << 20)) { return x; } else { return h * (2 * t - t * t / (t + 1)); } } // NOTE: |x| > log(0x1p26) + eps could be h * exp(x) return h * (t + t / (t + 1)); } // |x| > log(DBL_MAX) or nan return 2 * h * expo2(ax); }