The detailed explanation is in the attached file.
The problem is that I get the following error while calculating Laplace operator RuntimeWarning: overflow encountered in multiply Laplace= (k * k + l * l).real
My question is how to avoid the error or is there another/better way to calculate wavenumbers on a shifted gird points.