rint 0.1.0

A pure Rust library for the numerical integration of real or complex valued functions of real variables in multiple dimensions.
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

// Test the basic Gauss-Kronrod integration rules with a smooth positive function.
// Ported from gsl-2.6/integration/test.c

use num_complex::Complex;
use num_complex::ComplexFloat;
use rint::quadrature::Basic;
use rint::quadrature::Rule;
use rint::Limits;

mod util;

basic_test_complex_equal! {
        name: basic_gauss_kronrod_smooth_function_equal_re_im_15_point,
        function: util::ComplexF1F1,
        alpha: f64 => 2.6,
        rule: gk15,
        lower: 0.0,
        upper: 1.0,
        exp_result: 7.716049357767090777E-02,
        exp_error: 2.990224871000550874E-06,
        abs_tolerance: 1.0e-15,
        rel_tolerance: 1.0e-7,
        "ComplexF1F1 alpha1 = alpha2 Basic 15-point"
}

basic_test_complex_equal! {
        name: basic_gauss_kronrod_smooth_function_equal_re_im_21_point,
        function: util::ComplexF1F1,
        alpha: f64 => 2.6,
        rule: gk21,
        lower: 0.0,
        upper: 1.0,
        exp_result: 7.716049379303084599E-02,
        exp_error: 9.424302194248481445E-08,
        abs_tolerance: 1.0e-15,
        rel_tolerance: 1.0e-7,
        "ComplexF1F1 alpha1 = alpha2 Basic 21-point"
}

basic_test_complex_equal! {
        name: basic_gauss_kronrod_smooth_function_equal_re_im_31_point,
        function: util::ComplexF1F1,
        alpha: f64 => 2.6,
        rule: gk31,
        lower: 0.0,
        upper: 1.0,
        exp_result: 7.716049382494900855E-02,
        exp_error: 1.713503193600029893E-09,
        abs_tolerance: 1.0e-15,
        rel_tolerance: 1.0e-7,
        "ComplexF1F1 alpha1 = alpha2 Basic 31-point"
}

basic_test_complex_equal! {
        name: basic_gauss_kronrod_smooth_function_equal_re_im_41_point,
        function: util::ComplexF1F1,
        alpha: f64 => 2.6,
        rule: gk41,
        lower: 0.0,
        upper: 1.0,
        exp_result: 7.716049382681375302E-02,
        exp_error: 9.576386660975511224E-11,
        abs_tolerance: 1.0e-15,
        rel_tolerance: 1.0e-7,
        "ComplexF1F1 alpha1 = alpha2 Basic 41-point"
}

basic_test_complex_equal! {
        name: basic_gauss_kronrod_smooth_function_equal_re_im_51_point,
        function: util::ComplexF1F1,
        alpha: f64 => 2.6,
        rule: gk51,
        lower: 0.0,
        upper: 1.0,
        exp_result: 7.716049382708510540E-02,
        exp_error: 1.002079980317363772E-11,
        abs_tolerance: 1.0e-15,
        rel_tolerance: 1.0e-7,
        "ComplexF1F1 alpha1 = alpha2 Basic 51-point"
}

basic_test_complex_equal! {
        name: basic_gauss_kronrod_smooth_function_equal_re_im_61_point,
        function: util::ComplexF1F1,
        alpha: f64 => 2.6,
        rule: gk61,
        lower: 0.0,
        upper: 1.0,
        exp_result: 7.716049382713800753E-02,
        exp_error: 1.566060362296155616E-12,
        abs_tolerance: 1.0e-15,
        rel_tolerance: 1.0e-5,
        "ComplexF1F1 alpha1 = alpha2 Basic 61-point"
}