Let $f \in C^k_b([0,\infty); \mathbb{R})$. We want to extend $f$ to the space $\mathbb{R}$ without loosing regularity and boundedness (thanks zhw. for this observation), i.e. we want some $F \in C^k_b(\mathbb{R}; \mathbb{R})$ with $F|_{[0,\infty)} = f$. First idea to come into ones mind is to simply reflect, i.e. to set $$ F(x) = \begin{cases} f(x) &\quad x \geq 0 \\ f(-x) &\quad x < 0\end{cases} $$
This gives us $F(x) \in C(\mathbb{R}; \mathbb{R})$ but higher derivatives needn't be continuous. The idea of higher order reflections is to write $$ F(x) = \sum_{j = 1}^{k + 1} \lambda_j \; f(-jx), \quad (x < 0) $$ and determine the linear factors $\lambda_j$ such that $D^mF$ is continuous at $x= 0$ for all $m \in \{0,\dots,k\}$. This gives the following set of equations \begin{align} f(0) &= F(0) = \sum_{j = 1}^{k + 1} \lambda_j \; f(0) \\ f'(0) &= F'(0) = \sum_{j = 1}^{k + 1} (-j) \lambda_j \; f'(0) \\ &\dots \\ f^k(0) &= F^k(0) = \sum_{j = 1}^{k + 1} (-j)^k \lambda_j \; f^k(0) \end{align}
How would you motivate this procedure? The only references I found so far were this blog post and this other post which are a good start but not completely satisfying. In both posts it is said that to preserve the $C^k$ smoothness, the extension process must preserve polynomials of x of degrees up to $k$. Could someone elaborate on this fact?
Bottom line so far: An extension operator which is a convex combination of "scaling" operators", i.e. some operator $$ E(u) := \sum_{j = 1}^{k + 1} \lambda_j \; \delta^{-j}(u), $$ where $\delta^{-j}(u)(x) := u(-jx)$ seems to be an apt candidate for an extension operator that maintains differentiability. Is there some intuition to this? Why would one think this is a good idea, other than making a direct calculation.