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