interpn 0.11.1

N-dimensional interpolation/extrapolation methods, no-std and no-alloc compatible.
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

use num_traits::Float;

use crate::mul_add;

#[inline]
pub(crate) fn dot4<T: Float>(weights: [T; 4], vals: [T; 4]) -> T {
    let lo = mul_add(weights[1], vals[1], weights[0] * vals[0]);
    let hi = mul_add(weights[3], vals[3], weights[2] * vals[2]);
    lo + hi
}

#[inline]
pub(crate) fn hermite_basis<T: Float>(t: T) -> [T; 4] {
    let one = T::one();
    let two = one + one;
    let three = two + one;
    let t2 = t * t;
    let t3 = t2 * t;

    [
        mul_add(two, t3, mul_add(-three, t2, one)),
        mul_add(t, mul_add(t, t - two, one), T::zero()),
        mul_add(-two, t3, three * t2),
        mul_add(t2, t - one, T::zero()),
    ]
}