rsl-interpolation 0.1.18

A re-write of GSL's Interpolation Routines 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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48

//! For 2d interpolation profiling

use std::hint::black_box;
use std::time::Instant;

use ndarray::Array1;
use rsl_interpolation::*;

const N: usize = 10_000_000;
const TYP: Bilinear = Bilinear; // keep to "Bilinear", otherwise needs BLAS to build and publish

// ===============================================================================================

pub fn main() {
    let mut xacc = Accelerator::new();
    let mut yacc = Accelerator::new();
    let mut cache = Cache::new();
    let xa: Vec<f64> = (0..100).map(f64::from).collect();
    let ya: Vec<f64> = (0..200).map(f64::from).collect();
    let za: Vec<f64> = (0..(100 * 200)).map(f64::from).collect();
    let interp2d = TYP.build(&xa, &ya, &za).unwrap();

    let xs = Array1::linspace(1.0, 99.0, N);
    let ys = Array1::linspace(1.0, 199.0, N);

    let start = Instant::now();
    for i in 0..N {
        let x = black_box(xs[i]);
        let y = black_box(ys[i]);
        black_box(
            interp2d
                .eval(&xa, &ya, &za, x, y, &mut xacc, &mut yacc, &mut cache)
                .unwrap(),
        );
        black_box(
            interp2d
                .eval_deriv_x(&xa, &ya, &za, x, y, &mut xacc, &mut yacc, &mut cache)
                .unwrap(),
        );
        black_box(
            interp2d
                .eval_deriv_y(&xa, &ya, &za, x, y, &mut xacc, &mut yacc, &mut cache)
                .unwrap(),
        );
    }
    let end = start.elapsed();
    println!("{end:?}")
}