fn expo2d(x: f64) f64 x: f6428
29
30
31
32
33
34
35fn expo2d(x: f64) f64 {
const k: u32 = 2043;
const kln2 = 0x1.62066151ADD8BP+10;
const u = (0x3FF + k / 2) << 20;
const scale = @as(f64, @bitCast(@as(u64, u) << 32));
return @exp(x - kln2) * scale * scale;
}