Consider the following integral: $$ I(x) \ = \ \int_{0}^{\infty} dp \int_{-\infty}^{\infty} dr \int_{-\infty}^{\infty} ds\ \sinh(\pi p) \left[ K_{ip} \left( \sqrt{ r^{2} + s^{2} } \ x \right) \right]^{2} $$
Where $K_{ip}$ is the modified Bessel function of the second kind, or imaginary order $\nu = i p$.
I have a strong feeling that this integral should diverge, but I am unsure how to show this explicitly. Is there a way to see this?
I feel like maybe I could use the identity $$ K_{\nu}(x) K_{\nu}(y) \ = \ \frac{1}{2} \int_{0}^{\infty} \frac{1}{2a} \exp\left( - \frac{a}{2} - \frac{x^{2} + y^{2}}{2a}\right)\ K_{\nu}\left( \frac{xy}{a} \right) \ dp $$
Should I make a variable change $(r,s) \mapsto$ to some kind of polar coordinates? in this case, the "angle" wouldn't appear in the integrand. I'm unsure what is the best way to approach such a problem.