quad-rs 0.5.0

Adaptive Gauss-Kronrod Integration in Rust
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

use super::policy::RangeSplitError;
use crate::storage::SegmentHeapError;

#[derive(thiserror::Error, Debug, PartialEq)]
/// Error type for integrator
pub enum IntegratorError<T, E> {
    #[error("exceeded the maximum number of function evaluations before convergence")]
    ExceededMaxFunctionEvaluations,
    #[error("non-finite error estimate")]
    NonFiniteErrorEstimate,
    #[error("no segments found")]
    NoSegments,
    #[error("empty integration segment")]
    EmptySegment,
    #[error("non-finite integrand at {point:?}")]
    NonFiniteIntegrand { point: T },
    #[error("possible singularity at {singularity:?}")]
    PossibleSingularity { singularity: T },
    #[error("refined contour piece smaller than minimum")]
    PieceTooSmall,
    #[error("error in user code")]
    User(E),
}

impl<T, E> From<SegmentHeapError> for IntegratorError<T, E> {
    fn from(value: SegmentHeapError) -> Self {
        match value {
            SegmentHeapError::NonFiniteErrorEstimate => Self::NonFiniteErrorEstimate,
        }
    }
}

impl<T: nalgebra::ComplexField, E> From<RangeSplitError<T>> for IntegratorError<T, E> {
    fn from(value: RangeSplitError<T>) -> Self {
        match value {
            RangeSplitError::SegmentTooSmall { center, .. } => Self::PossibleSingularity {
                singularity: center,
            },
        }
    }
}