We need tests for the Rayleigh integral (harmonic formulation) in jwave.acoustics.time_harmonic.rayleigh_integral
It could be tested against angularSpectrum of k-Wave or using some convergence plot against a known solution (e.g. propagation of a single-frequency gaussian beam)