Abstract
We present a direct, adaptive solver for the Poisson equation which can achieve any prescribed order of accuracy. It is based on a domain decomposition approach using local spectral approximation, as well as potential theory and the fast multipole method. In two space dimensions, the algorithm requires O(NK) work, where N is the number of discretization points and K is the desired order of accuracy.
| Original language | English |
|---|---|
| Pages (from-to) | 415-424 |
| Number of pages | 10 |
| Journal | Journal of Computational Physics |
| Volume | 125 |
| Issue number | 2 |
| DOIs | |
| State | Published - May 1996 |
Bibliographical note
Funding Information:* The authors were supported by the Applied Mathematical Sciences Program of the U.S. Department of Energy under Contract DEFGO288ER25053, by the Office of Naval Research under Contract N00014-91-J-1312, by a NSF Presidential Young Investigator Award to L.G. and by a Packard Foundation Fellowship to L.G.