When dealing with coupled variables, one can uncouple them using the trick of representing them as a Gaussian integral:
$$\exp \left\{\frac{b^{2}}{4a}\right\}=\int \frac{\mathrm{d} x}{\sqrt{2 \pi / a}} \exp \left\{-a x^{2} \pm b x\right\}$$
Similarly, when $b$ is complex, we can still linearise using a complex Gaussian:
$$\exp \left(\frac{|b|^{2}}{a}\right)=\int \frac{\mathrm{d} z \mathrm{d} \bar{z}}{2 \mathrm{i} \pi a} \exp \left \{-a \bar{z} z-\bar{b} z-b \bar{z}\right \}$$
However $\text{Re}[a]>0$. I have: $\exp\left(-|b|^2\right)$. Is it impossible to represent it as a gaussian integral?
(Note, I am trying to reproduce the results of Sommers et al. But it seems very unclear how they performed their integrals...cf Sommers, H. J., Crisanti, A., Sompolinsky, H., & Stein, Y. (1988). Spectrum of large random asymmetric matrices. Physical review letters, 60(19), 1895.)