a repository of mathematical know-how

Mapped Quadrature

Quick description

Many quadrature rules are derived from others by first making a change-of-variable in the integral and then applying an existing quadrature rule.

 I=\int_a^b f(x) \, dx =\int_{g(a)}^{g(b)} f(g(u)) g'(u) \, du \approx Q = \sum w_i f(g(u_i))

Quadrature error can be represented by a contour integral over a curve enclosing the interval. This error can be estimated using residues around poles, integrals along branch cuts, or saddle points.

 I-Q_n = \frac{1}{2\pi i} \int_C G(w) f(g(w)) g'(w) dw

where C encloses (g(a),g(b)) and G(w) has poles at the quadrature points u_i.

My idea is to unify the existing cases by defining a Riemann surface by the change-of-variable and then the error is represented by an integral over the Riemann surface.


  • Riemann Surfaces

  • Integral Representation of Error - eg Donaldson and Elliott

  • Residues and Steepest Descents

  • Riemann-Hilbert Problems - Percy Deift

Example 1

Simple x=e^u, x=\sqrt{u}

Example 3

Trapezoidal Rule - Sinc Stenger

Example 4

Trapezoidal Rule - Double Exponential Mori

Example 6

Rational Basis Functions


Example 7

Elliott and Johnson


Gauss Legendre

General discussion