Re: [MATH] Comment on MATH-352

2010-03-10 Thread Dimitri Pourbaix
Luc, So I would suggest to put a customizable threshold with a default value of MathUtils.SAFE_MIN and a setQRRankingThreshold(double threshold) method. This will preserve compatibility. With this setting, one could check all existing tests should pass. The solution proposed by the reporter does

Re: [MATH] Comment on MATH-352

2010-03-10 Thread Luc Maisonobe
Dimitri Pourbaix a écrit : > Hi, > > I looked at the suggestion related to MATH-352 regarding the robustness of > the determination of the rank of the Jacobian in Levenberg-Marquardt. It > does make sense indeed to replace ak2 == 0 with something more robust like > ak2<= treshold. However, if on

[MATH] Comment on MATH-352

2010-03-10 Thread Dimitri Pourbaix
Hi, I looked at the suggestion related to MATH-352 regarding the robustness of the determination of the rank of the Jacobian in Levenberg-Marquardt. It does make sense indeed to replace ak2 == 0 with something more robust like ak2<= treshold. However, if one uses the epsilon machine for the thr