Green's function of a differential operator contains a lot of information of that operator. In particular, if we have a differential operator on a compact manifold with discrete spectrum, then Green's function satisfies
$$G(x,x',\lambda)=\sum_n\frac{\psi_n^*(x)\psi_n(x)}{\lambda-\lambda_n}$$
where $\lambda_n$ are the eigenvalues of the differential operator, $\psi_n$ the corresponding eigenfunctions, and $\psi_n^*$ the complex conjugate.
In theory, we can look where the poles are in a Green's function and find the spectrum, and then calculate the residue of the poles to get the eigenfunctions. However, when an eigenvalue has multiplicity, then the residue might not give us the eigenfunctions in a straight-forward way. In this situation, is there any way to work around this difficulty and calculate the eigenfucntions?
More specifically, how can we write down a basis of the eigenspace associated to an eigenvalue, say 6, of the Laplacian on $S^3$ in terms of the Green's function of this operator?