Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Failing tests related to fminsearch #7

Open
jbect opened this issue Aug 23, 2020 · 1 comment
Open

Failing tests related to fminsearch #7

jbect opened this issue Aug 23, 2020 · 1 comment

Comments

@jbect
Copy link
Member

jbect commented Aug 23, 2020

There currently two failing tests under Matlab (tested with R2016b) when fminsearch is used as last recourse for box-constrained problems.

The first one is in stk_optim_testmin_box.m and simply indicates that the solution not not satisfy the box constraint:

>> stk_test stk_optim_testmin_box
***** assert (stk_optim_testmin_box (algo));
!!!!! test failed
Error using stk_test>eval_test_code (line 505)
Assertion failed.

The second one is in stk_param_estim, and is probably a symptom of a more severe problem:

>> stk_test stk_param_estim
***** test  % noiseless
 zi = stk_feval (f, xi);
 param1 = stk_param_estim (model, xi, zi, param0);
 param2 = stk_param_estim (model, xi, zi, param0, [], @stk_param_relik);
 % We cannot assume a DETERMINISTIC optimization algorithm
 % (for some reason, Octave's sqp is not exactly deterministic)
 assert (stk_isequal_tolrel (param1, param2, 1e-2))
!!!!! test failed
Error using stk_cholcov (line 66)
A contains NaNs or Infs.

Perhaps should we try to provide a better default optimizer for Matlab users that have neither the Optimization toolbox nor MOSEK?

(ref original ticket on Sourceforge: https://sourceforge.net/p/kriging/tickets/65/)

@jbect
Copy link
Member Author

jbect commented Aug 23, 2020

We already have a copy of Octave's sqp inside STK. Therefore, a solution would be to make it possible to use a different QP solvers when neither MOSEK nor the Optimization toolbox are available. (For instance provide a copy of Octave's qp solver as well & recommend to install GLPKmex?)

Another option woulld be to rely on NLopt++.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant