pub fn integrate(
f: impl FnMut(Scalar) -> GeopResult<Vec<Scalar>>,
breaks: &[Scalar],
components: usize,
quadrature: &Quadrature,
) -> GeopResult<Integral>Expand description
The integral of the components-valued f from breaks[0] to the last
of breaks, each consecutive pair a panel to start with — where the
integrand is not smooth, say, at the knots of a spline.
Global adaptive: the panel whose error is largest, relative to its
component’s scale, is halved until every component has converged or
the budget is spent (see Quadrature). Where to halve a panel is a
free choice, so its midpoint is sharpened.