let pi_other = reach[1 - p.into_idx()];

            let pi_other = reach.iter().product::<f32>() / pi_self;