Investigate failures of block decomposition when computing steady states

The file solve_v8.mod fails computing the terminal steady state with trust region, while it works with fsolve. Part of the problem is that dogleg returns NaN due to 0*Inf. parameter.mat steadystate.mat

This seems to be caused by blocks computed from the Dulmage-Mendelsohn decomposition having singular Jacobians.