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


Post new comment

(Note: commenting is not possible on this snapshot.)

Before posting from this form, please consider whether it would be more appropriate to make an inline comment using the Turn commenting on link near the bottom of the window. (Simply click the link, move the cursor over the article, and click on the piece of text on which you want to comment.)