pub fn forward_hessian_const<const N: usize, F>( x: &[F; N], grad: &'_ dyn Fn(&[F; N]) -> Result<[F; N], Error>, ) -> Result<[[F; N]; N], Error>where F: Float + FromPrimitive,