quad 0.1.2

Quadrature integration
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66

use crate::constants::{HeapItem, Myf64};
use std::collections::{BinaryHeap, HashMap};

#[derive(Debug, Clone)]
pub struct QagIntegrationResult {
    pub result: Vec<f64>,
    pub abserr: f64,
    pub more_info: Option<MoreInfo>,
}

impl QagIntegrationResult {
    pub fn new_more_info(
        result: Vec<f64>,
        abserr: f64,
        neval: i32,
        last: usize,
        hash: HashMap<(Myf64, Myf64), Vec<f64>>,
        heap: BinaryHeap<HeapItem>,
    ) -> Self {
        Self {
            result,
            abserr,
            more_info: Some(MoreInfo::new(neval, last, hash, heap)),
        }
    }

    pub fn new(result: Vec<f64>, abserr: f64) -> Self {
        Self {
            result,
            abserr,
            more_info: None,
        }
    }

    pub fn new_error() -> Self {
        Self {
            result: vec![0.0],
            abserr: 0.0,
            more_info: None,
        }
    }
}

#[derive(Debug, Clone)]
pub struct MoreInfo {
    pub neval: i32,
    pub last: usize,
    pub hash: HashMap<(Myf64, Myf64), Vec<f64>>,
    pub heap: BinaryHeap<HeapItem>,
}

impl MoreInfo {
    pub fn new(
        neval: i32,
        last: usize,
        hash: HashMap<(Myf64, Myf64), Vec<f64>>,
        heap: BinaryHeap<HeapItem>,
    ) -> Self {
        Self {
            neval,
            last,
            hash,
            heap,
        }
    }
}