English

A $p$-adaptive treecode solution of the Poisson equation in the general domain

Numerical Analysis 2024-12-31 v1 Numerical Analysis

Abstract

Raising the order of the multipole expansion is a feasible approach for improving the accuracy of the treecode algorithm. However, a uniform order for the expansion would result in the inefficiency of the implementation, especially when the kernel function is singular. In this paper, a pp-adaptive treecode algorithm is designed to resolve the efficiency issue for problems defined on a general domain. Such a pp-adaptive implementation is realized through i). conducting a systematical error analysis for the treecode algorithm, ii). designing a strategy for a non-uniform distribution of the order of multipole expansion towards a given error tolerance, and iii). employing a hierarchy geometry tree structure for coding the algorithm. The proposed pp-adaptive treecode algorithm is validated by a number of numerical experiments, from which the desired performance is observed successfully, i.e., the computational complexity is reduced dramatically compared with the uniform order case, making our algorithm a competitive one for bottleneck problems such as the demagnetizing field calculation in computational micromagnetics.

Keywords

Cite

@article{arxiv.2412.20434,
  title  = {A $p$-adaptive treecode solution of the Poisson equation in the general domain},
  author = {Zixuan Cui and Lei Yang},
  journal= {arXiv preprint arXiv:2412.20434},
  year   = {2024}
}
R2 v1 2026-06-28T20:51:04.812Z