I have to show the identity I wrote in the title: it should be $\partial_z\partial_{\bar z}=\frac14\left(\partial_r^2+\frac1r\partial_r+\frac1{r^2}\partial_{\theta}^2\right)$ but some computation doesn't match with my book with these operators.
We know that $\partial_z\partial_{\bar z}=\frac14(\partial_x^2+\partial_y^2)$. Then we pass to polar coordinates writing $x=r\cos\theta$ and $y=r\sin\theta$, from which, obviously we have $r=\sqrt{x^2+y^2}$ and $\theta=\arctan\frac yx$. We use than the chain rule to write $$ \partial_x=\frac{\partial r}{\partial x}\partial_r+\frac{\partial\theta}{\partial x}\partial_{\theta}=\cos\theta\partial_r-\frac{\sin\theta}{r}\partial_{\theta} $$ and $$ \partial_y=\frac{\partial r}{\partial y}\partial_r+\frac{\partial\theta}{\partial y}\partial_{\theta}=\sin\theta\partial_r+\frac{\cos\theta}{r}\partial_{\theta} $$
and till here all matches with my book. My problem comes now: I used the above cited $\partial_z\partial_{\bar z}=\frac14(\partial_x^2+\partial_y^2)$ but this gives me $$ \partial_z\partial_{\bar z}=\frac14\left(\partial_r^2+\frac1{r^2}\partial_{\theta}^2\right) $$ because the "crossed" terms are equal and opposite: squaring $\partial_x$ the "crossed" term is $-\frac2r\cos\theta\sin\theta\partial_r\partial_{\theta}$, squaring $\partial_y$ the "crossed" term is $\frac2r\cos\theta\sin\theta\partial_r\partial_{\theta}$ hence in the sum they should vanish.
Where is the error?
EDIT: following Andrea's suggest, I wrote \begin{align*} \partial_x^2+\partial_y^2 &=\left(\cos\theta\partial_r-\frac{\sin\theta}{r}\partial_{\theta}\right)\left(\cos\theta\partial_r-\frac{\sin\theta}{r}\partial_{\theta}\right)+\\ &\;\;\;\;\;\,\left(\sin\theta\partial_r+\frac{\cos\theta}{r}\partial_{\theta}\right) \left(\sin\theta\partial_r+\frac{\cos\theta}{r}\partial_{\theta}\right)\\ &=\cos^2\theta\partial_r^2+\frac{\sin\theta}{r}\cos\theta\partial_{\theta} +\frac{\sin^2}{r}\partial_r+\cos\theta\frac{\sin\theta}{r^2}\partial_\theta+\\ &\;\;\;\;\;\;\sin^2\theta\partial_r^2-\frac{\cos\theta}{r}\sin\theta\partial_{\theta} +\frac{\cos^2}{r}\partial_r-\sin\theta\frac{\cos\theta}{r^2}\partial_\theta\\ =&\partial_r^2+\frac1r\partial_r \end{align*}
which is false again! I checked several times but I cannot find where my mistake is!
Many thanks!