1,801 research outputs found
Conformal Maps to Multiply-Slit Domains and Applications
By exploiting conformal maps to vertically slit regions in the complex plane, a recently developed rational spectral method [Tee and Trefethen, 2006] is able to solve PDEs with interior layer-like behaviour using significantly fewer collocation points than traditional spectral methods. The conformal maps are chosen to 'enlarge the region of analyticity' in the solution: an idea which can be extended to other numerical methods based upon global polynomial interpolation. Here we show how such maps can be rapidly computed in both periodic and nonperiodic geometries, and apply them to some challenging differential equations
On the computation of Gaussian quadrature rules for Chebyshev sets of linearly independent functions
We consider the computation of quadrature rules that are exact for a
Chebyshev set of linearly independent functions on an interval . A
general theory of Chebyshev sets guarantees the existence of rules with a
Gaussian property, in the sense that basis functions can be integrated
exactly with just points and weights. Moreover, all weights are positive
and the points lie inside the interval . However, the points are not the
roots of an orthogonal polynomial or any other known special function as in the
case of regular Gaussian quadrature. The rules are characterized by a nonlinear
system of equations, and earlier numerical methods have mostly focused on
finding suitable starting values for a Newton iteration to solve this system.
In this paper we describe an alternative scheme that is robust and generally
applicable for so-called complete Chebyshev sets. These are ordered Chebyshev
sets where the first elements also form a Chebyshev set for each . The
points of the quadrature rule are computed one by one, increasing exactness of
the rule in each step. Each step reduces to finding the unique root of a
univariate and monotonic function. As such, the scheme of this paper is
guaranteed to succeed. The quadrature rules are of interest for integrals with
non-smooth integrands that are not well approximated by polynomials
How fast do radial basis function interpolants of analytic functions converge?
The question in the title is answered using tools of potential theory. Convergence and divergence rates of interpolants of analytic functions on the unit interval are analyzed. The starting point is a complex variable contour integral formula for the remainder in RBF interpolation. We study a generalized Runge phenomenon and explore how the location of centers and affects convergence. Special attention is given to Gaussian and inverse quadratic radial functions, but some of the results can be extended to other smooth basis functions. Among other things, we prove that, under mild conditions, inverse quadratic RBF interpolants of functions that are analytic inside the strip , where is the shape parameter, converge exponentially
A fast and well-conditioned spectral method for singular integral equations
We develop a spectral method for solving univariate singular integral
equations over unions of intervals by utilizing Chebyshev and ultraspherical
polynomials to reformulate the equations as almost-banded infinite-dimensional
systems. This is accomplished by utilizing low rank approximations for sparse
representations of the bivariate kernels. The resulting system can be solved in
operations using an adaptive QR factorization, where is
the bandwidth and is the optimal number of unknowns needed to resolve the
true solution. The complexity is reduced to operations by
pre-caching the QR factorization when the same operator is used for multiple
right-hand sides. Stability is proved by showing that the resulting linear
operator can be diagonally preconditioned to be a compact perturbation of the
identity. Applications considered include the Faraday cage, and acoustic
scattering for the Helmholtz and gravity Helmholtz equations, including
spectrally accurate numerical evaluation of the far- and near-field solution.
The Julia software package SingularIntegralEquations.jl implements our method
with a convenient, user-friendly interface
The exponentially convergent trapezoidal rule
It is well known that the trapezoidal rule converges geometrically when applied to analytic functions on periodic intervals or the real line. The mathematics and history of this phenomenon are reviewed and it is shown that far from being a curiosity, it is linked with computational methods all across scientific computing, including algorithms related to inverse Laplace transforms, special functions, complex analysis, rational approximation, integral equations, and the computation of functions and eigenvalues of matrices and operators
Regularity Theory and Superalgebraic Solvers for Wire Antenna Problems
We consider the problem of evaluating the current distribution that is induced on a straight wire antenna by a time-harmonic incident electromagnetic field. The scope of this paper is twofold. One of its main contributions is a regularity proof for a straight wire occupying the interval . In particular, for a smooth time-harmonic incident field this theorem implies that , where is an infinitely differentiable function—the previous state of the art in this regard placed in the Sobolev space , . The second focus of this work is on numerics: we present three superalgebraically convergent algorithms for the solution of wire problems, two based on Hallén's integral equation and one based on the Pocklington integrodifferential equation. Both our proof and our algorithms are based on two main elements: (1) a new decomposition of the kernel of the form , where and are analytic functions on the real line; and (2) removal of the end-point square root singularities by means of a coordinate transformation. The Hallén- and Pocklington-based algorithms we propose converge superalgebraically: faster than and for any positive integer , where and are the numbers of unknowns and the number of integration points required for construction of the discretized operator, respectively. In previous studies, at most the leading-order contribution to the logarithmic singular term was extracted from the kernel and treated analytically, the higher-order singular derivatives were left untreated, and the resulting integration methods for the kernel exhibit convergence at best. A rather comprehensive set of tests we consider shows that, in many cases, to achieve a given accuracy, the numbers of unknowns required by our codes are up to a factor of five times smaller than those required by the best solvers previously available; the required number of integration points, in turn, can be several orders of magnitude smaller than those required in previous methods. In particular, four-digit solutions were found in computational times of the order of four seconds and, in most cases, of the order of a fraction of a second on a contemporary personal computer; much higher accuracies result in very small additional computing times
- …