Here is yet another solution. After the change of variable $x = -y$ with $y \to 1^-$, we have
$$ f(x) = \sum_{n=0}^{\infty} y^{4n^2}(1 - y^{4n+1}). \tag{1}$$
Now the idea is to identify $\text{(1)}$ as sort of Riemann sum for the integral $\int_{0}^{\infty} e^{-2t} \, dt = \frac{1}{2}$.
To this end, we first fix arbitrary $\epsilon \in (0, 1)$ and choose $N = N(\epsilon, y) \in \mathbb{Z}_{> 0}$ such that $y^{4N^2} \to \epsilon$ as $y \to 1^-$. For instance, choose $N$ to be any nearest integer to $(\log \epsilon / 4 \log y)^{1/2}$. Then using the alternating behavior of $\text{(1)}$ we find that
$$ \left| \sum_{n=N+1}^{\infty} y^{4n^2}(1 - y^{4n+1}) \right| \leq y^{4N^2} \xrightarrow[y \to 1^-]{} \epsilon $$
which will be shown to be negligible as $\epsilon \to 0^+$. Next, we define $t_n = t_n(y)$ by
$$t_n = \sum_{k=0}^{n} (4k+1)(1-y) = (2n^2+3n+1)(1-y) $$
and $\Delta t_n = \Delta t_n (y) := t_n - t_{n-1}$. Then we have $y = 1 - \frac{t_n}{2n^2+3n+1}$. Now from the following simple inequality
$$ (4n+1)y^{4n+1}(1-y) \leq 1 - y^{4n+1} \leq (4n+1)(1-y), $$
we find that
\begin{align*}
&\sum_{n=0}^{N} \exp\left\{ 4(n+1)^2 \log\left( 1 - \frac{t_n}{2n^2+3n+1} \right) \right\} \Delta t_n \tag{2}\\
&\qquad \leq \sum_{n=0}^{N} y^{4n^2}(1 - y^{4n+1}) \\
&\qquad \leq \sum_{n=0}^{N} \exp\left\{ 4n^2 \log\left( 1 - \frac{t_n}{2n^2+3n+1} \right) \right\} \Delta t_n \tag{3}
\end{align*}
Since $t_N \to \frac{1}{2}\log\frac{1}{\epsilon}$ and $\max_{0 \leq n \leq N} \Delta t_n = (4N+1)(1-y) \to 0$ as as $y \to 1^-$, we can think of bounds $\text{(2)}$ and $\text{(3)}$ as Riemann sums of $e^{-2t}$ over the inteval $[0, \frac{1}{2}\log\frac{1}{\epsilon}]$. Taking limit as $y \to 1^-$, it is not hard to show (but slightly annoying to give a full argument) that both bounds converge to the same value, which is
$$ \lim_{y \to 1^-} \sum_{n=0}^{N} y^{4n^2}(1 - y^{4n+1}) = \int_{0}^{\frac{1}{2}\log\frac{1}{\epsilon}} e^{-2t} \, dt = \frac{1 - \epsilon}{2}. $$
Then by the usual 3-epsilon argument, we have
$$ \limsup_{y \to 1^-} \left| f(-y) - \frac{1}{2} \right| \leq \frac{3 \epsilon}{2}
\xrightarrow[\epsilon \to 0^+]{} 0 $$
and hence the claim follows.