poly-cool 0.3.4

Polynomial root-finding
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

use criterion::{Criterion, criterion_group, criterion_main};
use polycool::Poly;
use std::hint::black_box;

pub fn eval(c: &mut Criterion) {
    let mut group = c.benchmark_group("eval");

    let coeffs_2 = [1.0, 2.0, 3.0];
    let coeffs_3 = [1.0, 2.0, 3.0, 4.0];
    let coeffs_4 = [1.0, 2.0, 3.0, 4.0, 5.0];

    group.bench_with_input("eval 2", &2, |b, _| {
        let p2 = Poly::new(coeffs_2);
        b.iter(|| black_box(p2).eval(black_box(1.0)))
    });
    group.bench_with_input("eval 3", &2, |b, _| {
        let p3 = Poly::new(coeffs_3);
        b.iter(|| black_box(p3).eval(black_box(1.0)))
    });
    group.bench_with_input("eval 3 alt", &2, |b, _| {
        let p3 = Poly::new(coeffs_3);
        b.iter(|| black_box(p3).eval_opt(black_box(1.0)))
    });
    group.bench_with_input("eval 2 and 3", &2, |b, _| {
        let p3 = Poly::new(coeffs_3);
        let p2 = Poly::new(coeffs_2);
        b.iter(|| {
            (
                black_box(p3).eval_opt(black_box(1.0)),
                p2.eval(black_box(1.0)),
            )
        })
    });
    group.bench_with_input("eval 2 and 3 alt", &2, |b, _| {
        let p3 = Poly::new(coeffs_3);
        let p2 = Poly::new(coeffs_2);
        b.iter(|| black_box(p3).eval_with_deriv_opt(&p2, black_box(1.0)))
    });
    group.bench_with_input("eval 4", &2, |b, _| {
        let p4 = Poly::new(coeffs_4);
        b.iter(|| black_box(p4).eval(black_box(1.0)))
    });
    // for coeffs in [coeffs_2, coeffs_3, coeffs_4] {
    //     let poly = Poly::new(coeffs
    //     group.bench_with_input(
    //         format!("roots_between {accuracy:?}"),
    //         &accuracy,
    //         |b, accuracy| b.iter(|| black_box(poly).roots_between(-1.0, 4.0, *accuracy)),
    //     );
    // }
}

criterion_group!(benches, eval);
criterion_main!(benches);