LAPACK routine to calculate sqrt(x^2 + y^2 + z^2) in a manner which doesn't cause unnecessary overflow. - View it on GitHub