39,879 research outputs found
A Lanczos Method for Approximating Composite Functions
We seek to approximate a composite function h(x) = g(f(x)) with a global
polynomial. The standard approach chooses points x in the domain of f and
computes h(x) at each point, which requires an evaluation of f and an
evaluation of g. We present a Lanczos-based procedure that implicitly
approximates g with a polynomial of f. By constructing a quadrature rule for
the density function of f, we can approximate h(x) using many fewer evaluations
of g. The savings is particularly dramatic when g is much more expensive than f
or the dimension of x is large. We demonstrate this procedure with two
numerical examples: (i) an exponential function composed with a rational
function and (ii) a Navier-Stokes model of fluid flow with a scalar input
parameter that depends on multiple physical quantities
Solving polynomial eigenvalue problems by means of the Ehrlich-Aberth method
Given the matrix polynomial , we
consider the associated polynomial eigenvalue problem. This problem, viewed in
terms of computing the roots of the scalar polynomial , is treated
in polynomial form rather than in matrix form by means of the Ehrlich-Aberth
iteration. The main computational issues are discussed, namely, the choice of
the starting approximations needed to start the Ehrlich-Aberth iteration, the
computation of the Newton correction, the halting criterion, and the treatment
of eigenvalues at infinity. We arrive at an effective implementation which
provides more accurate approximations to the eigenvalues with respect to the
methods based on the QZ algorithm. The case of polynomials having special
structures, like palindromic, Hamiltonian, symplectic, etc., where the
eigenvalues have special symmetries in the complex plane, is considered. A
general way to adapt the Ehrlich-Aberth iteration to structured matrix
polynomial is introduced. Numerical experiments which confirm the effectiveness
of this approach are reported.Comment: Submitted to Linear Algebra App
Detailed ultraviolet asymptotics for AdS scalar field perturbations
We present a range of methods suitable for accurate evaluation of the leading
asymptotics for integrals of products of Jacobi polynomials in limits when the
degrees of some or all polynomials inside the integral become large. The
structures in question have recently emerged in the context of effective
descriptions of small amplitude perturbations in anti-de Sitter (AdS)
spacetime. The limit of high degree polynomials corresponds in this situation
to effective interactions involving extreme short-wavelength modes, whose
dynamics is crucial for the turbulent instabilities that determine the ultimate
fate of small AdS perturbations. We explicitly apply the relevant asymptotic
techniques to the case of a self-interacting probe scalar field in AdS and
extract a detailed form of the leading large degree behavior, including closed
form analytic expressions for the numerical coefficients appearing in the
asymptotics.Comment: v2: 19 pages, expanded version accepted to JHE
A numerical method for oscillatory integrals with coalescing saddle points
The value of a highly oscillatory integral is typically determined
asymptotically by the behaviour of the integrand near a small number of
critical points. These include the endpoints of the integration domain and the
so-called stationary points or saddle points -- roots of the derivative of the
phase of the integrand -- where the integrand is locally non-oscillatory.
Modern methods for highly oscillatory quadrature exhibit numerical issues when
two such saddle points coalesce. On the other hand, integrals with coalescing
saddle points are a classical topic in asymptotic analysis, where they give
rise to uniform asymptotic expansions in terms of the Airy function. In this
paper we construct Gaussian quadrature rules that remain uniformly accurate
when two saddle points coalesce. These rules are based on orthogonal
polynomials in the complex plane. We analyze these polynomials, prove their
existence for even degrees, and describe an accurate and efficient numerical
scheme for the evaluation of oscillatory integrals with coalescing saddle
points
- …