rsl-interpolation 0.1.19

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

//! For 2d interpolation profiling.

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

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

const N: usize = 100_000_000;

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

pub fn main() -> Result<(), InterpolationError> {
    let acc = &mut Accelerator2d::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 = BicubicInterpolator::build(&xa, &ya, &za)?;

    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, acc)?);
        // acc.reset();
        black_box(interp2d.eval_deriv_x(&xa, &ya, &za, x, y, acc)?);
        // acc.reset();
        black_box(interp2d.eval_deriv_y(&xa, &ya, &za, x, y, acc)?);
        // acc.reset();
    }
    let end = start.elapsed();
    println!("{end:?}");
    Ok(())
}