pub fn barrier_method( f: fn(&[f64]) -> f64, inequalities: &[fn(&[f64]) -> f64], x: &[f64], t: f64, ) -> f64