provide an option to handle the zero Jacobian case

[tpapp]

It can happen that we are at a (local) minimum wrt to the sum of squared residuals, but don't want to stop the optimizer if those are not close to zero. Currently this errors in finding the Cauchy point.