Fourier–Galerkin method for 2D solitons of Boussinesq equation

We develop a Fourier-Galerkin spectral technique for computing the stationary solutions of 2D generalized wave equations. To this end a special complete orthonormal system of functions in L^2(-~,~) is used for which product formula is available. The exponential rate of convergence is shown. As a featuring example we consider the Proper Boussinesq Equation (PBE) in 2D and obtain the shapes of the stationary propagating localized waves. The technique is thoroughly validated and compared to other numerical results when possible.