quad-rs 0.5.0

Adaptive Gauss-Kronrod Integration in Rust
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

use num_complex::Complex;
use quad_rs::{Contour, Integrable, IntegratorConfig, integrate_complex};

struct Identity;

impl Integrable for Identity {
    type Float = f64;
    type Input = Complex<f64>;
    type Output = Complex<f64>;

    fn integrand(&self, z: &Complex<f64>) -> Complex<f64> {
        *z
    }
}

fn main() -> Result<(), Box<dyn std::error::Error>> {
    let contour = Contour::piecewise_linear(vec![Complex::new(0.0, 0.0), Complex::new(1.0, 1.0)]);

    let result = integrate_complex(Identity, contour, IntegratorConfig::default())?;

    println!("Integral estimate: {}", result.integral);
    println!("Expected:          {}", Complex::new(0.0, 1.0));
    println!("Error estimate:    {}", result.error);
    println!("Evaluations:       {}", result.evaluations);
    println!("Refinements:       {}", result.refinements);

    Ok(())
}