I think the easiest way to follow through with the calculations in a case like this is to use the well-known vector calculus identities,
$\nabla \cdot (fX) = \nabla f \cdot X + f \nabla \cdot X, \tag{1}$
$\nabla \times (fX) = \nabla f \times X + f \nabla \times X, \tag{2}$
where $f$ is a scalar function and $X$ is a vector field, both assumed to be sufficiently differentiable.
We also have, for $\mathbf r = (x, y, z)^T$,
$\nabla \cdot \mathbf r = 3, \tag{3}$
and
$\nabla \times \mathbf r = 0; \tag{4}$
the first of these, (3), is very east to calculate; the second, (4), follows from the fact that the $x$ component of $\mathbf r$, $x$, exhibits no dependence on $y$ or $z$ and so forth; see my answer to this question. Now with $f(x, y, z)$ as in the text of the present problem,
$f(x, y, z) = (x^2+y^2) \log(1-z), \tag{5}$
we have
$\nabla f = \begin{pmatrix} 2x \log(1-z) \\ 2y \log(1-z) \\ -\dfrac{x^2 + y^2}{1 - z} \end{pmatrix}; \tag{6}$
it follows from (1) and (3) that
$\nabla \cdot (f \mathbf r) = \nabla f \cdot \mathbf r + 3 f$
$= 2x^2\log(1-z) + 2y^2 \log(1-z) - \dfrac{(x^2 + y^2)z}{1 - z} + 3(x^2 + y^2) \log(1 - z)$
$= (5\log(1 - z) - \dfrac{z}{1 - z})(x^2 + y^2); \tag{7}$
as for $\nabla \times f \mathbf r$, from (2) and (4),
$\nabla \times f \mathbf r = \nabla f \times \mathbf r; \tag{9}$
in (9), we have reduced the problem to an ordinary vector cross product; I trust my readers can finish this calculation if they need to know the result.