Function newton_rootfinder::solver::compute_jacobian_from_finite_difference

pub fn compute_jacobian_from_finite_difference<M, D>(
    model: &mut M,
    perturbations: &OVector<f64, D>,
    update_residuals: &ResidualsConfig<'_>
) -> Result<OMatrix<f64, D, D>, ModelError<M, D>>where
    M: Model<D>,
    D: Dim,
    DefaultAllocator: Allocator<f64, D>,
    DefaultAllocator: Allocator<f64, D, D>,

Evaluate a jacobian per forward finite difference when perturbation step eps is provided

This function has been made public for testing purpose only