pub fn sqrt_atanh2(x: i32, y: i32, z: i32) -> (i32, i32)
Returns G*sqrt(x**2 - y**2), z + atanh2(y, x)
G*sqrt(x**2 - y**2), z + atanh2(y, x)