pub fn min_componentwise_quadratic_solve<const R: usize, const C: usize>( a: &SMatrix<Real, R, C>, b: &SMatrix<Real, R, C>, c: Real, sol_half_range: (Real, Real), ) -> Real