The question starts with the well known facts that: if $f$ is a smooth function on $S^1$, then its Fourier series converges to it in smooth topology.
This must be true in more general setting. I have the following generalizations in mind.
let $(M,g)$ be a closed Riemannian manifold, let $\phi_n$ be the eigenfunctions of the Laplace operator, counting multiplicity and normalized. If $c_n=\int_M f \phi_n dV_g$, does the series $\sum_n c_n \phi_n$ converges to $f$ in smooth topology?
Let $E$ be some vector bundle over $(M,g)$ with a reasonable metric and connection. There are some elliptic operators whose eigensections form an orthonormal basis of $L^2$. Say, $\phi_n$ is a basis, normalized. Repeat the above question for a smooth section $f$ of $E$. In particular, what about differential forms, i.e. $E=\Lambda^k(M)$, and the Hodge Lapalce operator?
My guess is that these are true, but can not find any reference. A proof would require good estimates on variaous derivatives of the eigenfunctions (and sections in the second case), which I do not know how to prove.