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.


[2] <http://math.stanford.edu/~akshay/research/bsv.pdf>



Hejhal's algorithm and computational methods for non-classical Maass wave forms
I don't think it is practical to directly find higher rank Maass forms along the lines of what Hejhal did, because they are functions of several variables and their Fourier expansions involve multiple sums. Not to mention the need to implement the appropriate special functions that appear. Even if you have the Fourier expansion on the computer (which has been done for SL(3,Z) (by Boris Mezhericher) but not Sp(4,Z)), using the relations to make equations in the coefficients and/or the analog of integrating along horocycles is probably computationally prohibitive. 


