Hacker News Reader

Comment by xscott | Hacker News Reader

The hypot() function avoids many overflow and underflow cases:

    return hypot(hypot(q.w, q.x), hypot(q.y, q.z));

Interestingly, with clang and -O3 I get identical assembly for hypot() and the "naive" implementation.