Computes square root of a*a + b*b without destructive overflow or underflow.