A simple spectral algorithm for solving large-scale Poisson equation in 2D

摘要

We show that it is possible with easy-to-program algorithms to reach spatial resolutions of the order of 10~8 grid points for computing the eclectic potential on 2D periodic lattices, such as the Si (111)7 * 7 surface. We have used a spectral Fourier technique and parallelized FFTs with OPEN_Mp on SGI machines. This method can be easily extended to 3D.