fn atanh_32(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
54
55
56fn atanh_32(x: f32) f32 {
const u = @as(u32, @bitCast(x));
const i = u & 0x7FFFFFFF;
const s = u >> 31;
var y = @as(f32, @bitCast(i)); // |x|
if (y == 1.0) {
return math.copysign(math.inf(f32), x);
}
if (u < 0x3F800000 - (1 << 23)) {
if (u < 0x3F800000 - (32 << 23)) {
// underflow
if (u < (1 << 23)) {
mem.doNotOptimizeAway(y * y);
}
}
// |x| < 0.5
else {
y = 0.5 * math.log1p(2 * y + 2 * y * y / (1 - y));
}
} else {
y = 0.5 * math.log1p(2 * (y / (1 - y)));
}
return if (s != 0) -y else y;
}
atanh(x) = log((1 + x) / (1 - x)) / 2 = log1p(2x / (1 - x)) / 2 ~= x + x^3 / 3 + o(x^5)