Let $g$ be a smooth Riemannian metric on the closed $n$-dimensional unit disk $\mathbb{D}^n$. Let $f$ be a harmonic function w.r.t $g$.
Is it true that $f$ must be real-analytic?
I think that this is true if we assume that $g$ is real-analytic, but I am not sure. Is it true in that case? I would like to find a reference.
This should be related to whether or not the Riemannian laplacian $\Delta_g$ is "analytically hypoelliptic".
Not if $g$ is just smooth: take any nonconstant harmonic function $f_0$ on $\mathbb{D}^n$ (in the usual sense, i.e. with respect to the Euclidean metric) and let $T:\mathbb{D}^n\to\mathbb{D}^n$ be a diffeomorphism such that $f=f_0\circ T$ is not real-analytic. Transporting the Euclidean metric along $T$, we get a smooth metric $g$ such that $f$ is harmonic with respect to $g$.
I don't know anything about what you can say if $g$ must be real-analytic.