Details

Type: Bug

Status: Closed

Priority: Major

Resolution: Fixed

Affects Version/s: None

Fix Version/s: 3.3

Labels:None
Description
I think the existing LegendreGauss object fails for certain integrals. An example of failure and a solution that divides segments based on error is provided. Please let me know if I'm not using the LegendreGauss object correctly.
Issue Links
 is a clone of

MATH994 Could we have an infinite integral class?
 Resolved
Activity
 All
 Comments
 Work Log
 History
 Activity
 Transitions
The attached files show how to use AdaptiveQuadrature and how the existing method fails. I've wrapped the IterativeLegendreGaussIntegrator in an InfiniteIntegral object to show a specific instance of a failure of the class. The AdaptiveQuadrature object is more efficient at solving problems (in function evaluation counts) because it selectively increases resolution where the error is high.
This problem is not limited to infinite integrals because the underlying IterativeLegendreGaussIntegrator is integrating in the region [1,1].
The attached solution uses 1st and 2nd order polynomials, but it can be generalized to a higher order polynomial solutions.