pub fn compute_d(x_0: &Array1<f64>, grad: &Array1<f64>) -> Array1<f64>
Auxiliary function to calculate Delta, as defined in Boros2007