Hejhal's algorithm [1] was a little gadget invented in the 90's for calculating the Hecke eigenvalues and Fourier coefficients of Maass wave forms. Later, Booker, Strombergsson, and Venkatesh (BSV) [2] took Hejhal's paper and made it more viable to perform high-precision calculations.

The basic idea is to assume the Fourier expansion takes a certain form, and write down a finite linear system (where the finitude is possible due to assumptions on the final order of magnitude of desired precision). Now take a horocycle on the complex upper half plane, and use these points in the linear system, and you can almost solve it. In BSV's words, "to solve for r [the desired eigenvalue] the above linear system is repeatedly solved for two different Y-values [horocycles], successively adjusting r to make the two solution vectors as nearly equal as possible."

My question is whether such an algorithm or the germ of an idea for such an algorithm exists for non-classical Maass wave forms. What would take the analogue of horocycles, and would we have enough information to solve an analogous linear system? (Remember, Fourier expansions are a good deal more complicated in multiple dimensions)

[1] Dennis A. Hejhal. On eigenfunctions of the Laplacian for Hecke triangle groups. In Emerging applications of number theory (Minneapolis, MN, 1996), volume 109 of IMA Vol. Math. Appl., pages 291–315. Springer, New York, 1999.