I'm using NLEQ to estimate a distortion model for a camera. I already have a Scilab code that estimates the distortion so I know a good solution and alglib's nonlinear solver NLEQ also finds the known solution but returns report code, -4 ERROR: algorithm has converged to the stationary point Xf which is local minimum of f=F[0]^2+...+F[m-1]^2, but is not solution of nonlinear system.
In this problem, the number of equations m=224 and the dimension of the state is n=61. How does NLEQ know that it has only found a local minimum and that there is a global minimum?
|