ruvector-math 2.0.6

Advanced mathematics for next-gen vector search: Optimal Transport, Information Geometry, Product Manifolds
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

//! Benchmarks for spectral methods

use criterion::{black_box, criterion_group, criterion_main, BenchmarkId, Criterion};
use ruvector_math::spectral::{ChebyshevExpansion, ChebyshevPolynomial};

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

    for degree in [10, 20, 50, 100] {
        group.bench_with_input(
            BenchmarkId::new("degree", degree),
            &degree,
            |bench, &deg| {
                let poly = ChebyshevPolynomial::new(deg);
                bench.iter(|| poly.eval(black_box(0.5)));
            },
        );
    }

    group.finish();
}

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

    for degree in [10, 20, 50, 100] {
        group.bench_with_input(
            BenchmarkId::new("degree", degree),
            &degree,
            |bench, &deg| {
                bench.iter(|| ChebyshevPolynomial::eval_all(black_box(0.5), black_box(deg)));
            },
        );
    }

    group.finish();
}

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

    for degree in [10, 20, 50] {
        group.bench_with_input(
            BenchmarkId::new("from_function", degree),
            &degree,
            |bench, &deg| {
                bench.iter(|| ChebyshevExpansion::from_function(|x| x.sin(), black_box(deg)));
            },
        );
    }

    group.finish();
}

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

    for degree in [10, 20, 50] {
        for t in [0.1, 1.0, 10.0] {
            group.bench_with_input(
                BenchmarkId::new(format!("deg={}_t={}", degree, t), degree),
                &(degree, t),
                |bench, &(deg, t)| {
                    bench.iter(|| ChebyshevExpansion::heat_kernel(black_box(t), black_box(deg)));
                },
            );
        }
    }

    group.finish();
}

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

    for degree in [10, 20, 50, 100] {
        let expansion = ChebyshevExpansion::from_function(|x| x.sin(), degree);

        group.bench_with_input(
            BenchmarkId::new("degree", degree),
            &expansion,
            |bench, exp| {
                bench.iter(|| exp.eval(black_box(0.5)));
            },
        );
    }

    group.finish();
}

criterion_group!(
    benches,
    bench_chebyshev_eval,
    bench_chebyshev_eval_all,
    bench_chebyshev_expansion,
    bench_heat_kernel,
    bench_clenshaw_eval,
);
criterion_main!(benches);