CGO arbfMIP

This page is part of the CGO Manual. See CGO Manual.

Purpose

Solve general constrained mixed-integer global black-box optimization problems with costly objective functions.

The optimization problem is of the following form

${\begin{array}{cccccc}\min \limits _{x}&f(x)&&&&\\s/t&x_{L}&\leq &x&\leq &x_{U}\\&b_{L}&\leq &Ax&\leq &b_{U}\\&c_{L}&\leq &c(x)&\leq &c_{U}\\&~x_{j}\in \mathbb {N} \ ~~\forall j\in \mathbb {I} ,\\\end{array}}$

where $f(x)\in \mathbb {R}$; $x_{L},~x,~x_{U}\in \mathbb {R} ^{d}$; the $m_{1}$ linear constraints are defined by $A\in \mathbb {R} ^{m_{1}\times d}$, $b_{L},~b_{U}\in \mathbb {R} ^{m_{1}}$; and the $m_{2}$ nonlinear constraints are defined by $c_{L},~c(x),~c_{U}\in ~\mathbb {R} ^{m_{2}}$. The variables $x_{I}$ are restricted to be integers, where $\mathbb {I}$ is an index subset of $\{1,\ldots ,d\},$ possibly empty. It is assumed that the function $f(x)$ is continuous with respect to all variables, even if there is a demand that some variables only take integer values. Otherwise it would not make sense to do the surrogate modeling of $f(x)$ used by all CGO solvers.

f (x) is assumed to be a costly function while c(x) is assumed to be cheaply computed. Any costly constraints can be treated by adding penalty terms to the objective function in the following way:

${\begin{array}{ll}\min \limits _{x}&p(x)=f(x)+\sum \limits _{j}w_{j}\max \left(0,~c^{j}(x)-c_{U}^{j},~~c_{L}^{j}-c^{j}(x)\right),\end{array}}$

where weighting parameters w_j have been added. The user then returns p(x) instead of f (x) to the CGO solver.

Calling Syntax

Result = arbfMIP(Prob,varargin) 
Result = tomRun('arbfMIP', Prob);

Description of Inputs

Problem structure

The following fields are used in the problem description structure Prob:

Input	Description
Name	See Common input for all CGO solvers
FUNCS.f
FUNCS.c
x_L
x_U
b_L
b_U
A
c_L
c_U
WarmStart
MaxCPU
user
PriLevOpt
f_Low
optParam
CGO	See the table below but also this table for input common to all CGO solvers
GO	See common input for all CGO solvers
MIP	See common input for all CGO solvers
varargin	Additional parameters to arbfmip are sent to the costly f(x)

- Special ARBF algorithm parameters in Prob.CGO -
rbfType	Selects type of radial basis function Value Type 1 Thin Plate Spline 2 Cubic Spline (default) 3 Multiquadric 4 Inverse multiquadric 5 Gaussian 6 Linear.
infStep	If =1, add search step with target value -inffirst in cycle. Default 0
TargetMin	Which minimum of several to pick in target value problem: Value Minimum picked 0 Use global minimum. 1 Use best interior local minima, if none use global minimum. 2 Use best interior local minima, if none use RBF interior minimum. 3 Use best minimum with lowest number of coefficients on bounds. Default is TargetMin = 3.
fStarRule	Global-Local search strategy. N = cycle length. Define min_sn as the global minimum on surface. Value fStar target value 1 min_sn - ((N - (n - nInit))/N )2 * Deltan (Default) 2 min_sn - (N - (n - nInit))/N * Deltan. Strategy 1 and 2 depends on Deltan estimate (see DeltaRule). 3 -inf-step, min_sn-k *0.1*|min_sn| k = N,...,0. If infStep true, addition of -inf-step first in cycle. Strategy names in Gutmanns thesis: III, II, I
DeltaRule	1 = Skip large f(x) when computing f(x) interval Delta. 0 = Use all points. If objType > 0, default DeltaRule = 0, otherwise default is 1.
eps_sn	Relative tolerance used to test if the minimum of surface, min_sn, is sufficiently lower than the best point (fMin) found. Default is eps_sn = 10-7.

Description of Outputs

Result structure

The output structure Result contains results from the optimization.
The following fields are set:

Field	Description
x_k	See Common output for all CGO solvers for details.
f_k
Iter
FuncEv
ExitText
ExitFlag	Always 0
Inform	Information parameter. Value Signification 0 Normal termination. 1 Function value f(x) is less than fGoal. 2 Error in function value f (x), |f - fGoal| <= fTol, fGoal = 0. 3 Relative Error in function value f (x) is less than fTol, i.e. |f - fGoal|/|fGoal| <= fTol. 6 All sample points same as previous point for the last 11 iterations. 7 All feasible integers tried. 9 Max CPU Time reached.
CGO	Subfield WarmStartInfo saves warm start information, the same information as in cgoSave.mat, see Common output for all CGO solvers#WSInfo.

Description

arbfMIP implements the Adaptive Radial Basis Function (ARBF) algorithm. The ARBF method handles linear equality and inequality constraints, and nonlinear equality and inequality constraints, as well as mixed-integer problems.

M-files Used

daceInit.m, iniSolve.m, endSolve.m, conAssign.m, glcAssign.m, snSolve.m, gnSolve.m, expDesign.m.

MEX-files Used

tomsol

See Also

rbfSolve.m and ego.m

Warnings

Observe that when cancelling with CTRL+C during a run, some memory allocated by arbfMIP will not be deallocated. To deallocate, do:

''>> ''clear cgolib