pub fn find_face_edges(
x: &[Float],
y: &[Float],
z: &[Float],
pmax: usize,
qmax: usize,
) -> Vec<FaceDiff>Expand description
Compute forward and backward differences along each direction for a 2D face.
x, y, z are flat arrays of length pmax * qmax, stored row-major
(p varies fastest within each row of q).