fn phase_linearity_edge(
phase: &[f64],
i: usize,
j: usize,
k: usize,
ni: usize,
nj: usize,
nk: usize,
_dim: usize,
nx: usize,
ny: usize,
nz: usize,
) -> f64Expand description
Phase linearity for an edge between two voxels.
Computes the product of two triplet linearities: one looking “behind” the edge
and one looking “ahead”. Matches ROMEO.jl phaselinearity(P, i, j).