pub fn forward_hessian_nograd_sparse_const<const N: usize, F>( x: &[F; N], f: &'_ dyn Fn(&[F; N]) -> Result<F, Error>, indices: Vec<[usize; 2]>, ) -> Result<[[F; N]; N], Error>where F: Float + FromPrimitive + AddAssign,