Researcher profile

Lloyd N. Trefethen

Lloyd N. Trefethen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
7works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Spectacularly large expansion coefficients in Müntz's theorem

Müntz's theorem asserts, for example, that the even powers $1, x^2, x^4,\dots$ are dense in $C([0,1])$. We show that the associated expansions are so inefficient as to have no conceivable relevance to any actual computation. For example, approximating $f(x)=x$ to accuracy $\varepsilon = 10^{-6}$ in this basis requires powers larger than $x^{280{,}000}$ and coefficients larger than $10^{107{,}000}$. We present a theorem establishing exponential growth of coefficients with respect to $1/\varepsilon$.

preprint2020arXiv

Exponential node clustering at singularities for rational approximation, quadrature, and PDEs

Rational approximations of functions with singularities can converge at a root-exponential rate if the poles are exponentially clustered. We begin by reviewing this effect in minimax, least-squares, and AAA approximations on intervals and complex domains, conformal mapping, and the numerical solution of Laplace, Helmholtz, and biharmonic equations by the "lightning" method. Extensive and wide-ranging numerical experiments are involved. We then present further experiments showing that in all of these applications, it is advantageous to use exponential clustering whose density on a logarithmic scale is not uniform but tapers off linearly to zero near the singularity. We give a theoretical explanation of the tapering effect based on the Hermite contour integral and potential theory, showing that tapering doubles the rate of convergence. Finally we show that related mathematics applies to the relationship between exponential (not tapered) and doubly exponential (tapered) quadrature formulas. Here it is the Gauss--Takahasi--Mori contour integral that comes into play.

preprint2016arXiv

Multivariate polynomial approximation in the hypercube

A theorem is proved concerning approximation of analytic functions by multivariate polynomials in the $s$-dimensional hypercube. The geometric convergence rate is determined not by the usual notion of degree of a multivariate polynomial, but by the {\it Euclidean degree,} defined in terms of the 2-norm rather than the 1-norm of the exponent vector $\bf k$ of a monomial $x_1^{k_1}\cdots \kern .8pt x_s^{k_s}$.

preprint2016arXiv

Trigonometric Interpolation and Quadrature in Perturbed Points

The trigonometric interpolants to a periodic function $f$ in equispaced points converge if $f$ is Dini-continuous, and the associated quadrature formula, the trapezoidal rule, converges if $f$ is continuous. What if the points are perturbed? With equispaced grid spacing $h$, let each point be perturbed by an arbitrary amount $\le αh$, where $α\in [\kern .5pt 0,1/2)$ is a fixed constant. The Kadec 1/4 theorem of sampling theory suggests there may be be trouble for $α\ge 1/4$. We show that convergence of both the interpolants and the quadrature estimates is guaranteed for all $α<1/2$ if $f$ is twice continuously differentiable, with the convergence rate depending on the smoothness of $f$. More precisely it is enough for $f$ to have $4α$ derivatives in a certain sense, and we conjecture that $2α$ derivatives is enough. Connections with the Fejér--Kalmár theorem are discussed.

preprint2015arXiv

Chopping a Chebyshev Series

Chebfun and related software projects for numerical computing with functions are based on the idea that at each step of a computation, a function $f(x)$ defined on an interval $[a,b]$ is &#34;rounded&#34; to a prescribed precision by constructing a Chebyshev series and chopping it at an appropriate point. Designing a chopping algorithm with the right properties proves to be a surprisingly complex and interesting problem. We describe the chopping algorithm introduced in Chebfun Version 5.3 in 2015 after many years of discussion and the considerations that led to this design.