fn irls_weights( ws: &mut Ws, x: &[f32], mx: &[f32], my: &[f32], mz: &[f32], vr: &mut [f32], )
IRLS weights Vr = 1/sqrt(|M·∇χ|² + eps) for one TV term.