摘要
We present a fast Poisson solver on spherical shells.With a special change of variable,the radial part of the Laplacian transforms to a constant coefficient differential operator.As a result,the Fast Fourier Transform can be applied to solve the Poisson equation with O(N^(3) logN)operations.Numerical examples have confirmed the accuracy and robustness of the new scheme.
基金
The research of Liu was supported by the NSF grant DMS 10-11738
The research of Wang was supported by National Science Council of Taiwan under grant 97-2115-M-007-005
In addition,this work is also supported in part by National Center for Theoretical Sciences of Taiwan.
