The numerical manifold method for harmonic wave propagation in unbounded domains. (December 2022)