quad_gk 1.0.0

Pure rust numerical integration library based on Gauss Kronrod quadrature rule.
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
42

//! Configuration in numerical integration.

use crate::GK11;

pub struct GKConfig<const N: usize> {
    // required relative tolerance
    pub rel_tol: f64,
    // maximum number of interval divisions
    pub max_interval_count: usize,
    // integral points and coefficients
    pub coef: [(f64, f64, Option<f64>); N],
}

impl<const N: usize> GKConfig<N> {
    // Default configuration.
    pub fn new() -> GKConfig<11> {
        GKConfig {
            rel_tol: 1e-6,
            max_interval_count: 100,
            coef: GK11,
        }
    }

    pub fn rel_tol(self, rel_tol: f64) -> Self {
        Self { rel_tol, ..self }
    }

    pub fn max_interval_count(self, max_interval_count: usize) -> Self {
        Self {
            max_interval_count,
            ..self
        }
    }

    pub fn coef<const M: usize>(self, coef: [(f64, f64, Option<f64>); M]) -> GKConfig<M> {
        GKConfig {
            rel_tol: self.rel_tol,
            max_interval_count: self.max_interval_count,
            coef,
        }
    }
}