Log In Sign Up

AAA-least squares rational approximation and solution of Laplace problems

by   Stefano Costa, et al.

A two-step method for solving planar Laplace problems via rational approximation is introduced. First complex rational approximations to the boundary data are determined by AAA approximation, either globally or locally near each corner or other singularity. The poles of these approximations outside the problem domain are then collected and used for a global least-squares fit to the solution. Typical problems are solved in a second of laptop time to 8-digit accuracy, all the way up to the corners, and the conjugate harmonic function is also provided. The AAA-least squares combination also offers a new method for avoiding spurious poles in other rational approximation problems, and for greatly speeding them up in cases with many singularities. As a special case, AAA-LS approximation leads to a powerful method for computing the Hilbert transform or Dirichlet-to-Neumann map.


page 1

page 2

page 3

page 4


Solving Laplace problems with corner singularities via rational functions

A new method is introduced for solving Laplace problems on 2D regions wi...

Solving Laplace problems with the AAA algorithm

We present a novel application of the recently developed AAA algorithm t...

Numerical conformal mapping with rational functions

New algorithms are presented for numerical conformal mapping based on ra...

An algorithm for real and complex rational minimax approximation

Rational minimax approximation of real functions on real intervals is an...

Zolotarev's fifth and sixth problems

In an influential 1877 paper, Zolotarev asked and answered four question...

Nonlinear approximation of functions based on non-negative least squares solver

In computational practice, most attention is paid to rational approximat...

Modifying iterated Laplace approximations

In this paper, several modifications are introduced to the functional ap...