Solving Laplace problems with the AAA algorithm

Stefano Costa
DOI: https://doi.org/10.48550/arXiv.2001.09439
2020-01-26
Abstract:We present a novel application of the recently developed AAA algorithm to the solution of Laplace 2D problems; an application to conformal mapping is also shown as a particular case. These classes of problems have also been addressed by means of dedicated software implementations based on rational function approximation, with unprecedented speed and accuracy. Still we show how the AAA algorithm, conceived to be a flexible and domain independent method, is capable of addressing these situations within a unified framework.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to use the AAA algorithm to solve two - dimensional Laplace problems and conformal mapping problems. Specifically, the author proposes a new method by using the flexible and adaptive Antoulas - Anderson (AAA) algorithm to handle these mathematical problems and demonstrates its high efficiency and accuracy under complex boundary conditions such as corners. ### Problem Background Laplace problems and conformal mapping problems are widely used in fields such as physics and engineering. Traditional solution methods usually rely on specific software implementations, such as methods based on rational function approximation (such as the "Lightning Laplace Solver" LLS). These methods perform well when dealing with problems such as corners, but they need to perform exponential clustering on sample points to achieve root - exponential convergence. ### Main Contributions of the Paper 1. **General Framework**: Compared with existing specialized software, the method proposed in this paper provides a unified framework that can flexibly deal with various boundary conditions without strict requirements for parameterized boundaries or explicit function definitions. 2. **High Efficiency**: Although not specifically designed for speed optimization, this method can also complete calculations within a few seconds on a modern laptop and has good accuracy. 3. **Flexibility**: The flexibility of the AAA algorithm enables it to return reasonable approximate results when dealing with any reasonably discretized contours, even without exponential clustering sampling near corners. ### Specific Applications - **Laplace Problems**: Include analytical solutions for internal and external regions. - **Conformal Mapping**: Mapping from multiply - connected domains to annular domains, and mapping from L - shaped domains to disks. ### Technical Details The core idea of this method is to decompose the problem into two parts: - **Smooth Part**: Use power series approximation to ensure analyticity within the required region. - **Singular Part**: Use the AAA algorithm for approximation and exclude certain poles as needed to maintain analyticity. In this way, the author shows how to solve different types of Laplace problems and conformal mapping problems within a unified framework while maintaining high computational efficiency and accuracy. ### Summary This paper aims to show the potential of the AAA algorithm in solving Laplace problems and conformal mapping, especially for cases where there are only discrete value sets. Through this method, researchers can obtain high - precision results without relying on complex boundary conditions.