An FFT Based Fast Poisson Solver on Spherical Shells
Yin-Liang Huang 1*, Jian-Guo Liu 2, Wei-Cheng Wang 31 Department of Applied Mathematics, National University of Tainan, Tainan 70005, Taiwan.
2 Department of Physics and Department of Mathematics, Duke University, Durham, NC 27708, USA.
3 Department of Mathematics, National Tsing Hua University, Hsinchu 30013, Taiwan.
Received 6 May 2009; Accepted (in revised version) 8 June 2009
Available online 17 September 2010
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\log N)$ operations. Numerical examples have confirmed the accuracy and robustness of the new scheme.AMS subject classifications: 35Q86, 65N06, 65N15, 65N22, 65T50
Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Poisson equation, spherical coordinate, FFT, spectral-finite difference method, fast diagonalization, high order accuracy, error estimate, trapezoidal rule, Euler-Maclaurin formula, Bernoulli numbers.
Email: email@example.com (Y.-L. Huang), firstname.lastname@example.org (J.-G. Liu), email@example.com (W.-C. Wang)