A fast 3D Poisson solver of arbitrary order accuracy

被引:43
作者
Braverman, E [1 ]
Israeli, M
Averbuch, A
Vozovoi, L
机构
[1] Technion Israel Inst Technol, Dept Comp Sci, IL-32000 Haifa, Israel
[2] Tel Aviv Univ, Sch Math Sci, IL-69978 Tel Aviv, Israel
关键词
3D Poisson solver for Dirichlet problem; Fourier method; corner and edge singularities;
D O I
10.1006/jcph.1998.6001
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present a direct solver for the Poisson and Laplace equations in a 3D rectangular box. The method is based on the application of the discrete Fourier transform accompanied by a subtraction technique which allows reducing the errors associated with the Gibbs phenomenon and achieving any prescribed rate of convergence. The algorithm requires O (N-3 log N) operations, where N is the number of grid points in each direction. We show that our approach allows accurate treatment of singular cases which arise when the boundary function is discontinuous or incompatible with the differential equation. (C) 1998 Academic Press.
引用
收藏
页码:109 / 136
页数:28
相关论文
共 20 条
[1]  
Aharoni G., 1993, Journal of Mathematical Imaging and Vision, V3, P7, DOI 10.1007/BF01248401
[2]  
[Anonymous], 1989, DOMAIN DECOMPOSITION
[3]   On a fast direct elliptic solver by a modified Fourier method [J].
Averbuch, A ;
Vozovoi, L ;
Israeli, M .
NUMERICAL ALGORITHMS, 1997, 15 (3-4) :287-313
[4]  
AVERBUCH A, 1998, S NEAM WORKSH SIGN I, P341
[5]  
AVERBUCH A, IN PRESS SIAM J SCI
[6]  
AVERBUCH A, UNPUB ADAPTIVE SOLUT
[7]  
BRAVERMAN E, UNPUB PARALLEL ADAPT
[8]  
BRAVERMAN E, UNPUB FAST SPECTRAL
[9]  
Canuto C., 1989, SPECTRAL METHODS FLU
[10]  
Gottlieb D., 1977, NUMERICAL ANAL SPECT