let light_cones =
    let total_depth = depth + perturbation_range +1 in
    Array.init length
    (fun i -> if discard_cell i then [||]
     else Measure.Light_cones.build_lc_index neighbourhoods total_depth i)