zbMATH — the first resource for mathematics

Why Gaussian quadrature in the complex plane? (English) Zbl 0977.41013
Numer. Algorithms 26, No. 3, 251-280 (2001); addendum ibid. 27, No. 2, 215-217 (2001).
Authors’ abstract: This paper synthesizes formally orthogonal polynomials, Gaussian quadrature in the complex plane and the bi-conjugate gradient method together with an application. Classical Gaussian quadrature approximates an integral over (a region of) the real line. We present an extension of Gaussian quadrature over an arc in the complex plane, which we call complex Gaussian quadrature. Since there has not been any particular interest in the numerical evaluation of integrals over the long history of complex function theory, complex Gaussian quadrature is in need of motivation. Gaussian quadrature in the complex plane yields approximations of certain sums connected with the bi-conjugate gradient method. The scattering amplitude \(c^TA^{-1 }b\) is an example where \(A\) is a discretization of a differential integral operator corresponding to the scattering problem and \(b\) and \(c\) are given vectors. The usual method to estimate this is to use \(c^Tx^{(k)}\). A result of Warnick is that this is identically equal to the complex Gaussian quadrature estimate of \(1/\lambda\). Complex Gaussian quadrature thereby replaces this particular inner product in the estimate of the scattering amplitude.

41A55 Approximate quadratures
65D32 Numerical quadrature and cubature formulas
65F10 Iterative numerical methods for linear systems
Full Text: DOI