nathanrooy / bwo

Black Widow Optimization implemented in pure Python.

Black Widow Optimization

From the abstract:

...a novel meta-heuristic algorithm suitable for continuous nonlinear optimization problems. The proposed method, Black Widow Optimization Algorithm (BWO), is inspired by the unique mating behavior of black widow spiders. This method includes an exclusive stage, namely, cannibalism. Due to this stage, species with inappropriate fitness are omitted from the circle, thus leading to early convergence. BWO algorithm is evaluated on 51 various benchmark functions to verify its efficiency in obtaining the optimal solutions for the problems. The obtained results indicate that the proposed algorithm has numerous advantages in different aspects such as early convergence and achieving optimized fitness value compared to other algorithms.

Installation

pip install bwo

or

pip install git+https://github.com/nathanrooy/bwo

Usage

As a simple example, let's find the minimum of the single objective sphere function available in the Landscapes Python package.

pip install landscapes

Next, let's import everything we need.

>>> from bwo import minimize
>>> from landscapes.single_objective import sphere

Now, we just need to call the minimize function. For this particular example, let's optimize across 5 degrees of freedom.

>>> fbest, xbest = minimize(sphere, dof=5)

Where fbest is the best function value achieved during optimization, and xbest is the function input which yielded fbest.

You can also minimize a function given boundary constraints, represented by a list of tuples. Each tuple must follow the (min, max) format.

>>> bounds = [(-1,1),(-2,2), (-3,3), (-4,4), (-5,5)]
>>> fbest, xbest = minimize(sphere, bounds=bounds)

To display the convergence, simply set the disp flag to true as such:

>>> fbest, xbest = minimize(sphere, dof=5, maxiter=20, disp=True)

Which should yield something like the following:

> ITER:  1 | GBEST: 0.561091
> ITER:  2 | GBEST: 0.358288
> ITER:  3 | GBEST: 0.103173
> ITER:  4 | GBEST: 0.008892
> ITER:  5 | GBEST: 0.008892
> ITER:  6 | GBEST: 0.008887
> ITER:  7 | GBEST: 0.005540
> ITER:  8 | GBEST: 0.005540
> ITER:  9 | GBEST: 0.001907
> ITER: 10 | GBEST: 0.001221
> ITER: 11 | GBEST: 0.001147
> ITER: 12 | GBEST: 0.001056
> ITER: 13 | GBEST: 0.000874
> ITER: 14 | GBEST: 0.000874
> ITER: 15 | GBEST: 0.000760
> ITER: 16 | GBEST: 0.000760
> ITER: 17 | GBEST: 0.000731
> ITER: 18 | GBEST: 0.000731
> ITER: 19 | GBEST: 0.000719
> ITER: 20 | GBEST: 0.000719

Options

minimize(func, x0=None, dof=None, bounds=None, pp=0.6, cr=0.44, pm=0.4, npop=10, disp=False, maxiter=50)

• func (callable) : The objective function to be minimized.
• x0 (list) : Initial guess (optional).
• dof (int) : degrees of freedom (optional)
• bounds(list of tuples) : boundary constraints as specified as a list of (xi_min, xi_max) tuples.
• pp (float) : procreating percentage [0, 1].
• cr (float) : cannibalism rate [0, 1].
• pm (float) : mutation rate [0, 1].
• maxiter (int) : maximum number of iterations.
• disp (bool) : output intermediate results at each iteration.

References

@article{article,
author = {Hayyolalam, Vahideh and Pourhaji Kazem, Ali Asghar},
year = {2019},
month = {10},
pages = {103249},
title = {Black Widow Optimization Algorithm: A novel meta-heuristic approach for solving engineering optimization problems ✩},
volume = {87},
journal = {Engineering Applications of Artificial Intelligence},
doi = {10.1016/j.engappai.2019.103249}
}