fn compute_relative_change(chi: &[f32], chi_prev: &[f32], mask: &[u8]) -> f32
Compute relative change for convergence check