pub fn forward_jacobian_vec_prod<const N: usize, const M: usize, F>( f: &'_ dyn Fn(&[F; N]) -> Result<[F; M], Error>, ) -> impl Fn(&[F; N], &[F; N]) -> Result<[F; M], Error> + '_where F: Float + FromPrimitive,