Abstract:Vandermonde matrices are usually exponentially ill-conditioned and often result in unstable approximations. In this paper, we introduce and analyze the \textit{multivariate Vandermonde with Arnoldi (V+A) method}, which is based on least-squares approximation together with a Stieltjes orthogonalization process, for approximating continuous, multivariate functions on $d$-dimensional irregular domains. The V+A method addresses the ill-conditioning of the Vandermonde approximation by creating a set of discrete orthogonal basis with respect to a discrete measure. The V+A method is simple and general. It relies only on the sample points from the domain and requires no prior knowledge of the domain. In this paper, we first analyze the sample complexity of the V+A approximation. In particular, we show that, for a large class of domains, the V+A method gives a well-conditioned and near-optimal $N$-dimensional least-squares approximation using $M=\mathcal{O}(N^2)$ equispaced sample points or $M=\mathcal{O}(N^2\log N)$ random sample points, independently of $d$. We also give a comprehensive analysis of the error estimates and rate of convergence of the V+A approximation. Based on the multivariate V+A approximation, we propose a new variant of the weighted V+A least-squares algorithm that uses only $M=\mathcal{O}(N\log N)$ sample points to give a near-optimal approximation. Our numerical results confirm that the (weighted) V+A method gives a more accurate approximation than the standard orthogonalization method for high-degree approximation using the Vandermonde matrix.

Robust Padé Approximation via SVD

Algebraic properties of robust Padé approximants

High Order Numerical Methods To Approximate The Singular Value Decomposition

A Power Method for Computing Singular Value Decomposition

$(n D+1)$ Padé approximation

Simultaneous Universal Padé Approximation

Rational Hermite-Padé Approximants vs Padé Approximants. Numerical Results

A mixed precision Jacobi SVD algorithm

PRIMME_SVDS: A High-Performance Preconditioned SVD Solver for Accurate Large-Scale Computations

A Padé Approximate Linearization Algorithm for Solving the Quadratic Eigenvalue Problem with Low‐rank Damping

Close to optimal column approximations with a single SVD

Row-aware Randomized SVD with applications

Theoretical error estimates for computing the matrix logarithm by Padé-type approximants

Convergence and Near-optimal Sampling for Multivariate Function Approximations in Irregular Domains via Vandermonde with Arnoldi

Analysis of the rSVDdpd Algorithm: A Robust Singular Value Decomposition Method using Density Power Divergence

Singular Value Decomposition Approximation via Kronecker Summations for Imaging Applications

Kronecker Product Approximation of Operators in Spectral Norm via Alternating SDP

Efficient Bounds and Estimates for Canonical Angles in Randomized Subspace Approximations

Error Analysis of Padé Iterations for Computing Matrix Invariant Subspaces

Robust SVD Made Easy: A fast and reliable algorithm for large-scale data analysis

Efficient Semidefinite Programming with approximate ADMM