4.10. pyopus.optimizer.de — Box constrained differential evolution optimizer

Box constrained differential evolution global optimizer (PyOPUS subsystem name: DEOPT)

Published first in [de].

[de]	Storn R., Price K.: Differential evolution - a simple and efficient heuristic for global optimization over continuous spaces. Journal of Global Optimization vol. 11, pp. 341–359, 1997.

class pyopus.optimizer.de.DifferentialEvolution(function, xlo, xhi, debug=0, fstop=None, maxiter=None, maxSlaves=None, minSlaves=0, maxGen=None, spawnerLevel=1, populationSize=100, w=0.5, pc=0.3)

Box constrained differential evolution global optimizer class

If debug is above 0, debugging messages are printed.

The lower and upper bound (xlo and xhi) must both be finite.

maxGen is the maximal number of generations. It sets an upper bound on the number of iterations to be the greater of the following two values: maxiter and .

populationSize is the number of inidividuals (points) in the population.

w is the differential weight between 0 and 2.

pc is the crossover probability between 0 and 1.

If spawnerLevel is not greater than 1, evaluations are distributed across available computing nodes (that is unless task distribution takes place at a higher level).

Infinite bounds are allowed. This means that the algorithm behaves as an unconstrained algorithm if lower bounds are and upper bounds are . In this case the initial population must be defined with the reset() method.

See the BoxConstrainedOptimizer for more information.

check()

Checks the optimization algorithm’s settings and raises an exception if something is wrong.

generatePoint(i)

Generates a new point through mutation and crossover.

initialPopulation()

Constructs and returns the initial population.

Fails if any bound is infinite.

reset(x0=None)

Puts the optimizer in its initial state.

x is eiter the initial point (which must be within bounds and is ignored) or the initial population in the form of a 2-dimensional array or list where the first index is the population member index while the second index is the component index. The initial population must lie within bounds xlo and xhi and have populationSize members.

If x is None the initial population is generated automatically. In this case xlo and xhi determine the dimension of the problem. See the initialPopulation() method.

run()

Run the algorithm.