Compute distance induced by riemannian metric

Not that you are doing anything wrong, it is just that this is way too complicated and it is unclear if you will ever arrive to an answer within a reasonable timeframe if you follow that direction. Here are the steps to solve this efficiently:

1. Suppose that $M=M_1\times M_2$ is the direct product of two Riemannian manifolds with the product metric, where $M_1$ is uniquely geodesic: Any two points are connected by a unique (up to an affine reparameterization) geodesic. Show that the distance function on $M$ is given by the Pythagorean formula (maybe you already proved this in your class/textbook you are using): $$ d^2((x_1,x_2), (y_1, y_2))= d^2(x_1, y_1) + d^2(x_2, y_2). $$ (Use the fact that the LC connection on $M$ is the sum of two connections on the factors to show that geodesics in $M$ are given by $c(t)= (c_1(t), c_2(t))$, where $c_1, c_2$ are geodesics in $M_1, M_2$.)

2. Compute (or copy from somewhere, say, wikipedia) a formula for the distance function in the upper half space model of the hyperbolic 3-space; the best one I know uses cross-ratios.

3. Show that your space is isometric to the Riemannian direct product $$ {\mathbb H}^3 \times S^1. $$ Hint: Consider first the subspace of your space given by $x_2=0, x_1>0$; then rotate this half-space around the 2-plane $x_1=x_2=0$. The circle $S^1$ will be parameterized by the angle of rotation. Equivalently, write your metric in cylindrical coordinates on $X$.