Suppose I'm in flat spacetime, and I have some $C^\infty$ $SU(N)$ gauge field coefficient function $A_\mu^a(x)$, and I'd like to know what the value of this function is at a spacetime point $y$ not very far from $x$. Does $$A^a_\mu(y)=e^{(y-x)\cdot\partial}A^a_\mu(x)$$ give the correct result? I.e. is the connection for the gauge index trivial like it is for the spacetime index?
More out of general interest, if the above is correct, how does it generalize? I'd imagine for a curved spacetime one has $$A^a_\mu(y)=e^{(y-x)\cdot\nabla}A^a_\mu(x),$$ where $\nabla_\mu$ is the covariant derivative with metric compatible Christoffel symbol; is that correct?
Is there a physical setup in which the gauge fibre bundle(?) is curved and then one has a covariant derivative involving the color index?
EDIT: Some more context. Suppose I have a quark field $\psi(x)$ that transforms under a local gauge transformation as $\psi(x)\rightarrow V(x)\psi(x)$, where $V(x) = \exp(i\alpha^a(x)t^a)$, $\alpha^a(x)\in\mathbb R$ and $t^a\in SU(N)$. The QCD fermionic action generalizes from the free fermionic action in the way one would expect: the partial derivative becomes a covariant derivative, $\partial_\mu\rightarrow D_\mu$, where $D_\mu \equiv \partial_\mu -igA_\mu$.
Now if I want $\psi(x+\epsilon)$, I'd figure that the usual Taylor expansion in terms of partial derivatives would be generalized to covariant derivatives. To wit, I'd expect that $$\psi(x+\epsilon) = \psi(x)+\epsilon^\mu \partial_\mu \psi(x) + \mathcal O(\epsilon^2) \rightarrow \psi(x+\epsilon) = \psi(x)+\epsilon^\mu D_\mu \psi(x) + \mathcal O(\epsilon^2).$$ The latter equation is consistent with Eq. 7.25 in Greiner's QCD book.
However, let's consider the transformation properties of $\psi(x+\epsilon)$. Under a gauge transformation we have that $$\psi(x+\epsilon) \rightarrow V(x+\epsilon)\psi(x+\epsilon) = \bigl( V(x) + \epsilon^\mu\partial_\mu V(x)\bigr)\bigl(\psi(x) + \epsilon^\mu\partial_\mu \psi(x) \bigr) + \mathcal O (\epsilon^2) \\ = V(x)\psi(x) + \epsilon^\mu \bigl( \partial_\mu V(x)\bigr)\psi(x) + V(x)\epsilon^\mu\partial_\mu\psi(x) + \mathcal O (\epsilon^2).$$ Now consider the transformation of $$\psi(x+\epsilon) = \psi(x)+\epsilon^\mu\partial_\mu\psi(x) + \mathcal O(\epsilon^2) \\ \rightarrow V(x)\psi(x) + \epsilon^\mu\partial_\mu\bigl(V(x)\psi(x)\bigr) + \mathcal O(\epsilon^2) \\ = V(x)\psi(x) + \epsilon^\mu \bigl( \partial_\mu V(x)\bigr)\psi(x) + V(x)\epsilon^\mu\partial_\mu\psi(x) + \mathcal O (\epsilon^2).$$
On the other hand, if we take as Greiner says $\psi(x+\epsilon) = \psi(x)+\epsilon^\mu D_\mu\psi(x) + \mathcal O(\epsilon^2)$, then we'll pick up extra $A_\mu$ terms that don't cancel.
(Parenthetically, even if we decided to change all partial derivatives to covariant derivatives in the above, we still don't get the right transformation properties. To wit suppose that $$\psi(x+\epsilon)\rightarrow V(x+\epsilon)\psi(x+\epsilon) \\ = \bigl[ V(x)+\epsilon^\mu D_\mu V(x) \bigr] \bigl[ \psi(x) + \epsilon^\mu D_\mu\psi(x) \bigr] + \mathcal O(\epsilon^2) \\ = V(x)\psi(x) + \epsilon^\mu D_\mu\bigl( V(x)\psi(x) \bigr) + \mathcal O(\epsilon^2).$$ At the same time, $$\psi(x+\epsilon) = \psi(x) + \epsilon^\mu D_\mu\psi(x) + \mathcal O(\epsilon^2) \\ \rightarrow V(x)\psi(x) + \epsilon^\mu V(x)D_\mu\psi(x)+\mathcal O(\epsilon^2),$$ and clearly the above two expressions are different.)
Perhaps, then, a better way to phrase the above question is: why doesn't the Taylor expansion for a gauge theory field not generalize from partial derivatives to covariant derivatives as one would naively expect (and as written in at least one book)?