The method of separation of variables could be used as shown below:
$$\frac{\partial^2 g(x,\xi)}{\partial w \partial \xi}=A g(x,\xi),$$
where $w=\ln(1/x) \quad\to\quad dw=-\frac{1}{x}dx\quad\to\quad \frac{\partial g}{\partial w}=\frac{\partial g}{\partial x}\frac{dx}{dw}=-x\frac{\partial g}{\partial x}$
$$-x\frac{\partial^2 g(x,\xi)}{\partial x \partial \xi}=Ag(x,\xi),$$
$g=F(x)G(\xi) \quad\to\quad -x\frac{F'}{F}\frac{G'}{G}=A \quad\to\quad \begin{cases} -x\frac{F'}{F}=\lambda \quad\to\quad F=x^{-\lambda} \\ \frac{G'}{G}=\frac{A}{\lambda} \quad\to\quad G=e^{A\xi/\lambda} \end{cases}$
This leads to a family of particular solutions of the PDE : $\quad g_{\lambda}(x,\xi)=C_\lambda e^{A\xi/\lambda}x^{-\lambda}$
Combining all these particular solutions leads to a general form expressing the solution of the PDE:
$$g(x,\xi)=\int\Phi(\lambda) e^{A\xi/\lambda}x^{-\lambda}d\lambda =\int\Phi(\lambda) e^{A\xi/\lambda-\lambda\ln{x}} d\lambda $$
$$g(x,\xi)=\int\Phi(\lambda) e^{\frac{A}{\lambda}\xi+\lambda w} d\lambda $$
with any function $\Phi$ insofar the integral be convergent. This function $\Phi$ has to be determined according to the boundary conditions. Generally that is the most difficult part of the job, depending a lot of what are the boundary conditions.
To the question : 1.Why don't they use separation of variables for this PDE? , the answer might be : Because this form of function to express the PDE solution is complicated and not easy to approximate according to the context of the problem.
You wrote : "The book provides a solution only in the limit $w\xi≫1$ for some reason". Moreover, the boundary conditions are not explicitly given. All this makes difficult to derive an approximate from the above general formula. That is probably why a simpler approach (may be for didactic reason) was investigated and then presented in the book.