On Fri, May 14, 2010 at 09:45:55AM +0100, Gary wrote: > Hi Olaf, > > Thanks for responding so quickly. I'm sorry if I didn't make it clear > in my previous message, but the script calls an external program and > has a lot of code that handles that aspect of the procedure. In this > script, the bounds are only sometimes violated. Unfortunately, I > cannot supply you with a simple self-contained script that reproduces > the error. I could post the entire script, but due to licensing issues > I cannot post the external program.
I suggest you insert into the script, just after the call to leasqr, a check for what you call an error (violation of bounds?), interrupting the script in case of that error: if (check_for_the_error_here) keyboard; endif >From the debugging prompt you'll get, check the values of all arguments to leasqr and of all return-values of leasqr. Post all this information to the list, together with the exact command line which calls leasqr. Do you mean that the model function passed to leasqr calls an external program which is restrictively licensed? This might make it impossible to help. If possible, post the model function and information on any global variables set for the model function. > I have started rectifying the parameters before passing them to the > external program. This allows the previously failing optimisation run > to complete and seems to produce sensible results. Though I cannot say > the optimisation is working correctly, A failing optimization does not necessarily mean that the optimizer is buggy. I'm more interested in checking that bounds work correctly. Olaf ------------------------------------------------------------------------------ _______________________________________________ Octave-dev mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/octave-dev
