sequential-integration 1.0.2

Lightweight library for sequential 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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112

mod calculation_step;
pub mod helper_equation_traits;
mod integrators;
pub mod quadrature;
pub mod range_generator;
pub use calculation_step::CalculationStep;
mod calculation_result;
pub use calculation_result::CalculationResult;
mod bounds;
pub use bounds::Bounds;
pub mod function_types;

use fehler::throws;

use crate::errors::Error;

#[throws]
pub fn calculate_single_integral<Q: quadrature::QuadratureSingleIntegral>(
    quadrature: Q,
    first_integral_begin: f64,
    first_integral_end: f64,
) -> f64 {
    let result = integrators::Integrator::integrate::<Q, Q>(
        first_integral_begin,
        first_integral_end,
        quadrature.get_step_size(),
        &quadrature,
        &quadrature,
    )?;

    result
}

#[throws]
pub fn calculate_double_integral<
    Q: quadrature::QuadratureDoubleIntegral,
    F1: Fn(f64) -> f64,
    F2: Fn(f64) -> f64,
>(
    quadrature: Q,
    first_integral_begin: f64,
    first_integral_end: f64,
    second_integral_begin: F1,
    second_integral_end: F2,
) -> f64 {
    let second_integrator = integrators::SecondIntegrator::<Q, Q, F1, F2>::new(
        second_integral_begin,
        second_integral_end,
        quadrature.get_step_size().1,
        &quadrature,
    )?;

    let result =
        integrators::Integrator::integrate::<integrators::SecondIntegrator<Q, Q, F1, F2>, Q>(
            first_integral_begin,
            first_integral_end,
            quadrature.get_step_size().0,
            &second_integrator,
            &quadrature,
        )?;

    result
}

#[throws]
pub fn calculate_triple_integral<
    Q: quadrature::QuadratureTripleIntegral,
    F1: Fn(f64) -> f64,
    F2: Fn(f64) -> f64,
    F3: Fn(f64, f64) -> f64,
    F4: Fn(f64, f64) -> f64,
>(
    quadrature: Q,
    first_integral_begin: f64,
    first_integral_end: f64,
    second_integral_begin: F1,
    second_integral_end: F2,
    third_integral_begin: F3,
    third_integral_end: F4,
) -> f64 {
    let third_integrator = integrators::ThirdIntegrator::<Q, Q, F3, F4>::new(
        third_integral_begin,
        third_integral_end,
        quadrature.get_step_size().2,
        &quadrature,
    )?;

    let second_integrator = integrators::SecondIntegrator::<
        Q,
        integrators::ThirdIntegrator<Q, Q, F3, F4>,
        F1,
        F2,
    >::new(
        second_integral_begin,
        second_integral_end,
        quadrature.get_step_size().1,
        &third_integrator,
    )?;

    let result = integrators::Integrator::integrate::<
        integrators::SecondIntegrator<Q, integrators::ThirdIntegrator<Q, Q, F3, F4>, F1, F2>,
        Q,
    >(
        first_integral_begin,
        first_integral_end,
        quadrature.get_step_size().0,
        &second_integrator,
        &quadrature,
    )?;

    result
}