Let $f$ be a continuous real-valued function defined on an $n$ dimesional simplex $\Sigma\subset \mathbb{R}^n $. The classical simplicial approximation scheme provides a sequence $f\_k$ of piecewise affine functions converging uniformly to $f$: namely, on each simplex $S$ of the $k$-th baricentric subdivision of $\Sigma$, one takes ${f\_k} \_ {|S}$ to be the affine interpolation of the data of $f$ on the vertices of $S$. Precisely, if $\omega$ is a modulus of continuity for $f$, it follows easily that $$ \big\| f - f\_k \big\|\_{\infty,\Sigma}\, \le\, \omega\bigg( \Big(\frac{n}{n+1}\Big )^k \operatorname{diam}\Sigma \bigg)\, .$$


Notice that $\nabla f\_k(x) $ is defined and locally constant at any point $x\in\Sigma$ not in the $(n-1)$-skeleton of the $k$-th barycentric subdivisions; in particular, a.e. in $\Sigma$. However, understanding the behaviour of this sequence is not as simple as before, even for smooth $f$. For instance, the norm of $\nabla f\_k(x) $ is in general not bounded by the uniform norm of $\nabla f $.


**Questions**. Assuming that $f$ has bounded and continuous first order derivatives, can we conclude that $\nabla f\_k $ converges to $\nabla f $ in $L^1(\Sigma)$? How to dominate pointwise $|\nabla f\_k (x) |$ (independently of $k$), and how to bound $\|\nabla f\_k -\nabla f\|\_{1,\Sigma} $ ? 



$W^{1,1}$ simplicial approximation
You can represent on every simplex $\sigma$ of your triangulation the derivative $D f\_k$ by
$$
 D f\_k = \frac{1}{\lvert \sigma \rvert} \sum\_{i = 0}^n \int\_{\sigma} \Bigl(\frac{1}{(1 - \beta\_i (x))^n} - 1\Bigr)\, Df(x)[a\_i - x]\,dx D\beta\_i,
$$
where $\beta\_i$ are the maps that give barycentric coordinates 
$$
 \beta\_0 (x) + \dots + \beta\_n (x) = 1
 \quad \text{ and } \quad 
 \beta\_0 (x)a\_0 + \dots + \beta\_n (x)a\_n = x.
$$
For a proof, see J. Van Schaftingen, *Approximation in Sobolev spaces by piecewise affine interpolation*, [[arXiv:1312.5986]](http://arxiv.org/abs/1312.5986)).


The kernel appearing in this integral representation is absolutely integrable ($a\_i - x$ controls $\beta\_i (x)$), uniformly with respect to simplices that do not degenerate. This should imply in fact an $L^\infty$ approximation property under your assumptions. (If the function $f$ is merely in a Sobolev space, one can improve the convergence by translating the triangulation points.)


