The AAA algorithm for rational approximation

Abstract

We introduce a new algorithm for approximation by rational functions on a real or complex set of points, implementable in 40 lines of Matlab and requiring no user input parameters. Even on a disk or interval the algorithm may outperform existing methods, and on more complicated domains it is especially competitive. The core ideas are (1) representation of the rational approximant in barycentric form with interpolation at certain support points and (2) greedy selection of the support points to avoid exponential instabilities. The name AAA stands for "adaptive Antoulas--Anderson" in honor of the authors who introduced a scheme based on (1). We present the core algorithm with a Matlab code and nine applications and describe variants targeted at problems of different kinds. Comparisons are made with vector fitting, RKFIT, and other existing methods for rational approximation.

Cited by in corpus (90)

• Taming Quantum Noise for Efficient Low Temperature Simulations of Open Quantum Systems
• Efficient low temperature simulations for fermionic reservoirs with the hierarchical equations of motion method: Application to the Anderson impurity model
• Robust analytic continuation of Green's functions via projection, pole estimation, and semidefinite relaxation
• Non-Markovian Quantum Exceptional Points
• Non-Hermitian Pseudomodes for Strongly Coupled Open Quantum Systems: Unravelings, Correlations and Thermodynamics
• GW band structure of monolayer MoS2 using the SternheimerGW method and effect of dielectric environment
• Quasi-Lindblad pseudomode theory for open quantum systems
• Equivalence between a time-fractional and an integer-order gradient flow: The memory effect reflected in the energy
• Polydisperse Streaming Instability II. Methods for solving the linear stability problem
• Analytic continuation from limited noisy Matsubara data
• Estimating Dispersion Curves from Frequency Response Functions via Vector-Fitting
• High-performance finite elements with MFEM
• An Operator Inference Oriented Approach for Mechanical Systems
• Recovering the Potential and Order in One-Dimensional Time-Fractional Diffusion with Unknown Initial Condition and Source
• A fractional PDE model for turbulent velocity fields near solid walls
• Multivariate Rational Approximation
• Summation-by-parts operators for general function spaces
• Computation of 2D Stokes flows via lightning and AAA rational approximation
• Determining kernels in linear viscoelasticity
• Estimating Experimental Dispersion Curves from Steady-State Frequency Response Measurements
• Collapse vs. blow up and global existence in the generalized Constantin-Lax-Majda equation
• Managing Temperature in Open Quantum Systems Strongly Coupled with Structured Environments
• Numerical Evaluation of Mittag-Leffler Functions
• Stochastic pole expansion method
• Current-induced forces in nanosystems: A hierarchical equations of motion approach
• Efficient rational approximation of optical response functions with the AAA algorithm
• Windowed multipole representation of R-matrix cross sections
• NEP: a module for the parallel solution of nonlinear eigenvalue problems in SLEPc
• A reduced conjugate gradient basis method for fractional diffusion
• A Unifying Framework for Interpolatory ${\mathcal{L}}_2$-optimal Reduced-order Modeling
• Minimal pole representation for spectral functions
• Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments
• Orthogonal Rational Approximation of Transfer Functions for High-Frequency Circuits
• Toward a certified greedy Loewner framework with minimal sampling
• Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization
• Structured vector fitting framework for mechanical systems
• Least-squares methods for nonnegative matrix factorization over rational functions
• ${\mathcal{H}}_2$-Optimal Model Reduction Using Projected Nonlinear Least Squares
• Multi-dimensional summation-by-parts operators for general function spaces: Theory and construction
• Nanoptera in nonlinear woodpile chains with zero precompression
• Multivariate approximation by polynomial and generalised rational functions
• A Fast Randomized Algorithm for Computing an Approximate Null Space
• Penguin huddling: a continuum model
• A framework to compute resonances arising from multiple scattering
• Resonance modes in microstructured photonic waveguides: Efficient and accurate computation based on AAA rational approximation
• Summation-by-parts operators for general function spaces: The second derivative
• A generalisation of de la Vallée-Poussin procedure to multivariate approximations
• Nonadiabatic H-Atom Scattering Channels on Ge(111) Elucidated by the Hierarchical Equations of Motion
• Next-generation EDIpack: A Lanczos-based package for quantum impurity models featuring general broken-symmetry phases, flexible bath topologies and multi-platform interoperability
• Barycentric rational function approximation made simple: A fast analytic continuation method for Matsubara Green's functions
• Efficient Pseudomode Representation and Complexity of Quantum Impurity Models
• About the performance of perturbative treatments of the spin-boson dynamics within the hierarchical equations of motion approach
• Rational kernel-based interpolation for complex-valued frequency response functions
• Burgers' equation in the complex plane
• Uncovering hidden resonances in non-Hermitian systems with scattering thresholds
• Demonstration of Metaplectic Geometrical Optics for Reduced Modeling of Plasma Waves
• Structured barycentric forms for interpolation-based data-driven reduced modeling of second-order systems
• Adaptive approximation of nonlinear eigenproblems by minimal rational interpolation
• Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-Like Calculations
• Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description
• Fully-Adaptive and Semi-Adaptive Frequency Sweep Algorithm Exploiting Loewner-State Model for EM Simulation of Multiport Systems
• Nonlinear response theory for lossy superconducting quantum circuits
• Diffraction by a set of collinear cracks on a square lattice: an iterative Wiener-Hopf method
• Adaptive choice of near-optimal expansion points for interpolation-based structure-preserving model reduction
• Physically interpretable approximations of many-body spectral functions
• Stokes flows in a two-dimensional bifurcation
• Evaluation of resonances: adaptivity and AAA rational approximation of randomly scalarized boundary integral resolvents
• Regarding the extension of metaplectic geometrical optics to modelling evanescent waves in ray-tracing codes
• Numerical simulation of dilute polymeric fluids with memory effects in the turbulent flow regime
• Constrained Random Phase Approximation: the spectral method
• Unitarity of some barycentric rational approximants
• Structured model order reduction for vibro-acoustic problems using interpolation and balancing methods
• Extended Non-Markovian Stochastic Schrödinger Equation with Complex Frequency Modes for General Basis Functions
• Approximating a flexible beam model in the Loewner framework
• Entanglement dynamics and performance of two-qubit gates for superconducting qubits under non-Markovian effects
• Approximating maps into manifolds with lower curvature bounds
• Purified pseudomode model for nonlinear system-bath interactions
• Surrogate modeling of resonant behavior in scattering problems through adaptive rational approximation and sketching
• Exponential Asymptotics using Numerical Rational Approximation in Linear Differential Equations
• Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting
• Testing bath correlation functions for open quantum dynamics simulations
• Making Non-Markovian master equations accessible with approximate environments
• Walk based Laplacians for Modeling Diffusion on Complex Networks
• Exact description of fermionic reservoirs via purified damped ancillary fermions
• Robust Bilinear-Noise-Optimal Control for Gravitational-Wave Detectors: A Mixed LQG/$H_{\infty }$ Approach
• Barycentric rational approximation for learning the index of a dynamical system from limited data
• Residues in Partial Fraction Decomposition Applied to Pole Sensitivity Analysis and Root Locus Construction
• A quasi-local inhomogeneous dielectric tensor for arbitrary distribution functions
• Computing Stokes flows in periodic channels via rational approximation
• Rational function approximation with normalized positive denominators