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