fn asinh32(x: f32) f32 x: f3229
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54fn asinh32(x: f32) f32 {
const u = @as(u32, @bitCast(x));
const i = u & 0x7FFFFFFF;
const s = u >> 31;
var rx = @as(f32, @bitCast(i)); // |x|
// |x| >= 0x1p12 or inf or nan
if (i >= 0x3F800000 + (12 << 23)) {
rx = @log(rx) + 0.69314718055994530941723212145817656;
}
// |x| >= 2
else if (i >= 0x3F800000 + (1 << 23)) {
rx = @log(2 * rx + 1 / (@sqrt(rx * rx + 1) + rx));
}
// |x| >= 0x1p-12, up to 1.6ulp error
else if (i >= 0x3F800000 - (12 << 23)) {
rx = math.log1p(rx + rx * rx / (@sqrt(rx * rx + 1) + 1));
}
// |x| < 0x1p-12, inexact if x != 0
else {
mem.doNotOptimizeAway(rx + 0x1.0p120);
}
return if (s != 0) -rx else rx;
}
asinh(x) = sign(x) * log(|x| + sqrt(x * x + 1)) ~= x - x^3/6 + o(x^5)