fn k2_kernel(dims: (usize, usize, usize), vs: (f64, f64, f64)) -> Vec<f64>
STI’s spectral Laplacian kernel k² = kx² + ky² + kz², k in cycles/mm, in FFT order.
k² = kx² + ky² + kz²
k