gkquad 0.0.4

Numerical integration library for Rust
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

use crate::error::IntegrationResult;
use crate::single::algorithm::*;
use crate::single::common::{Integrand, IntegrationConfig, Range};

/// Automatically select algorithm based on configuration
#[derive(Clone)]
pub struct AUTO {
    qags: QAGS<'static>,
    qagp: QAGP<'static>,
}

impl AUTO {
    #[inline(always)]
    pub fn new() -> Self {
        Self {
            qags: QAGS::new(),
            qagp: QAGP::new(),
        }
    }
}

impl<F: Integrand> Algorithm<F> for AUTO {
    fn integrate(
        &mut self,
        f: &mut F,
        range: &Range,
        config: &IntegrationConfig,
    ) -> IntegrationResult {
        if config.points.is_empty() {
            self.qags.integrate(f, range, config)
        } else {
            self.qagp.integrate(f, range, config)
        }
    }
}

extra_traits!(AUTO);