Let $\Omega = [0, 1]^d$ and suppose that $f \colon \Omega \to \mathbb{R}$ lies in order $m > d/2$ Sobolev space; i.e., 
$$
\|f\|_{H^m(\Omega)}^2 = \sum_{|\alpha| \leq m} \|D^\alpha f\|_{L^2(\Omega)}^2 < \infty.
$$
Let $x_1, \dots, x_n\in \Omega$ and define the distance 
$$
h \equiv h_\Omega(x_1, \dots, x_n) = \sup_{y \in \Omega} \min_{i = 1, \dots, n} |y - x|, 
$$
where above $|\cdot|$ denotes the usual Euclidean norm on $[0, 1]^d$. 

We can then consider the minimum Sobolev norm interpolant: 
$$
s_{m, \Omega}(x_1, \dots, x_n) = 
\mathrm{argmin}_{g \in L^2(\Omega)} \Big\{\|g\|_{H^m(\Omega)} : g(x_i) = f(x_i), ~~\mbox{for}~i = 1,\dots n\Big\}
$$

I have seen the following theorem quoted in various papers/sources (Matveev, the book by Anatoly Yu. Bezhaev, Vladimir A. Vasilenko, etc.), but none containing a full proof:
 
> **Theorem:** Suppose that $f \in H^m(\Omega)$ where $\Omega = [0, 1]^d$ and $m > d/2$. Let $x_1, \dots, x_n \in \Omega$ and let $\widehat{f}$ denote the interpolant $s_{m, \Omega}(x_1, \dots, x_n)$. Then there exists a constant $C$ depending only on $\Omega, m$ such that $$
\|f - \widehat{f}\|_{L^2(\Omega)} \leq C~ h^{m}~\|f\|_{H^m(\Omega)}.
$$

Is there any source (in English) which provides a full (ideally self-contained) proof of this result?