Hyperbolic Fourier series and the Klein-Gordon equation

In an effort to extend classical Fourier theory, Hedenmalm and Montes-Rodríguez (2011) found that the function system $$e_m(x)=e^{i\pi mx},e_n^{†}(x)=e_n(-1\mathrm{/}x)=e^{-i\pi n\mathrm{/}x}$$ is weak-star complete in $L^{\mathrm{\infty }}(\mathbb{R})$ when $m,n$ range over the integers with $n\mathrm{\ne }0$. It turns out that the system can be used to provide unique representation of functions and more generally distributions on the real line $\mathbb{R}$. For instance, we may represent uniquely the unit point mass at a point $x\in \mathbb{R}$: $${\delta }_x(t)=A_0(x)+\sum _{n\mathrm{\ne }0}(A_n(x)e^{i\pi nt}+B_n(x)e^{-i\pi n\mathrm{/}t}),$$ with at most polynomial growth of the coefficients, so that the sum converges in the sense of distribution theory. In a natural sense, the system $\{A_n,B_n{\}}_n$ is biorthogonal to the initial system $\{e_n,e_n^{†}{\}}_n$ on the real line. More generally, for a distribution $f$ on the compactified real line, we may decompose it in a \emph{hyperbolic Fourier series} $$f(t)=a_0(f)+\sum _{n\mathrm{\ne }0}(a_n(f)e^{i\pi nt}+b_n(f)e^{-i\pi n\mathrm{/}t}),$$ understood to converge in the sense of distribution theory. Such hyperbolic Fourier series arise from two different considerations. One is the Fourier interpolation problem of recovering a radial function $\phi$ on ${\mathbb{R}}^d$ from partial information on $\phi$ and its Fourier transform $\widehat{\phi }$, studied by Radchenko and Viazovska (2019). Another consideration is the interpolation theory of the Klein-Gordon equation ${\mathrm{\partial }}_x{\mathrm{\partial }}_{y}u+u=0$. For instance, the biorthogonal system $\{A_n,B_n{\}}_n$ leads to a collection of solutions that vanish along the lattice-cross of points $(\pi k,0)$ and $(0,\pi l)$ save for one of these points. These interpolating solutions allow for restoration of a given solution $u$ from its values on the lattice-cross.