Why aren't these two integrals equal? $$\int_0^\pi \int_{-1}^{1} e^r \,\mathrm dr\,\mathrm d\theta \ne \int_0^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta\tag1$$
The common integrand $e^r$ does not vary the same way over the two different integration domains (let's call them $S_1$ and $S_2,$ respectively), which merely have the same measure and geometric representation. Consequently, the two integrals are not guaranteed to be equal. Indeed, taking their difference: \begin{align}&\int_0^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta
-\int_0^\pi \int_{-1}^{1} e^r \,\mathrm dr\,\mathrm d\theta\\
={}& \left(\int_0^\pi \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta+\int_\pi^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta\right)
-\left(\int_0^\pi \int_{-1}^{0} e^r \,\mathrm dr\,\mathrm d\theta+\int_0^\pi \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta\right)\\
={}&\int_\pi^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta-\int_0^\pi \int_{-1}^{0} e^r \,\mathrm dr\,\mathrm d\theta\\
={}& \pi\left(\int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta
-\int_{-1}^{0} e^r \,\mathrm dr\,\mathrm d\theta\right)\\
={}&\pi\,(1.72-0.63)\\
\ne{}&0.\\\end{align}
So in theory, it seems like we should be able to integrate in polar coordinates using this region $-1 \leq r \leq 1, 0 \leq \theta < \pi$
Only $S_2,$ but not $S_1,$ is in polar coordinates.
We can consider the entire inequation $(1)$ to be residing in a “coordinate system” that is simply not isomorphic to $\mathbb R^2$ (and thus maps the given geometric region to multiple integration domains). As such, notice that to fill in this blank $$\int_0^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta = \int_0^\pi \int_{-1}^{1} \fbox{$\phantom{filler}$}\,\mathrm dr\,\mathrm d\theta,$$ just replace each instance of $r$ with $|r|.$
why does this not violate the change of variables theorem?
This theorem isn't necessary here, but can be invoked via \begin{align}x&\color{red}=r\left|\cos\theta\right|,\\y&\color{red}=r\sin\theta,\\&f\Big(g(r,\theta),h(r,\theta)\Big)\det \left| \frac{\partial(x,y)}{\partial(r,\theta)} \right|\color{red}=e^{|r|}\ne e^r;\end{align} then $$\int_0^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta
\\=\int_0^{\frac\pi2} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta
+\int_{\frac\pi2}^{\frac{3\pi}2} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta
+\int_{\frac{3\pi}2}^{2\pi} \int_{0}^{1} e^r \,\mathrm dr\,\mathrm d\theta
\\\color{red}=\int_0^{\frac\pi2} \int_{0}^{1} e^{|r|} \,\mathrm dr\,\mathrm d\theta
+\left(\int_{\frac\pi2}^\pi \int_{0}^{1} e^{|r|} \,\mathrm dr\,\mathrm d\theta
+\int_{0}^{\frac{\pi}2} \int_{-1}^{0} e^{|r|} \,\mathrm dr\,\mathrm d\theta\right)
+\int_{\frac{\pi}2}^{\pi} \int_{-1}^{0} e^{|r|} \,\mathrm dr\,\mathrm d\theta
\\=\int_0^\pi \int_{-1}^{1} e^{|r|} \,\mathrm dr\,\mathrm d\theta.$$