pub fn forward_jacobian_pert<F>( f: &'_ dyn Fn(&Vec<F>) -> Result<Vec<F>, Error>, ) -> impl Fn(&Vec<F>, &PerturbationVectors) -> Result<Vec<Vec<F>>, Error> + '_where F: Float + FromPrimitive + AddAssign,