pub fn central_jacobian_vec_prod_vec<F>( x: &[F], fs: &'_ dyn Fn(&Vec<F>) -> Result<Vec<F>, Error>, p: &[F], ) -> Result<Vec<F>, Error>where F: Float + FromPrimitive,