If not set, compute epsilon based on half_range

compute_v_rel(v_rel, half_range, n)