Let $U\subseteq \Bbb R^n$ open and $f:U \to \Bbb R$ be a continuously differentiable function, such that additionally the second partial derivatives exist and are continuous in the variables $(x_r , \ldots , x_n)$.
Is there a form of mixed Taylor formula such that
$$f(x+h) = f(x) + \sum_{k=1}^{n} \frac{\partial f}{\partial x_k} (x) h_k + \frac 1 2\sum_{k,\ell = r}^n \frac{\partial ^2 f}{\partial x_k \partial x_\ell} (x) h_k h_\ell + r(h;x)$$
such that $$\frac{r(h;x)}{\sum_{k=1}^{r-1} \vert h_k\vert +\sum_{k=r}^n h_k^2} \to 0$$ for $h\to0$?
Especially I want to have some more explicit information on $r(h;x)$ such I can deduce something about the behaviour on compacta, namely that the convergence to $0$ above holds uniform in $x$ on compacta.
I think denominator can be in this form because of the intuition that the error in the direction of the first $r-1$ variables behaves as $o(\vert h\vert $) while the error in the last variables behaves as $o(\vert h\vert ^2)$.
Matching to this I tried to apply the usual Taylor approximation for the two cases separately, but for example in the case $(x_r , \ldots , x_n)$ one can either approximate around $x + (h_1 , \ldots , h_{r-1}, 0, \ldots , 0)$ or $x$ with adding $(0, \ldots , 0 , h_r , \ldots , h_n)$ but the mixed term does not lead to the desired.
Ok here is half of an answer: if $f \in C^2(U)$ your formula is correct. From the usual Taylor formula you get $r(h, x)=p(x,h) + o (\|h\|_2^2)$ where $$p(h, x) = \sum_{\substack{1 \le i < r\\ r \le j<n}}\frac{\partial^2 f}{\partial x_i \partial x_j} (x) h_ih_j + \frac{1}{2} \sum_{\substack{1 \le i < r\\ 1 \le j < r}}\frac{\partial^2 f}{\partial x_i \partial x_j} (x) h_ih_j$$ Doing some standard calculation you get
$$ |p(x,h)| \le \left(\sum_{i=1}^{r-1} |h_i| \right) \left(\sum_{j=r}^n |h_j| \max_{1 \le i <r} \left|\frac{\partial^2 f}{\partial x_i \partial x_j} (x)\right| + \frac{1}{2} \sum_{j=1}^{r-1} |h_j| \max_{1 \le i \le r} \left|\frac{\partial^2 f}{\partial x_i \partial x_j} (x)\right| \right)$$
from which you get $$\frac{r(x,h)}{\sum_{i=1}^{r-1} |h_i| + \sum_{i=r}^n |h_i|^2} \rightarrow 0 \quad \text{for} \quad h \to 0$$ as desired.