Usually in textbooks, $$\int_{\mathbb{R}^d} \delta(\mathbf{x}-\mathbf{y})f(\mathbf{x}) = f(\mathbf{y})$$ holds given $f$ is continuous. On the other hand, the definition of Lebesugue point $\mathbf{y}$ is that the following limit exists $$ \lim_{\epsilon\rightarrow0}\frac{1}{|B_{\epsilon}|}\int_{B_{\epsilon}}f(\mathbf{x}), $$ where $B_{\epsilon}$ is the ball of radius $\epsilon$ centered around $\mathbf{y}.$ Lebesgure differentiation theorem says for $L^1(\mathbb{R}^d)$ function almost every point is a Lebesgue point.
All these facts let me think the delta function is also well-defined almost everywhere for $f\in L^1(\mathbb{R}^d)$ because the definition of Lebesgure point is exactly an example of approximation of the delta function. More precisely, can the definition of the delta functional be extended to the domain $L^1$ in this way?
For delta function, we take this definition: The delta function is a generalized function that can be defined as the limit of a class of delta sequences. My question is to extend this functional to $L^1$ such that the extended delta function is a linear bounded functional on $L^1.$
An element of $L^{1}$ is an equivalence class of functions which are equal a.e.. Elements of $L^{1}$ don't have pointwise values. Point values can make sense for specific cases, such as when there is a function in the equivalence class which happens to be continuous--that's because if there is such a function in the equivalence class, then there cannot be a different continuous function in the same equivalence class. However, point values are not determined in general.