Let $G(t, x) := \frac{1}{\sqrt{4 \pi t}} \exp\left( -\frac{x^2}{4 t }\right)$ for all $(t, x) \in (0, T) \times \mathbb{R}$ be the fundamental solution to the heat equation $\partial_tu = \partial_{xx}u$.
Evans PDE book (not verbatim but its clear from the context) states that $$ \lim_{\delta \downarrow 0} \int_{\mathbb{R}} G(\delta,y) f(t-\delta, x-y)~\mathrm{d}y = f(x, y), $$ if $f \in C_c((0, T) \times \mathbb{R})$. I am wondering if this holds true a.e. if $f \in L^2((0, T); L^2(\mathbb{R})) \cap L^\infty((0, T), L^\infty(\mathbb{R}))$. My first impulse was to approximate $f$ by some $g_\varepsilon \in C_c^\infty((0, T) \times \mathbb{R})$ and to use triangle inequality a couple of times. Then, the only term that I am not sure about yet is $$ \int_{\mathbb{R}} G(\delta, y) \big(f(t-\delta, x-y) - g_\varepsilon(t-\delta, x-y) \big)~\mathrm{d}y. $$ So of course, for fixed $\delta$, this gets small as $\varepsilon \downarrow 0$. But does it get small uniformly in $\delta$? Remember, that afterwards I want to let $\delta \downarrow 0$?
I would be grateful for any input.