Function block_update 

pub fn block_update<F>(
    jac: &mut Array2<F>,
    delta_y: &Array1<F>,
    delta_f: &Array1<F>,
    block_size: usize,
)
where
    F: IntegrateFloat,

Performs a block-update of the Jacobian using block structure