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
A Direct Adaptive Poisson Solver of Arbitrary Order Accuracy
โ Scribed by Leslie Greengard; June-Yub Lee
- Publisher
- Elsevier Science
- Year
- 1996
- Tongue
- English
- Weight
- 351 KB
- Volume
- 125
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
โฆ Synopsis
Almost all currently available methods are based on iterative techniques using multigrid [7,23], domain decom-We present a direct, adaptive solver for the Poisson equation which can achieve any prescribed order of accuracy. It is based on position [11], or some other preconditioning strategy. Una domain decomposition approach using local spectral approximafortunately, while multilevel iterations can achieve optimal tion, as well as potential theory and the fast multipole method. In efficiency in theory, they require an appropriate hierarchy two space dimensions, the algorithm requires O(NK ) work, where of coarse grids which are not provided in many practical N is the number of discretization points and K is the desired order situations. There has, however, been significant progress of accuracy.
๐ SIMILAR VOLUMES
We present a fourth order numerical solution method for the singular Neumann boundary problem of Poisson equations. Such problems arise in the solution process of incompressible Navier-Stokes equations and in the time-harmonic wave propagation in the frequence space with the zero wavenumber. The equ
The approximate solution of the differential equation d'+d? + Q'(z)+ = 0 by a general modification of certain phase-integral approximations of arbitrary order is considered. A consistent modification of these higher-order phase-integral approximations is derived on the assumption that one has found
Following the approach of Jones, Low, and Young, a generalized O(2,l) expansion is developed for amplitudes that have a power bounded growth asymptotically. The expansion, set up in an 0(1, 1) basis, holds in a new kinematical region, where all the incoming and outgoing clusters have space-like 0(2,