BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas Trogdon

Thomas Trogdon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.PRmath-phmath.MPNumerical Analysisnlin.SImath.APmath.CVmath.FA

Catalog footprint

What is connected

19works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

19 published item(s)

preprint2024arXiv

The ultraspherical rectangular collocation method and its convergence

We develop the ultraspherical rectangular collocation (URC) method, a collocation implementation of the sparse ultraspherical method of Olver \& Townsend for two-point boundary-value problems. The URC method is provably convergent, the implementation is simple and efficient, the convergence proof motivates a preconditioner for iterative methods, and the modification of collocation nodes is straightforward. The convergence theorem applies to all boundary-value problems when the coefficient functions are sufficiently smooth and when the roots of certain ultraspherical polynomials are used as collocation nodes. We also adapt a theorem of Krasnolsel'skii et al.~to our setting to prove convergence for the rectangular collocation method of Driscoll \& Hale for a restricted class of boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

An artificially-damped Fourier method for dispersive evolution equations

Computing solutions to partial differential equations using the fast Fourier transform can lead to unwanted oscillatory behavior. Due to the periodic nature of the discrete Fourier transform, waves that leave the computational domain on one side reappear on the other and for dispersive equations these are typically high-velocity, high-frequency waves. However, the fast Fourier transform is a very efficient numerical tool and it is important to find a way to damp these oscillations so that this transform can still be used. In this paper, we accurately model solutions to four nonlinear partial differential equations on an infinite domain by considering a finite interval and implementing two damping methods outside of that interval: one that solves the heat equation and one that simulates rapid exponential decay. Heat equation-based damping is best suited for small-amplitude, high-frequency oscillations while exponential decay is used to damp traveling waves and high-amplitude oscillations. We demonstrate significant improvements in the runtime of well-studied numerical methods when adding in the damping method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Numerical Unified Transform Method for Initial-boundary Value Problems on the Half-line

We implement the Unified Transform Method of Fokas as a numerical method to solve linear partial differential equations on the half-line. The method computes the solution at any x and t without spatial discretization or time stepping. With the help of contour deformations and oscillatory integration techniques, the method's complexity does not increase for large x,t and the method is more accurate as x,t increase. Our goal is to make no assumptions on the functional form of the initial or boundary functions while maintaining high accuracy in a large region of the (x,t) plane.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universality for the conjugate gradient and MINRES algorithms on sample covariance matrices

We present a probabilistic analysis of two Krylov subspace methods for solving linear systems. We prove a central limit theorem for norms of the residual vectors that are produced by the conjugate gradient and MINRES algorithms when applied to a wide class of sample covariance matrices satisfying some standard moment conditions. The proof involves establishing a four moment theorem for the so-called spectral measure, implying, in particular, universality for the matrix produced by the Lanczos iteration. The central limit theorem then implies an almost-deterministic iteration count for the iterative methods in question.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Inverse Scattering for the Toda Lattice

We present a method to compute the inverse scattering transform (IST) for the famed Toda lattice by solving the associated Riemann--Hilbert (RH) problem numerically. Deformations for the RH problem are incorporated so that the IST can be evaluated in $\mathcal O(1)$ operations for arbitrary points in the $(n,t)$-domain, including short- and long-time regimes. No time-stepping is required to compute the solution because $(n,t)$ appear as parameters in the associated RH problem. The solution of the Toda lattice is computed in long-time asymptotic regions where the asymptotics are not known rigorously.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Smoothed Analysis for the Conjugate Gradient Algorithm

The purpose of this paper is to establish bounds on the rate of convergence of the conjugate gradient algorithm when the underlying matrix is a random positive definite perturbation of a deterministic positive definite matrix. We estimate all finite moments of a natural halting time when the random perturbation is drawn from the Laguerre unitary ensemble in a critical scaling regime explored in Deift et al. (2016). These estimates are used to analyze the expected iteration count in the framework of smoothed analysis, introduced by Spielman and Teng (2001). The rigorous results are compared with numerical calculations in several cases of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An extension of the Derrida-Lebowitz-Speer-Spohn equation

Derrida, Lebowitz, Speer and Spohn have proposed a simplified model to describe the low temperature Glauber dynamics of an anchored Toom interface. We show how the derivation of the Derrida-Lebowitz-Speer-Spohn equation can be prolonged to obtain a new equation, generalizing the models obtained in the paper by these authors. We then investigate its properties from both an analytical and numerical perspective. Specifically, a numerical method is presented to approximate solutions of the prolonged equation. Using this method, we investigate the relationship between the solutions of the prolonged equation and the Tracy--Widom GOE distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gibbs phenomenon for dispersive PDEs on the line

We investigate the Cauchy problem for linear, constant-coefficient evolution PDEs on the real line with discontinuous initial conditions (ICs) in the small-time limit. The small-time behavior of the solution near discontinuities is expressed in terms of universal, computable special functions. We show that the leading-order behavior of the solution of dispersive PDEs near a discontinuity of the ICs is characterized by Gibbs-type oscillations and gives exactly the Wilbraham-Gibbs constant.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the condition number of the critically-scaled Laguerre Unitary Ensemble

We consider the Laguerre Unitary Ensemble (aka, Wishart Ensemble) of sample covariance matrices $A = XX^*$, where $X$ is an $N \times n$ matrix with iid standard complex normal entries. Under the scaling $n = N + \lfloor \sqrt{ 4 c N} \rfloor$, $c > 0$ and $N \rightarrow \infty$, we show that the rescaled fluctuations of the smallest eigenvalue, largest eigenvalue and condition number of the matrices $A$ are all given by the Tracy--Widom distribution ($β= 2$). This scaling is motivated by the study of the solution of the equation $Ax=b$ using the conjugate gradient algorithm, in the case that $A$ and $b$ are random: For such a scaling the fluctuations of the halting time for the algorithm are empirically seen to be universal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast computation of Gauss quadrature nodes and weights on the whole real line

A fast and accurate algorithm for the computation of Gauss-Hermite and generalized Gauss-Hermite quadrature nodes and weights is presented. The algorithm is based on Newton's method with carefully selected initial guesses for the nodes and a fast evaluation scheme for the associated orthogonal polynomial. In the Gauss-Hermite case the initial guesses and evaluation scheme rely on explicit asymptotic formulas. For generalized Gauss-Hermite, the initial guesses are furnished by sampling a certain equilibrium measure and the associated polynomial evaluated via a Riemann-Hilbert reformulation. In both cases the $n$-point quadrature rule is computed in $\mathcal{O}(n)$ operations to an accuracy that is close to machine precision. For sufficiently large $n$, some of the quadrature weights have a value less than the smallest positive normalized floating-point number in double precision and we exploit this fact to achieve a complexity as low as $\mathcal{O}(\sqrt{n})$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the application of GMRES to oscillatory singular integral equations

We present a new method for the numerical solution of singular integral equations on the real axis. The method's value stems from an explicit formula for the Cauchy integral of a complex exponential multiplied by a rational function. Additionally, the inner product of such functions is computed explicitly. With these tools, the GMRES algorithm is applied to both non-oscillatory and oscillatory singular integral equations. Ideas from Fredholm theory and Riemann--Hilbert problems are used to motivate preconditioners for these singular integral equations. A dramatic acceleration in convergence is realized. This presents a strong link between the theory of singular integral equations and the numerical analysis of such equations. Furthermore, this method presents a first step towards a solver for the inverse scattering transform that does not require the deformation of a Riemann--Hilbert problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rational approximation, oscillatory Cauchy integrals and Fourier transforms

We develop the convergence theory for a well-known method for the interpolation of functions on the real axis with rational functions. Precise new error estimates for the interpolant are de- rived using existing theory for trigonometric interpolants. Estimates on the Dirichlet kernel are used to derive new bounds on the associated interpolation projection operator. Error estimates are desired partially due to a recent formula of the author for the Cauchy integral of a specific class of so-called oscillatory rational functions. Thus, error bounds for the approximation of the Fourier transform and Cauchy integral of oscillatory smooth functions are determined. Finally, the behavior of the differentiation operator is discussed. The analysis here can be seen as an extension of that of Weber (1980) and Weideman (1995) in a modified basis used by Olver (2009) that behaves well with respect to function multiplication and differentiation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sampling unitary invariant ensembles

We develop an algorithm for sampling from the unitary invariant random matrix ensembles. The algorithm is based on the representation of their eigenvalues as a determinantal point process whose kernel is given in terms of orthogonal polynomials. Using this algorithm, statistics beyond those known through analysis are calculable through Monte Carlo simulation. Unexpected phenomena are observed in the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universality in Numerical Computations with Random Data. Case Studies

The authors present evidence for universality in numerical computations with random data. Given a (possibly stochastic) numerical algorithm with random input data, the time (or number of iterations) to convergence (within a given tolerance) is a random variable, called the halting time. Two-component universality is observed for the fluctuations of the halting time, i.e., the histogram for the halting times, centered by the sample average and scaled by the sample variance, collapses to a universal curve, independent of the input data distribution, as the dimension increases. Thus, up to two components, the sample average and the sample variance, the statistics for the halting time are universally prescribed. The case studies include six standard numerical algorithms, as well as a model of neural computation and decision making. A link to relevant software is provided in for the reader who would like to do computations of his'r own.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A numerical dressing method for the nonlinear superposition of solutions of the KdV equation

In this paper we present the unification of two existing numerical methods for the construction of solutions of the Korteweg-de Vries (KdV) equation. The first method is used to solve the Cauchy initial-value problem on the line for rapidly decaying initial data. The second method is used to compute finite-genus solutions of the KdV equation. The combination of these numerical methods allows for the computation of exact solutions that are asymptotically (quasi-)periodic finite-gap solutions and are a nonlinear superposition of dispersive, soliton and (quasi-)periodic solutions in the finite (x,t)-plane. Such solutions are referred to as superposition solutions. We compute these solutions accurately for all values of x and t.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Riemann--Hilbert approach to Jacobi operators and Gaussian quadrature

The computation of the entries of Jacobi operators associated with orthogonal polynomials has important applications in numerical analysis. From truncating the operator to form a Jacobi matrix, one can apply the Golub--Welsh algorithm to compute the Gaussian quadrature weights and nodes. Furthermore, the entries of the Jacobi operator are the coefficients in the three-term recurrence relationship for the polynomials. This provides an efficient method for evaluating the orthogonal polynomials. Here, we present an $\mathcal O(N)$ method to compute the first $N$ rows of Jacobi operators from the associated weight. The method exploits the Riemann--Hilbert representation of the polynomials by solving a deformed Riemann--Hilbert problem numerically. We further adapt this computational approach to certain entire weights that are beyond the reach of current asymptotic Riemann--Hilbert techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dispersive and soliton perturbations of finite-genus solutions of the KdV equation: computational results

All solutions of the Korteweg -- de Vries equation that are bounded on the real line are physically relevant, depending on the application area of interest. Usually, both analytical and numerical approaches consider solution profiles that are either spatially localized or (quasi)periodic. In this paper, we discuss a class of solutions that is a nonlinear superposition of these two cases: their asymptotic state for large $|x|$ is (quasi)periodic, but they may contain solitons, with or without dispersive tails. Such scenarios might occur in the case of localized perturbations of previously present sea swell, for instance. Such solutions have been discussed from an analytical point of view only recently. We numerically demonstrate different features of these solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear steepest descent and the numerical solution of Riemann-Hilbert problems

The effective and efficient numerical solution of Riemann-Hilbert problems has been demonstrated in recent work. With the aid of ideas from the method of nonlinear steepest descent for Riemann-Hilbert problems, the resulting numerical methods have been shown numerically to retain accuracy as values of certain parameters become arbitrarily large. Remarkably, this numerical approach does not require knowledge of local parametrices; rather, the deformed contour is scaled near stationary points at a specific rate. The primary aim of this paper is to prove that this observed asymptotic accuracy is indeed achieved. To do so, we first construct a general theoretical framework for the numerical solution of Riemann-Hilbert problems. Second, we demonstrate the precise link between nonlinear steepest descent and the success of numerics in asymptotic regimes. In particular, we prove sufficient conditions for numerical methods to retain accuracy. Finally, we compute solutions to the homogeneous Painlevé II equation and the modified Korteweg-de Vries equations to explicitly demonstrate the practical validity of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical solution of Riemann--Hilbert problems: random matrix theory and orthogonal polynomials

In recent developments, a general approach for solving Riemann--Hilbert problems numerically has been developed. We review this numerical framework, and apply it to the calculation of orthogonal polynomials on the real line. Combining this numerical algorithm with an approach to compute Fredholm determinants, we are able to calculate level densities and gap statistics for general finite-dimensional unitary ensembles. We also include a description of how to compute the Hastings--McLeod solution of the homogeneous Painlevé II equation.

0 citations0 reviewsNo code linkedOpen access