inline fn hypotUnfused(comptime F: type, x: F, y: F) F F: typex: Fy: F71
72
73
74
75
76
77
78
79
80
81
82inline fn hypotUnfused(comptime F: type, x: F, y: F) F {
const r = @sqrt(x * x + y * y);
if (r <= 2 * y) { // 30deg or steeper
const dx = r - y;
const z = x * (2 * dx - x) + (dx - 2 * (x - y)) * dx;
return r - z / (2 * r);
} else { // shallower than 30 deg
const dy = r - x;
const z = 2 * dy * (x - 2 * y) + (4 * dy - y) * y + dy * dy;
return r - z / (2 * r);
}
}