plot3d 0.1.10

Utilities for reading, writing, and manipulating NASA PLOT3D structured grids.
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

//! Point matching utilities for structured grid faces.

use crate::Float;

/// Find the indices of the closest point in a 2D face grid to `(x, y, z)`.
///
/// The face grid is represented by flat arrays `x2`, `y2`, `z2` of length
/// `nu * nv`, stored row-major (u varies fastest).
///
/// Returns `Some((u, v))` if the minimum distance is within `tol`,
/// otherwise `None`.
#[allow(clippy::too_many_arguments)]
pub fn point_match(
    x: Float,
    y: Float,
    z: Float,
    x2: &[Float],
    y2: &[Float],
    z2: &[Float],
    nu: usize,
    _nv: usize,
    tol: Float,
) -> Option<(usize, usize)> {
    debug_assert_eq!(x2.len(), y2.len());
    debug_assert_eq!(x2.len(), z2.len());

    let mut best_dist = Float::INFINITY;
    let mut best_idx = 0usize;

    for (idx, ((&px, &py), &pz)) in x2.iter().zip(y2.iter()).zip(z2.iter()).enumerate() {
        let dx = x - px;
        let dy = y - py;
        let dz = z - pz;
        let d = (dx * dx + dy * dy + dz * dz).sqrt();
        if d < best_dist {
            best_dist = d;
            best_idx = idx;
        }
    }

    if best_dist < tol {
        let u = best_idx % nu;
        let v = best_idx / nu;
        Some((u, v))
    } else {
        None
    }
}