1
$\begingroup$

For integration of polynomials, the number of qradrature points can be determined easily according to the $(2N-1)$ rule. While in most cases, it's not necessary to use Gaussian quadrature for simple polynomials since they can be calculated analytically, what we encounter in reality are generally complex integrands, such as rational functions.

In this case, how can we determin an appropriate number of quadrature points so that we can get a somewhat acceptible solution, and at the same time will not lead to a waste of computational efforts?

Improve this question
$\endgroup$
2
  • 3
    $\begingroup$ Possible duplicate of Method selection for numeric quadrature $\endgroup$ – GertVdE Nov 14 '15 at 17:30
  • 1
    $\begingroup$ Quadrature is used for integrating simple polynomials in finite element codes since it is easier than symbolically evaluating them. $\endgroup$ – James Nov 14 '15 at 19:53
4
$\begingroup$

There is a formula for the error term in Gaussian quadrature, so if you bound the derivatives of your integrand, you can bound the error.

If you can't bound the derivatives of your integrand you can use a Gauss- Kronrod rule in which you first use n points to estimate the integral and then reuse those points with n+1 additional points to get a better estimate. The difference of the two estimates is an estimate of the error.

This technique of reusing points in a higher order approximation is also used in Clenshaw-Curtis integration.

This is standard material in introductory courses in numerical analysis that you can find in many textbooks.

Improve this answer
$\endgroup$
0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.