a repository of mathematical know-how

Revision of Mapped Quadrature from Mon, 20/07/2009 - 07:15

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 [\maths] 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 Donaldson and Elliott

Residues and Steepest Descents

Riemann-Hilbert Problems

Percy Deift

Example 1

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

Example 2



Example 3

Trapezoidal Rule - Sinc, Double Exponential

Stenger, Mori

Example 4



Example 5

Rational Basis Functions


Example 6

Elliott and Johnson


Gauss Legendre

General discussion