GlobalOptimization[GlobalSolve]  find a global solution to a nonlinear program

Calling Sequence


GlobalSolve(obj, constr, bd, opts)
GlobalSolve(obj, bd, opts)
GlobalSolve(opfobj, ineqcon, eqcon, opfbd, opts)
GlobalSolve(opfobj, opfbd, opts)


Parameters


obj



algebraic; objective function

constr



set(relation) or list(relation); constraints

bd



sequence of name = range; bounds for all variables

opfobj



procedure; objective function

ineqcon



set(procedure) or list(procedure); inequality constraints

eqcon



set(procedure) or list(procedure); equality constraints

opfbd



sequence of ranges; bounds for all variables

opts



(optional) equation(s) of the form option = value where option is one of avgstopstepwidth, evaluationlimit, iterationlimit, initialpoint, invcrossprobability, maximize, method, numexperiments, numsigma, nugget, objectivetarget, optimalitytolerance, optsearch, populationsize, randomseed, thetamethod, theta, targetweight, or timelimit; specify options for the GlobalSolve command





Description


•

The GlobalSolve command computes a global solution to a nonlinear program (NLP) over a bounded region. An NLP involves the minimization (or maximization) of an objective function, possibly subject to inequality and equality constraints. For a more detailed explanation of the solution obtained, see the following Notes section.

•

The global solver minimizes a merit function, incorporating a penalty term for the constraints, and provides both differential evolution and adaptive stochastic search methods. The global search phase is followed by a local search phase to refine the solution. No derivatives are required.


The solver is designed to search the specified region for a global solution to a nonconvex optimization problem. If the optimization problem is
convex
(for example, a linear program) or a local solution is acceptable, it is recommended that you use the commands for local optimization in the Optimization package. The Optimization package commands, which are more efficient, can compute global solutions to convex problems.

•

The first and second calling sequences use the algebraic form of input. The first parameter obj is the objective function, which must be an algebraic expression.


Constraints can be provided using the constr parameter. This is a set or list of relations (of type `<=` or `=`) involving the problem variables. The problem variables are the indeterminates of type name found in obj and constr. They can also be specified using the variables option. For unconstrained problems, omit the constr parameter.


Bounds on the variables must be given as additional arguments, each of the form varname= varrange where varname is a variable name and varrange is its range. There must be exactly one argument for each problem variable, and the endpoints of each range must evaluate to finite numeric values.

•

The third and fourth calling sequences use the operator form of input. The objective function opfobj must be a procedure that accepts floatingpoint parameters representing the problem variables x1, x2, ..., xn and returns a float.


Inequality and equality constraints are provided using the ineqcon and eqcon parameters. An inequality constraint is specified by a procedure v in ineqcon that has the same form as opfobj and returns the lefthandside value of the constraint. Similarly, an equality constraint is specified by a procedure w in eqcon. Either ineqcon or eqcon can be an empty list or set. For unconstrained problems, omit both of these parameters.


Bounds on the variables must be provided. These must be a sequence of exactly ranges corresponding in order to x1, x2, ..., xn. The endpoints of each range must evaluate to finite numeric values.

•

Maple returns the solution as a list containing the final minimum (or maximum) value and a point (the extremum). If the input is in algebraic form, the point is a list containing elements of the form where varname is a problem variable and value is its final value. If the input is in operator form, the point is a Vector containing the values of the problem variables.



Options



The opts argument can contain one or more of the following options not specific to a particular solving method. The options specific to a particular solving method are described in more detail in the GlobalOptimization/Options help page.

•

evaluationlimit = posint  Set the maximum number of iterations performed by the Differential Evolution method. The global search phase terminates if this limit is reached.

•

feasibilitytolerance = positive and numeric  Set the allowed violation of constraints.

•

initialpoint = set(equation), list(equation), or list(numeric)  Use the provided initial point, which is a set or list of equations (for algebraic form input) or a list of exactly values (for operator form input). This is supported by the Differential Evolution method.

•

maximize or maximize = truefalse  Maximize the objective function when m is 'true' and minimize when m is 'false'. The option 'maximize' is equivalent to 'maximize'='true'. The default is 'maximize'='false'.

•

method = diffevol, ego  Set the global search algorithm: Differential Evolution (method = diffevol), or Efficient Global Optimization (method = ego). The default is method = diffevol.

•

objectivetarget = numeric  Set an acceptable target value for the objective function. If the objective function achieves this value, the search terminates.

•

timelimit = posint  Set the maximum computation time, in seconds, for the global solver.

•

variables = list(name) or set(name)  Specify the problem variables when the objective function is in algebraic form.



Notes


•

For more information on the methods used by the global solver, with suggestions for achieving best performance, see the GlobalOptimization/Computation help page.

•

The global solver searches for the optimal solution until one of the termination criteria is met. Then either the best available solution is returned or an error message is displayed stating that a solution could not be obtained. The termination criteria can be set using options. Otherwise, default values for these options are applied. In particular, the evaluationlimit option must be set to a sufficiently high value for difficult optimization problems or unexpected answers may be produced.

•

In the case that an error occurs because the time limit has been exceeded, the last solution computed may be retrieved using the GetLastSolution command.

•

The computation is performed in floatingpoint. Therefore, all data provided must have type realcons and all returned solutions are floatingpoint, even if the problem is specified with exact values. The solver uses externally called code that works with hardware floats, but it is possible to evaluate the objective function and the constraints in Maple with higher precision. For details, see the GlobalOptimization/Computation help page.



Examples


>


Find the global solution to an unconstrained nonlinear minimization problem.
>


 (1) 
Find the global solution to a constrained twovariable minimization problem.
>


 (2) 
Find the global solution to a constrained leastsquares minimization problem.
>


 (3) 
Find the global minimum using the operator input form.
>


 (4) 
Initial points can be provided to find the minimum.
>


 (5) 
Use the maximize option to maximize the objective function.
>


 (6) 
GlobalSolve solves realvalued optimization problems. Trying to solve a complex problem can produce an error message.
>


However, you can apply the Re command to use the real portion of the complexvalued expressions in a complex optimization problem to define a different problem.
>


 (7) 
Using Re can create distinct problems with distinct solutions. Notice the difference between the previous example and the following.
>


 (8) 
Set infolevel to or higher to display details about the solution procedure.
>


>


GlobalSolve: calling NLP solver
GlobalSolve: calling global optimization solver
GlobalSolve: number of problem variables 1
GlobalSolve: number of nonlinear inequality constraints 0
GlobalSolve: number of nonlinear equality constraints 0
GlobalSolve: OptimusDEVOL OptimusDEVOL
GlobalSolve: maximum iterations 50
GlobalSolve: population size 20
GlobalSolve: NULL 100
GlobalSolve: trying evalf mode
GlobalSolve: performing local refinement
 
 (9) 

