Function forward_jacobian_pert

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