soilrust 0.7.10

A geotechnical engineering library for soil mechanics calculations 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

/// Performs linear interpolation for a given x value based on provided x and y vectors.
///
/// # Arguments
/// * `x_values` - Array of x-axis values (must be sorted)
/// * `y_values` - Array of y-axis values
/// * `x` - The x value for which to interpolate
///
/// # Returns
/// * Interpolated y value as f64
///
/// # Panics
/// If x_values and y_values lengths are not equal or x is out of range.
pub fn interp1d(x_values: &[f64], y_values: &[f64], x: f64) -> f64 {
    assert_eq!(
        x_values.len(),
        y_values.len(),
        "x_values and y_values must have the same length"
    );

    if x <= x_values[0] {
        return y_values[0];
    }
    if x >= x_values[x_values.len() - 1] {
        return y_values[y_values.len() - 1];
    }

    for i in 0..x_values.len() - 1 {
        let x0 = x_values[i];
        let x1 = x_values[i + 1];
        let y0 = y_values[i];
        let y1 = y_values[i + 1];

        if x >= x0 && x <= x1 {
            return y0 + (y1 - y0) * (x - x0) / (x1 - x0);
        }
    }

    panic!("Interpolation error: x-value out of interpolation range");
}