I have this relationship:
\begin{align} \frac{1}{|x|}=f(x)*f(x)\ , \end{align} where $*$ denotes the convolution. I want to solve for $f(x)$. My first instinct was to apply the convolution theorem:
\begin{align} \mathcal{F}\left\{\frac{1}{|x|}\right\}=\mathcal{F}\left\{f(x)\right\}^2\implies f(x)=\mathcal{F}^{-1}\left\{\sqrt{\mathcal{F}\left\{\frac{1}{|x|}\right\}}\right\}\ , \end{align}
where $\mathcal{F}\{\}$ is the Fourier transform. According to Mathematica, \begin{align} \sqrt{\mathcal{F}\left(\frac{1}{|x|}\right)}=\sqrt{\frac{-2 \log (\left| \omega\right| )-2 \gamma }{\sqrt{2 \pi }}}\ , \end{align} where $\omega$ is the Fourier conjugate variable to $x$ and $\gamma$ is Euler's constant. How should I compute the inverse Fourier transform of the above function? I was thinking of using contour integration, but I am not sure what contour I should choose.
Alternatively, is there an easier way to solve for $f(x)$, possibly using functional equation methods? I saw this question, which makes me concerned that there is not a unique $f(x)$ as solution.