fn drop_large_rim_values(lap: &mut [f64], rim: &[u8])
Set the 20 % of rim Laplacian values largest in magnitude to zero: keep |L| ≤ the sorted rim magnitude at 1-based index round(0.8 N).
|L| ≤
round(0.8 N)