Next Article in Journal
A Boundary Distance-Based Symbolic Aggregate Approximation Method for Time Series Data
Previous Article in Journal
Special Issue “Nonsmooth Optimization in Honor of the 60th Birthday of Adil M. Bagirov”: Foreword by Guest Editors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Evolution with Linear Bias Reduction in Parameter Adaptation

by
Vladimir Stanovov
*,†,
Shakhnaz Akhmedova
and
Eugene Semenkin
Department of System Analysis and Operation Research, Reshetnev Siberian State University of Science and Technology, Krasnoyarsk 660014, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Algorithms 2020, 13(11), 283; https://doi.org/10.3390/a13110283
Submission received: 19 October 2020 / Revised: 6 November 2020 / Accepted: 6 November 2020 / Published: 9 November 2020
(This article belongs to the Section Evolutionary Algorithms and Machine Learning)

Abstract

:
In this study, a new parameter control scheme is proposed for the differential evolution algorithm. The developed linear bias reduction scheme controls the Lehmer mean parameter value depending on the optimization stage, allowing the algorithm to improve the exploration properties at the beginning of the search and speed up the exploitation at the end of the search. As a basic algorithm, the L-SHADE approach is considered, as well as its modifications, namely the jSO and DISH algorithms. The experiments are performed on the CEC 2017 and 2020 bound-constrained benchmark problems, and the performed statistical comparison of the results demonstrates that the linear bias reduction allows significant improvement of the differential evolution performance for various types of optimization problems.

1. Introduction

The Computational Intelligence (CI) methods include a variety of approaches, such as Evolutionary Computation (EC), Fuzzy Logic (FL), and Neural Networks (NN). Despite the differences among these directions, one of the most important research areas that connects the parts if CI is the classical numerical optimization, and thus, the development of modern optimization techniques is one of the most important directions. The new heuristic optimization methods, including Evolutionary Algorithms (EA) and Swarm Intelligence (SI), are usually developed by introducing new algorithmic schemes [1] or parameter control or adaptation techniques [2]. Among the numerical optimization techniques, the Differential Evolution (DE) algorithm [3] has attracted the attention of a large number of researchers due to its simplicity and high performance [4]. The results of the annual competitions on numerical optimization, conducted within the IEEE Congress on Evolutionary Computation, show [5] that in the last few years, the winners have been mostly DE-based frameworks, which include hybridization or novel parameter tuning techniques.
Currently, one of the main directions of studies on differential evolution is the development of parameter adaptation techniques. The parameter tuning mechanism should rely on the search process quality, usually fitness values and their improvement, and for a tuning scheme to be efficient, it should be designed with respect to the algorithm properties, which leads to many different adaptation schemes developed for specific algorithms. The final goal of every tuning scheme is similar: provide the algorithm with the ability to tune parameters to suboptimal values at every stage of the search so that the best possible solution is reached at the end. Several recent surveys considered the existing variants of DE and their properties [6,7].
The parameter tuning techniques include parameter adaptation and parameter control methods: the former change parameter values based on feedback, while the latter follow a predetermined control strategy. Some of the recent advances in parameter adaptation include the success-history adaptation-based SHADE algorithm [8] and the distance-based DISH modification [9], while parameter control includes the linear population size reduction method in the L-SHADE approach [10], as well as rules for parameter control in the jSO algorithm [11]. All these modifications have shown that the potential of DE improvements is far from being exhausted, and thus, new investigations in this area are important.
In this study, the linear bias reduction technique is proposed, which changed the success-history-based parameter adaptation scheme by introducing the generalized Lehmer mean proposed in [12] and a rule for changing the coefficient in the Lehmer mean depending on the stage of the search. The experiments are performed on the benchmark functions of the CEC 2017 and CEC 2020 competitions, and the results are compared to the state-of-the-art approaches.
The paper is organized as follows: the second section contains the description of the modern DE variants and, then, proposes the linear bias reduction approach; the third section contains the experimental results and comparison; the fourth section contains the discussion; and the last section concludes the paper, providing the directions for further research.

2. Materials and Methods

2.1. Differential Evolution Algorithm

The original differential evolution algorithm was proposed by Storn and Price [13], and its main feature was the usage of differences between the vectors representing solutions, hence the name of the algorithm. It has received much attention from researchers mainly because it is easy to implement and it has high efficiency. The DE starts with a population of randomly initialized solutions x i , j = r a n d ( x m a x , j x m i n , j ) , where i = 1 , , N P , j = 1 , , D , x m i n , j and x m a x , j are the lower and upper bounds of the search area for variable j, NP is the initial population size, and D is the problem dimension. After initialization, the algorithm proceeds with the mutation, crossover, and selection operations. The theoretical considerations of all these steps were studied in [14].
There are several mutation schemes known for DE, and the basic strategy is called rand/1 and implemented as follows:
v i , j = x r 1 , j + F ( x r 2 , j x r 3 , j ) ,
where F is the scaling factor parameter and indexes r 1 , r 2 and r 3 are mutually different and randomly generated. There were several more advanced strategies proposed, including the current-to-pbest/1, introduced in JADE [15] and also used in the SHADE framework [8]:
v i , j = x i , j + F ( x p b e s t , j x i , j ) + F ( x r 1 , j x r 2 , j )
where p b e s t is the index of one of the p b 100 % best in terms of fitness individuals and p b e s t r 1 r 2 i .
The mutation strategies significantly influence the algorithm performance and require specific parameter tuning techniques; most of them randomly sample F and C r values. There were several modifications to the current-to-pbest/1 strategy, including the current-to-pbest-w/1 strategy introduced in jSO [11], which significantly improved the performance:
v i , j = x i , j + F w ( x p b , j x i , j ) + F ( x r 1 , j x r 2 , j ) ,
where F w is defined as follows:
F w = 0.7 F , if N F E < 0.2 N F E m a x , 0.8 F , if 0.2 N F E m a x N F E < 0.4 N F E m a x , 1.2 F , otherwise .
where N F E is the number of function evaluations and N F E m a x is the total computational resource. In L-SHADE-RSP (L-SHADE with Rank-based Selective Pressure) [16], the mutation strategy was modified with selective pressure resulting in current-to-pbest-w/r, which allowed faster convergence. In current-to-pbest-w/r, the last index r 2 is selected considering the fitness values of the individuals in the population, i.e., the rank values r a n k i = i , i = 1 N P are set, with the largest rank assigned to the best individual (smallest goal function value in the case of minimization). The probabilities of individuals being selected are calculated as follows:
p i s e l = r a n k i j = 1 N P r a n k j .
After the mutation step, the crossover is performed, where the newly generated donor vector v i is combined with the current vector x i . There are two types of crossover operations known for DE, namely the binomial and exponential crossover. The exponential crossover is mostly used when there are known connections between adjacent variables, while the binomial crossover is more popular and implemented as follows:
u i , j = v i , j , if r a n d ( 0 , 1 ) < C r or j = j r a n d x i , j , otherwise
The j r a n d index is randomly generated in [ 1 , D ] and is needed to make sure that at least one component is taken from the donor vector, and C r is the crossover rate parameter. The jSO algorithm additionally changes the C r generation, so that large values of the crossover rate are not allowed at the beginning of the search:
C r = m a x ( C r , 0.7 ) , if N F E < 0.25 N F E m a x , m a x ( C r , 0.6 ) , if 0.25 N F E m a x N F E and N F E < 0.5 N F E m a x .
After generating the trial solution u i , the bound constraint handling method is applied to make sure that all vectors are within the search bounds:
u i , j = x m i n j + x i , j 2 , if u i , j < x m i n j x m a x j + x i , j 2 , if u i , j > x m a x j .
After calculating the fitness of the trial solution f ( u i ) , the selection step is performed:
x i , j = u i , j , if f ( u i , j ) f ( x i , j ) x i , j , if f ( u i , j ) > f ( x i , j ) .
Most DE implementations use greedy selection, where the newly generated solution is accepted only if it improves its corresponding current solution x i .
In addition, JADE, SHADE, L-SHADE, jSO, L-SHADE-RSP, and other approaches use an external archive of solutions. The archive is initially empty and is filled with parent vectors, replaced during selection. The archive usually has the same size as the population, and after its size hits the maximum value, newly inserted solutions replace randomly chosen ones. The archived solutions are used in current-to-pbest/1 mutation strategies in the last random index: r 2 is either selected from the population or the archive.

2.2. Related Studies on Differential Evolution

The original proposal of the original DE algorithm [13] included the three main control parameters of the algorithm: the population size N P , scaling factor F, and crossover rate C r . Most of the further studies on DE were focused on either developing methods to control or adapt these parameters or proposing new mutation strategies and modifying the algorithm scheme. Some of the earliest theoretical works [17,18] considered the convergence properties of the DE and their connection to parameter settings on a set of benchmark problems, resulting in a set of recommendations for parameter values mainly in connection with the current population variance [19,20].
Although these findings were important for the field, the adaptive DE variants appeared to have more influence on the field. Starting with one of the earlier studies [21], which proposed the popular jDE approach, which proposed the random generation of F and C r with remembering successful values, there was a set of other adaptive DE variants proposed, including the SDE [22], SaDE (Self-adaptive Differential Evolution) [23], and the already mentioned JADE algorithm [15]. All these methods are based on sampling parameter values with random values distributed either with a uniform, normal, or Cauchy distribution. The most efficient scheme used guided adaptation, where more efficient parameter values were determined based on fitness values, like in the JADE [15] or SHADE [8] approach. In [24], the scale factor local search approach was proposed, where unlike the jDE algorithm, where the F parameter is generated randomly, the local search algorithms are used, such as the golden ratio or hill-climb. The idea of such a local search is similar to the steepest descent method, where the step size is optimized to get the best possible improvement.
In [25], the parameter control strategy was proposed, where the scaling factor was linearly decreased with the iteration number, allowing achieving better results compared to other methods. Several studies [26,27] have proposed different approach to the parameter adaptation, where a pool of fixed parameter values is maintained, and some of them are randomly chosen during the search process. A number of papers described the adaptive selection of mutation strategies used in the DE algorithm, for example adaptive DE with four mutation strategies [28] or a multi-population ensemble of mutations strategies [29]. In [30], the multi-population DE with scaling factor inheritance was proposed, where efficient solutions from one population were transferred to another together with the scaling factor.
Despite the demonstrated efficiency of these approaches, the annual competitions on bound-constrained single-objective optimization organized within the Congress on Evolutionary Computation demonstrate that the most promising approaches are DE with adaptive parameter values, such as SHADE [31] or jDE variants, or hybrid methods, such as LSHADE-SPACMA [32], where L-SHADE was hybridized with CMA-ES (Covariance-Matrix Adaptation Evolution Strategy) with dynamic resource redistribution, or [33], where the super-fit individual generated with CMA-ES was included in the population, although the multi-strategy methods also demonstrate highly competitive results [34].
Several surveys were published discussing the problem of parameter adaptation in DE and considering the existing methods’ advantages and drawbacks; a detailed taxonomy was presented in [35]. According to the results of this study, one of the most efficient schemes, the success-history adaptation, was initially introduced in the SHADE algorithm, the winner of CEC 2014 competition, and it represents the parameter tuning techniques introduced in JADE. SHADE maintains H memory cells, each keeping a pair ( M F , h , M C r , h ) , where h is the current memory index. For every mutation and crossover operation, new values are generated as follows:
F = r a n d c ( M F , k , 0.1 ) , C r = r a n d n ( M C r , k , 0.1 ) ,
where r a n d c is a Cauchy distributed random value, r a n d n is normally distributed, and k is chosen in the range [ 1 , D ] for every individual. The memory cells are used as location parameters of the distributions, and if the generated F value is below zero, it is generated again; however, if F > 1 , then it is set to one; also if C r < 0 , then C r = 0 , and if C r > 1 , then C r = 1 .
The memory cells are updated using the values of F and C r , which allow generating solutions with better fitness values. The successful F and C r are stored in S F and S C r , as well as Δ f j = | f ( u j ) f ( x j ) | . The update is performed as follows:    
M F , k g + 1 = 0.5 ( M F , k g + m e a n w L ( F ) ) ,
M C r , k g + 1 = 0.5 ( M C r , k g + m e a n w L ( C r ) ) ,
where g is the current generation number and m e a n w L is the weighted Lehmer mean calculated as follows:
m e a n w L = j = 1 | S | w j S j 2 j = 1 | S | w j S j ,
where w j = Δ f j k = 1 | S | Δ f k and S is either S C r or S F .
Except the success-history adaptation, there are several other popular parameter adaptation methods, for example the jDE algorithm [21], which has been shown to be very efficient in solving complex optimization problems if a large computational resource is given [36]. In jDE, the following mechanism was proposed: for every individual in the population, a memory cell keeping a pair of F and C r values is maintained. Initially, the values for all cells are set to M F , i = 0.5 and M C r , i = 0.9 , i = 1 , , N P , and new values are generated before the mutation and crossover operations:
F i , g + 1 = r a n d o m ( F l , F u ) , if r a n d o m ( 0 , 1 ) < τ 1 F i , g , otherwise .
In a similar manner, the C r values are generated:
C r i , g + 1 = r a n d o m ( 0 , 1 ) , if r a n d o m ( 0 , 1 ) < τ 2 C r i , g , otherwise .
The parameters τ 1 and τ 2 control the frequency of the F and C r change and usually are set to 0.1. If the offspring generated with these parameter values have better fitness compared to the parent, then these F and C r values are saved in the corresponding memory cells. The resulting jDE algorithm appeared to be highly competitive despite its simplicity, and its modifications jDE100 [37] and j2020 [36] have shown promising results during the CEC 2019 and CEC 2020 competitions on single-objective bound-constrained optimization with a large computational resource.
One of the three main parameters of DE is the population size N P . There have been several important studies on the DE population size, including [38]. Although several population size tuning algorithms exist, for example the structured population size reduction introduced in the SPSRDEMMS (Structured Population Size Redution Differential Evolution with Multiple Mutation Strategies) algorithm [39], the idea proposed in L-SHADE, the Linear Population Size Reduction (LPSR) is one of the most widely used. In LPSR, one or several individuals to be deleted is determined, and the worst in terms of fitness are removed at the end of each generation g. The new population size is calculated as follows:
N P g = r o u n d ( N P m i n N P m a x N F E m a x N F E ) + N P m a x
where N P m i n = 4 is the minimal population size and N P m a x is the initial population size. The population size reduction idea is to allow a wide search at the beginning to cover as much search space as possible, but later decrease the population size for better convergence to the best located optimum. The archive size N A is decreased in the same way as N P .

2.3. Proposed Approach: Linear Bias Reduction

One of the problems in parameter adaptation for DE is the existence of structural bias in the parameter adaptation process. For example, the problem of reaching different areas of the search space by DE was discussed in [40]. The JADE and SHADE algorithms used the Lehmer mean [41] instead of the classical arithmetical mean to avoid the structural bias in terms of parameter values, because it is much easier for the algorithm to generate better solutions with smaller F and C r values than larger F and C r , because smaller parameter values result in more local search, and thus easy improvements. However, such greedy behavior could only be beneficial from a short-term perspective, and only for local search. In [12], the generalized Lehmer mean was proposed, where a larger p parameter in the mean was allowed:
m e a n p , m , w L ( S ) = j = 1 | S | w j S j p j = 1 | S | w j S j m
The provided equation in its general form defines a group of means, and the values of p and m control the bias towards smaller or larger values. For example, m e a n 0 , 1 , w L ( S ) is the harmonic mean, m e a n 0.5 , 1 , w L ( S ) the geometric mean, m e a n 1 , 1 , w L ( S ) the arithmetic mean, and m e a n 2 , 1 , w L ( S ) the contraharmonic mean. Any biased mean could be obtained by setting the parameter values, and the graph of possible m e a n p , m , w L is presented in Figure 1.
The behavior of the parameter adaptation in DE could be significantly influenced by the newly generated ( M F , h , M C r , h ) values, so that the setting of p and m may influence the sampled F and C r  values, thus leading to a more explorative search, with large p and small m, or a more exploitative search with small p and large m. Further, in this study, the fixed setting of m = p 1 will be used, as changing only the p parameter in the Lehmer mean is sufficient to control the biased mean.
Inspired by the idea of linear population size reduction, the Linear Bias Reduction (LBR) technique is proposed. As long as in most cases, the computational resource limit is known, at the beginning of the search process, more exploration is required to find the most promising areas of the search space, while at the end, fast convergence is required. The linear bias reduction starts with a large p m a x value and gradually decreases the Lehmer mean parameter down to p m i n as follows:
p g = p m i n + ( p m a x p m i n ) N F E m a x N F E N F E m a x
In this study, the lower limit is set to p m i n = 1 , which corresponds to the arithmetic mean, while the upper limit p m a x is changed. The idea of LBR is inherited from the linear parameter reduction scheme L-SHADE algorithm [10] and other studies, such as [25]. The pseudocode of the L-SHADE-LBR is presented in Algorithm 1.
Algorithm 1 L-SHADE-LBR.
1:
Set N P m a x , N F E = 0 , N P = N P m a x D, N F E m a x ,
2:
H = 5 , A = , M F , r = 0.5 , M C r , r = 0.8 , r = 1 H
3:
M F , H + 1 = 0.9 , M C r , H + 1 = 0.9 , N A = 0.7 N P if jSO
4:
g = 0
5:
Initialize population P 0 = ( x 1 , j , , x N P , j ) randomly
6:
while N F E < N F E m a x do
7:
     S F = , S C r = , Δ f =
8:
    for i = 1 to N P do
9:
        Current memory index r = r a n d i n t [ 1 , H + 1 ]
10:
         C r = r a n d n ( M C r , r , 0.1 )
11:
         C r = m i n ( 1 , m a x ( 0 , C r ) )
12:
        repeat
13:
            F = r a n d c ( M F , r , 0.1 )
14:
        until F 0
15:
         F = m i n ( 1 , F )
16:
        Get F w if jSO
17:
        Limit C r if jSO
18:
        repeat
19:
           Randomly generate r 1 from p % best
20:
           Randomly generate r 2 from 1 , N P g
21:
           Randomly generate r 3 from 1 , N A g + N P g
22:
        until i r 1 r 2 r 3
23:
        for j=1 to D do
24:
           if jSO then
25:
                v i , j = x i , j + F w ( x r 1 , j x i , j ) + F ( x r 2 , j x r 3 , j )
26:
           else
27:
                v i , j = x i , j + F ( x r 1 , j x i , j ) + F ( x r 2 , j x r 3 , j )
28:
           end if
29:
        end for
30:
        Get u i with C r , and calculate f ( u i )
31:
        if f ( u i ) < f ( x i ) then
32:
            x i A , x i = u i , F S F , C r S C r , f ( x i ) f ( u i ) Δ f
33:
        end if
34:
    end for
35:
    Update population size N P g + 1
36:
    Recalculate maximum archive size N A g + 1
37:
    Recalculate LBR p g parameter
38:
    if | A | > N A g + 1 then
39:
        Remove random individuals from the archive
40:
    end if
41:
    if N P g > N P g + 1 then
42:
        Remove worst individuals from the population
43:
    end if
44:
    if S F and S C r then
45:
        if DISH then
46:
           Update M F , k and M C r , k with LBR and Δ f as the weight
47:
        else
48:
           Update M F , k and M C r , k with LBR and distance as the weight
49:
        end if
50:
    end if
51:
     k = k + 1 , g = g + 1
52:
    if k > H then
53:
         k = 1
54:
    end if
55:
end while
56:
Return best solution x b e s t
The next section contains the experimental setup and results of the LBR application to state-of-the-art DE algorithms on the CEC 2017 and 2020 competitions benchmark problems.

3. Results

The evaluation of the linear bias reduction approach was performed on two sets of benchmark problems, in particular the CEC 2017 [5] and CEC 2020 [42] competitions on bound-constrained single-objective optimization. These two sets of test problems were chosen because they represent two different scenarios: the CEC 2017 tests algorithms on 30 different 10, 30, 50, and 100D problems with a limited computational resource, and CEC 2020 has only 10 problems in 5, 10, 15, and 20 dimensions, but allows exponential growth of the computational resource available. For CEC 2017 problems, the resource was set to 10,000D fitness evaluations, and for CEC 2020 it was set to 50,000, 10 6 , 3 × 10 6 , and 10 7 for 5 D , 10 D , 15 D , and 20 D respectively. All test problems were shifted and rotated to avoid the cheating of the algorithms.
The LBR modification was tested with the L-SHADE, jSO, and DISH algorithms, and the following parameters were set for all algorithms: initial population size N P m a x = 25 · log ( D ) · D for CEC 2017, N P m a x = 30 · D · D for CEC 2020, final population size N P m i n = 4 , p b e s t = 0.17 , archive size N A = N P , number of memory cells H = 5 . The initial value of p m a x in LBR was set to 40 for F update and 32 for C r update; these values were chosen according to the experimental results reported in [12]. The algorithm was implemented in C++ with GCC and run on a PC with Ubuntu 20.04, an Intel Core i7 8700 k processor, and 32 GB RAM, with the results’ post-processing performed using Python 3.6.
The comparison of the results was performed for L-SHADE, jSO, and DISH with and without LBR. In the case that LBR was not used, the p-values in the Lehmer mean were set to two. The Mann–Whitney statistical test with normal approximation, tie-breaking, and the significance level set to 0.01 was applied to identify the difference between algorithm modifications. Table 1 presents the results of comparing the basic L-SHADE with L-SHADE-LBR on all dimensions of the CEC 2017 test problems. The values in the table show the number of wins (+), non-significant differences (=), and losses (−) summed over all functions.
The results presented in Table 1 show that the LBR gives up to 16 improvements for L-SHADE on 100-dimensional functions and 17 improvements for 50-dimensional functions. The jSO and DISH algorithms, which perform better than L-SHADE, and LBR give smaller improvements for these algorithms. Furthermore, the greater the dimension of the test problems, the more improvements are achieved, i.e., for example, for 10-dimensional problems, the LBR does not deliver any improvements, and for L-SHADE and jSO, there are even several performance deteriorations. Most of the improvements are for multimodal, hybrid, and composition functions.
Table 2, Table 3, Table 4 and Table 5 contain the results for all CEC 2017 functions for all dimensions. The results presented in Table 2, Table 3, Table 4 and Table 5 show that the linear bias reduction has the largest influence on 50- and 100-dimensional functions. Best values for each function are highlighted with bold. Table 6 contains the Mann–Whitney tests for the CEC 2020 benchmark problems.
The results presented in Table 6 shows that the LBR allows up to two improvements for the L-SHADE algorithm for 10 D , 15 D , and 20 D , and one decreased performance. The improvements were mainly for functions f5 and f6, i.e., hybrid functions. For the jSO algorithm, there were up to three improvements for 15D, and for the DISH algorithm, which also included jSO adaptations, there were up to two improvements, which were observed for f5, f6, and f7.
Table 7, Table 8, Table 9 and Table 10 contain the mean and standard deviation results for all functions and algorithms for 5-, 10-, 15-, and 20-dimensional problems. Best values for each function are highlighted with bold.

4. Discussion

Every evolutionary algorithm should be designed and tuned according to the specific purpose; for example, if the goal is to find the best possible function value in a very limited time period, the convergence speed of the algorithm is crucial; however, if the computational resource is relatively large, then the explorative properties of the algorithm become more important. Obviously, both the fast convergence and the high explorative properties are important, and the algorithm that is able to combine them and make use of both could be recommended for real-world usage.
The linear bias reduction scheme proposed in this study follows this idea: the experiments showed that it improves the performance of several state-of-the-art non-hybrid DE variants on both the CEC 2017 and CEC 2020 problems due to larger parameter values sampled at the beginning of the search, which promote exploration, and smaller closer to the end of the search. The LBR approach does not require any significant additional computational overhead and could be easily implemented in most modern DE-based optimizers, improving their performance.
Although the idea of exploring at the beginning of the search and exploiting at the end is not new to the field of evolutionary algorithms, it could be implemented in many ways, for example with linear population size reduction; however, as this study shows, not all of them are discovered and properly studied. The idea of LBR is relatively simple and could be used for other approaches with corresponding modifications.
For example, in a genetic algorithm for tuning mutation probability, a similar success-history-based adaptation could be used, as was shown in [43]. Implementing the LBR proposed in this study for this version of GA could significantly improve the performance. In general, any parameter adaptation mechanism used in any EA or SI algorithm, where the tuning of numerical parameters is required, could be modified with LBR with minor changes depending on the algorithm. The implementation of LBR for multiobjective algorithms is also possible; however, to the best of our knowledge, there have been only a few numerical parameter tuning variants proposed for multiobjective optimization. As long as LBR only relies on the amount of computational resource, if the resource is known in advance, it could be incorporated into the algorithm.

5. Conclusions

This paper proposes a parameter tuning modification for the differential evolution algorithm called linear bias reduction, which allows more exploration at the beginning of the search and faster convergence to the best located optimum at the end. The LBR is a relatively simple technique, which could be used not only for DE, but for other approaches. The performed experiments proved its efficiency, especially for complex hybrid functions, and at the same time, there were only a few performance losses for more simple problems. The further directions of study may include, but are not limited to experimenting with non-linear bias reduction, bias increase, and various initial and final parameter values in LBR, as well as implementing LBR-like techniques for other evolutionary or swarm intelligence algorithms.

Author Contributions

Conceptualization, V.S. and S.A.; methodology, V.S., S.A., and E.S.; software, V.S. and E.S.; validation, V.S., S.A., and E.S.; formal analysis, S.A.; investigation, V.S.; resources, E.S. and V.S.; data curation, E.S.; writing, original draft preparation, V.S. and S.A.; writing, review and editing, V.S.; visualization, S.A. and V.S.; supervision, E.S.; project administration, E.S. and S.A.; funding acquisition, S.A. and V.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Science and Higher Education of the Russian Federation within the limits of state contract No. FEFE-2020-0013.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
EAEvolutionary Algorithm
DEDifferential Evolution
L-SHADELinear population size reduction Success History-based Adaptive Differential Evolution
LBRLinear Bias Reduction
CECCongress on Evolutionary Computation

References

  1. Aleti, A.; Moser, I. A Systematic Literature Review of Adaptive Parameter Control Methods for Evolutionary Algorithms. ACM Comput. Surv. 2016, 49, 1–35. [Google Scholar] [CrossRef]
  2. Eiben, A.E.; Selmar, K.S. Evolutionary Algorithm Parameters and Methods to Tune Them. In Autonomous Search; Springer: Berlin/Heidelberg, Germany, 2012; pp. 15–36. [Google Scholar]
  3. Price, K.; Storn, R.M.; Lampinen, J.A. Differential Evolution: A Practical Approach to Global Optimization; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  4. Das, S.; Suganthan, P.N. Differential Evolution: A survey of the state-of-the-art. IEEE Trans. Evol. Comput. 2011, 15, 4–31. [Google Scholar] [CrossRef]
  5. Awad, N.H.; Ali, M.Z.; Liang, J.J.; Qu, B.Y.; Suganthan, P.N. Problem Definitions and Evaluation Criteria for the CEC 2017 Special Session and Competition on Single Objective Bound Constrained Real-Parameter Numerical Optimization; Technial Report; Nanyang Technological University: Singapore, 2016. [Google Scholar]
  6. Neri, F.; Tirronen, V. Recent advances in Differential Evolution: A survey and experimental analysis. Artif. Intell. Rev. 2010, 33, 61–106. [Google Scholar] [CrossRef]
  7. Das, S.; Mullick, S.S.; Suganthan, P.N. Recent advances in Differential Evolution–an updated survey. Swarm Evol. Comput. 2016, 15, 1–30. [Google Scholar] [CrossRef]
  8. Tanabe, R.; Fukunaga, A.S. Success-history based parameter adaptation for differential evolution. In Proceedings of the IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 71–78. [Google Scholar]
  9. Viktorin, A.; Senkerik, R.; Pluhacek, M.; Kadavy, T.; Zamuda, A. Distance based parameter adaptation for Success-History based Differential Evolution. Swarm Evol. Comput. 2019, 50, 100462. [Google Scholar] [CrossRef]
  10. Tanabe, R.; Fukunaga, A.S. Improving the search performance of SHADE using linear population size reduction. In Proceedings of the IEEE Congress on Evolutionary Computation, Beijing, China, 6–11 July 2014; pp. 1658–1665. [Google Scholar]
  11. Brest, J.; Maučec, M.S.; Boškovic, B. Single objective real-parameter optimization algorithm jSO. In Proceedings of the IEEE Congress on Evolutionary Computation, San Sebastian, Spain, 5–8 June 2017; pp. 1311–1318. [Google Scholar]
  12. Stanovov, V.; Akhmedova, S.; Semenkin, E.; Semenkina, M. Generalized Lehmer Mean for Success History based Adaptive Differential Evolution. In Proceedings of the IJCCI 2019: 11th International Joint Conference on Computational Intelligence, Vienna, Austria, 17–19 September 2019; pp. 93–100. [Google Scholar]
  13. Storn, R.M.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  14. Opara, K.R.; Arabas, J. Differential Evolution: A survey of theoretical analyses. Swarm Evol. Comput. 2019, 44, 546–558. [Google Scholar] [CrossRef]
  15. Zhang, J.; Snaderson, A. JADE: Adaptive Differential Evolution With Optional External Archive. IEEE Trans. Evol. Comput. 2009, 13, 945–958. [Google Scholar] [CrossRef]
  16. Stanovov, V.; Akhmedova, S.; Semenkin, E. Selective Pressure Strategy in differential evolution: Exploitation improvement in solving global optimization problems. Swarm Evol. Comput. 2019, 50, 100463. [Google Scholar] [CrossRef]
  17. Zaharie, D. On the explorative power of differential evolution. In Proceedings of the 3rd International Workshop Symbolic Numerical Algorithms Scientific Computing, Timişoara, Romania, 2–5 October 2001. [Google Scholar]
  18. Zaharie, D. Critical values for the control parameters of differential evolution algorithms. In Proceedings of the 8th International Mendel Conference Soft Computing, Brno, Czech Republic, 8–10 June 2002; pp. 62–67. [Google Scholar]
  19. Zaharie, D. Parameter adaptation in differential evolution by controlling the population diversity. In Proceedings of the 4th International Workshop Symbolic Numeric Algorithms Scientific Computing, Timisoara, Romania, 9–12 October 2002; pp. 385–397. [Google Scholar]
  20. Zaharie, D. Statistical properties of differential evolution and related random search algorithms. In Proceedings of the International Conference on Computational Statistics, Kraków, Poland, 23–25 June 2008; pp. 473–485. [Google Scholar]
  21. Brest, J.; Greiner, S.; Boškovic, B.; Mernik, M.; Žumer, V. Self-adapting control parameters in differential evolution: A comparative study on numerical benchmark problems. IEEE Trans. Evol. Comput. 2006, 10, 646–657. [Google Scholar] [CrossRef]
  22. Omran, M.G.H.; Salman, A.; Engelbrecht, A.P. Self-adaptive differential evolution. Comput. Intell. Secur. Lect. Notes Artif. Intell. 2005, 3801, 192–199. [Google Scholar]
  23. Qin, A.K.; Huang, V.L.; Suganthan, P.N. Differential evolution algorithm with strategy adaptation for global numerical optimization. IEEE Trans. Evol. Comput. 2009, 13, 398–417. [Google Scholar] [CrossRef]
  24. Neri, F.; Tirronen, V. Scale factor local search in differential evolution. Memet. Comput. 2009, 1, 153–171. [Google Scholar] [CrossRef]
  25. Das, S.; Konar, A.; Chakraborty, U.K. Two improved differential evolution schemes for faster global search. In Proceedings of the 2005 Conferenceon Genetic and Evolutionary Computation, Washington, DC, USA, 25–29 June 2005; pp. 991–998. [Google Scholar]
  26. Elsayed, S.M.; Sarker, R.A.; Ray, T. Differential evolution with automatic parameter configuration for solving the CEC2013 competition on real-parameter optimization. In Proceedings of the IEEE Congresson Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 1932–1937. [Google Scholar]
  27. Sarker, R.A.; Elsayed, S.M.; Ray, T. Differential evolution with dynamic parameters selection for optimization problems. IEEE Trans. Evol. Comput. 2014, 18, 689–707. [Google Scholar] [CrossRef]
  28. Gong, W.; Cai, Z.; Ling, C.X.; Li, H. Enhanced Differential Evolution With Adaptive Strategies for Numerical Optimization. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2011, 41, 397–413. [Google Scholar] [CrossRef] [PubMed]
  29. Wu, G.; Mallipeddi, R.; Suganthan, P.N.; Wang, R.; Chen, H. Differential evolution with multi-population based ensemble of mutation strategies. Inf. Sci. 2016, 329, 329–345. [Google Scholar] [CrossRef]
  30. Weber, M.; Tirronen, V.; Neri, F. Scale factor inheritance mechanism in distributed differential evolution. Soft Comput. 2010, 14, 1187–1207. [Google Scholar] [CrossRef]
  31. Awad, N.H.; Ali, M.Z.; Suganthan, P.N. Ensemble sinusoidal differential covariance matrix adaptation with euclidean neighborhood for solving CEC2017 benchmark problems. In Proceedings of the IEEE Congress on Evolutionary Computation (CEC), San Sebastian, Spain, 5–8 June 2017; pp. 372–379. [Google Scholar]
  32. Mohamed, A.W.; Hadi, A.A.; Fattouh, A.M.; Jambi, K.M. LSHADE with semi-parameter adaptation hybrid with CMA-ES for solving CEC 2017 benchmark problems. In Proceedings of the IEEE Congress on Evolutionary Computation (CEC), San Sebastian, Spain, 5–8 June 2017; pp. 145–152. [Google Scholar]
  33. Caraffini, F.; Neri, F.; Cheng, J.; Zhang, G.; Picinali, L.; Iacca, G.; Mininno, E. Super-fit Multicriteria Adaptive Differential Evolution. In Proceedings of the IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013. [Google Scholar]
  34. Sallam, K.M.; Elsayed, S.M.; Chakrabortty, R.K.; Ryan, M.J. Improved Multi-operator Differential Evolution Algorithm for Solving Unconstrained Problems. In Proceedings of the IEEE Congress on Evolutionary Computation (CEC), Glasgow, UK, 19–24 July 2020. [Google Scholar]
  35. Al-Dabbagh, R.D.; Neri, F.; Idris, N.; Baba, M.S. Algorithmic design issues in adaptive differential evolution schemes: Review and taxonomy. Swarm Evol. Comput. 2018, 44, 284–311. [Google Scholar] [CrossRef]
  36. Brest, J.; Maučec, M.S.; Boškovic, B. Differential Evolution Algorithm for Single Objective Bound-Constrained Optimization: Algorithm j2020. In Proceedings of the 2020 IEEE Congress on Evolutionary Computation (CEC), Glasgow, UK, 19–24 July 2020. [Google Scholar]
  37. Brest, J.; Maučec, M.S.; Boškovic, B. The 100-Digit Challenge: Algorithm jDE100. In Proceedings of the 2019 IEEE Congress on Evolutionary Computation (CEC), Wellington, New Zealand, 10–13 June 2019; pp. 19–26. [Google Scholar]
  38. Piotrowski, A. Review of Differential Evolution population size. Swarm Evol. Comput. 2019, 32, 1–24. [Google Scholar] [CrossRef]
  39. Zamuda, A.; Brest, J.; Mezura-Montes, E. Structured Population Size Reduction Differential Evolution with Multiple Mutation Strategies on CEC 2013 real parameter optimization. In Proceedings of the 2013 IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 1925–1931. [Google Scholar]
  40. Caraffini, F.; Kononova, A.V.; Corne, D. Infeasibility and structural bias in differential evolution. Inf. Sci. 2019, 496, 161–179. [Google Scholar] [CrossRef] [Green Version]
  41. Bullen, P.S. Handbook of Means and Their Inequalities; Springer: Dordrecht, The Netherlands, 2003. [Google Scholar]
  42. Yue, C.T.; Price, K.V.; Suganthan, P.N.; Liang, J.J.; Ali, M.Z.; Qu, B.Y.; Awad, N.H.; Biswas, P. Problem Definitions and Evaluation Criteria for the CEC 2020 Special Session and Competition on Single Objective Bound Constrained Numerical Optimization; Technial Report; Nanyang Technological University: Singapore, 2019. [Google Scholar]
  43. Stanovov, V.; Akhmedova, S.; Semenkin, E. Genetic Algorithm with Success History based Parameter Adaptation. In Proceedings of the IJCCI 2019: 11th International Joint Conference on Computational Intelligence, Vienna, Austria, 17–19 September 2019; pp. 180–187. [Google Scholar]
Figure 1. Generalized weighted Lehmer mean depending on parameters p and m.
Figure 1. Generalized weighted Lehmer mean depending on parameters p and m.
Algorithms 13 00283 g001
Table 1. Statistical comparison of algorithms with and without LBR, Mann–Whitney test, CEC 2017 problems.
Table 1. Statistical comparison of algorithms with and without LBR, Mann–Whitney test, CEC 2017 problems.
Algorithm 10 D 30 D 50 D 100 D
L-SHADE vs. L-SHADE-LBR0+/27 = /3−8+/21 = /1−17+/12 = /1−16+/10 = /4−
jSO vs. jSO-LBR0+/29 = /1−5+/24 = /1−10+/20 = /0−12+/16 = /2−
DISH vs. DISH-LBR0+/30 = /0−2+/28 = /0−4+/25 = /1−9+/16 = /5−
Table 2. Results for CEC 2017 problems, 10 D .
Table 2. Results for CEC 2017 problems, 10 D .
FunctionL-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f20.0 ± 0.00.0 ± 0.00.0 ± 0.0
f30.0 ± 0.00.0 ± 0.00.0 ± 0.0
f40.0 ± 0.00.0 ± 0.00.0 ± 0.0
f51.466332 × 10 0 ± 7.966326 × 10 1 1.820984 × 10 0 ± 8.296894 × 10 1 2.028990 × 10 0 ± 9.021237 × 10 1
f60.0 ± 0.00.0 ± 0.00.0 ± 0.0
f71.165531 × 10 1 ± 4.999618 × 10 1 1.206354 × 10 1 ± 5.092644 × 10 1 1.202930 × 10 1 ± 5.555746 × 10 1
f81.338034 × 10 0 ± 8.316749 × 10 1 1.878932 × 10 0 ± 8.513805 × 10 1 1.853360 × 10 0 ± 1.042962 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f102.224059 × 10 1 ± 4.032137 × 10 1 1.822742 × 10 1 ± 3.227775 × 10 1 3.377682 × 10 1 ± 5.852428 × 10 1
f111.181580 × 10 1 ± 4.216185 × 10 1 1.516416 × 10 1 ± 3.630057 × 10 1 0.0 ± 0.0
f121.526398 × 10 1 ± 3.919816 × 10 1 2.734428 × 10 1 ± 1.723809 × 10 1 1.673563 × 10 1 ± 4.110506 × 10 1
f132.441962 × 10 0 ± 2.352857 × 10 0 2.474151 × 10 0 ± 2.390764 × 10 0 3.043588 × 10 0 ± 2.362332 × 10 0
f144.638013 × 10 2 ± 1.575738 × 10 1 4.747525 × 10 2 ± 1.948063 × 10 1 3.902250 × 10 2 ± 1.931283 × 10 1
f151.346569 × 10 1 ± 1.999631 × 10 1 1.225807 × 10 1 ± 1.933628 × 10 1 3.401098 × 10 1 ± 2.058195 × 10 1
f164.510030 × 10 1 ± 1.916985 × 10 1 3.580149 × 10 1 ± 1.948515 × 10 1 5.429780 × 10 1 ± 3.146415 × 10 1
f171.572282 × 10 1 ± 1.571081 × 10 1 1.409888 × 10 1 ± 1.729228 × 10 1 9.310686 × 10 1 ± 2.782305 × 10 0
f182.218042 × 10 1 ± 1.984943 × 10 1 1.422543 × 10 1 ± 1.699201 × 10 1 1.992204 × 10 1 ± 1.953476 × 10 1
f199.994697 × 10 3 ± 1.044910 × 10 2 1.164565 × 10 2 ± 1.284938 × 10 2 1.266488 × 10 2 ± 1.275577 × 10 2
f206.121039 × 10 3 ± 4.328228 × 10 2 6.121039 × 10 3 ± 4.328228 × 10 2 3.672626 × 10 1 ± 1.476685 × 10 1
f211.450680 × 10 2 ± 5.038353 × 10 1 1.506000 × 10 2 ± 5.160565 × 10 1 1.302909 × 10 2 ± 4.692892 × 10 1
f221.000056 × 10 2 ± 3.965344 × 10 2 1.000078 × 10 2 ± 5.532071 × 10 2 1.000000 × 10 2 ± 0.0
f233.010949 × 10 2 ± 1.487577 × 10 0 3.024831 × 10 2 ± 1.330968 × 10 0 3.010731 × 10 2 ± 1.481965 × 10 0
f242.592233 × 10 2 ± 1.036631 × 10 2 2.666515 × 10 2 ± 1.025191 × 10 2 2.748786 × 10 2 ± 9.709894 × 10 1
f254.237633 × 10 2 ± 2.248348 × 10 1 4.095063 × 10 2 ± 1.979616 × 10 1 4.095055 × 10 2 ± 1.979735 × 10 1
f263.000000 × 10 2 ± 0.02.941176 × 10 2 ± 4.159452 × 10 1 3.000000 × 10 2 ± 0.0
f273.893875 × 10 2 ± 2.231327 × 10 1 3.892298 × 10 2 ± 5.642233 × 10 1 3.894176 × 10 2 ± 2.032786 × 10 1
f284.231326 × 10 2 ± 1.442053 × 10 2 3.578377 × 10 2 ± 1.172708 × 10 2 3.294695 × 10 2 ± 8.948919 × 10 1
f292.332684 × 10 2 ± 2.297889 × 10 0 2.334880 × 10 2 ± 1.961177 × 10 0 2.347158 × 10 2 ± 3.342209 × 10 0
f303.244280 × 10 4 ± 1.586213 × 10 5 3.945063 × 10 2 ± 2.293938 × 10 2 1.641868 × 10 4 ± 1.133010 × 10 5
FunctionL-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f20.0 ± 0.00.0 ± 0.00.0 ± 0.0
f30.0 ± 0.00.0 ± 0.00.0 ± 0.0
f40.0 ± 0.00.0 ± 0.00.0 ± 0.0
f51.931394 × 10 0 ± 8.221626 × 10 1 1.833883 × 10 0 ± 8.894004 × 10 1 1.775343 × 10 0 ± 6.914080 × 10 1
f60.0 ± 0.00.0 ± 0.00.0 ± 0.0
f71.208296 × 10 1 ± 5.567794 × 10 1 1.216231 × 10 1 ± 5.829416 × 10 1 1.201511 × 10 1 ± 6.008566 × 10 1
f82.185022 × 10 0 ± 8.815939 × 10 1 1.853397 × 10 0 ± 9.033848 × 10 1 2.263072 × 10 0 ± 9.046168 × 10 1
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f103.481015 × 10 1 ± 5.336775 × 10 1 5.487013 × 10 1 ± 8.016722 × 10 1 2.936753 × 10 1 ± 5.244595 × 10 1
f110.0 ± 0.00.0 ± 0.00.0 ± 0.0
f127.398462 × 10 0 ± 2.818943 × 10 1 1.674478 × 10 1 ± 4.117527 × 10 1 2.640647 × 10 0 ± 1.640653 × 10 1
f132.803025 × 10 0 ± 2.437121 × 10 0 2.675462 × 10 0 ± 2.342831 × 10 0 3.290499 × 10 0 ± 2.290796 × 10 0
f141.560720 × 10 1 ± 3.618380 × 10 1 1.462447 × 10 1 ± 3.453513 × 10 1 1.170540 × 10 1 ± 3.205656 × 10 1
f152.393091 × 10 1 ± 2.081327 × 10 1 2.935211 × 10 1 ± 2.028156 × 10 1 2.354370 × 10 1 ± 2.096280 × 10 1
f167.189243 × 10 1 ± 3.081263 × 10 1 5.107098 × 10 1 ± 2.956673 × 10 1 5.481845 × 10 1 ± 2.890193 × 10 1
f174.551964 × 10 1 ± 3.136878 × 10 1 5.855833 × 10 1 ± 3.998470 × 10 1 5.405732 × 10 1 ± 3.427532 × 10 1
f182.217356 × 10 1 ± 2.083894 × 10 1 5.850836 × 10 1 ± 2.756645 × 10 0 2.256478 × 10 1 ± 2.133737 × 10 1
f191.122558 × 10 2 ± 2.105432 × 10 2 1.145965 × 10 2 ± 1.292694 × 10 2 8.007165 × 10 3 ± 1.032611 × 10 2
f208.095347 × 10 1 ± 2.808347 × 10 0 3.672626 × 10 1 ± 1.716079 × 10 1 3.978678 × 10 1 ± 1.524131 × 10 1
f211.323348 × 10 2 ± 4.782543 × 10 1 1.363186 × 10 2 ± 4.917879 × 10 1 1.444566 × 10 2 ± 5.104734 × 10 1
f221.000079 × 10 2 ± 5.559800 × 10 2 1.000068 × 10 2 ± 4.797234 × 10 2 1.000268 × 10 2 ± 9.379501 × 10 2
f233.010237 × 10 2 ± 1.490490 × 10 0 2.955855 × 10 2 ± 4.183554 × 10 1 3.013980 × 10 2 ± 1.570153 × 10 0
f243.018687 × 10 2 ± 7.386128 × 10 1 2.656853 × 10 2 ± 1.019472 × 10 2 2.654399 × 10 2 ± 1.018634 × 10 2
f254.015055 × 10 2 ± 1.221117 × 10 1 4.059703 × 10 2 ± 1.730538 × 10 1 4.068576 × 10 2 ± 1.801739 × 10 1
f263.000000 × 10 2 ± 0.03.000000 × 10 2 ± 0.03.000000 × 10 2 ± 0.0
f273.894635 × 10 2 ± 1.538980 × 10 1 3.893975 × 10 2 ± 2.171814 × 10 1 3.894477 × 10 2 ± 1.761875 × 10 1
f283.294695 × 10 2 ± 8.948919 × 10 1 3.361344 × 10 2 ± 9.902322 × 10 1 3.411471 × 10 2 ± 1.032896 × 10 2
f292.350070 × 10 2 ± 2.765177 × 10 0 2.345793 × 10 2 ± 3.214256 × 10 0 2.337357 × 10 2 ± 2.522189 × 10 0
f303.244093 × 10 4 ± 1.586216 × 10 5 1.641771 × 10 4 ± 1.133012 × 10 5 3.945173 × 10 2 ± 4.152092 × 10 2
Table 3. Results for CEC 2017 problems, 30 D .
Table 3. Results for CEC 2017 problems, 30 D .
FunctionL-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f20.0 ± 0.00.0 ± 0.00.0 ± 0.0
f30.0 ± 0.00.0 ± 0.00.0 ± 0.0
f45.432730 × 10 1 ± 1.761079 × 10 1 5.856160 × 10 1 ± 0.05.856160 × 10 1 ± 0.0
f57.093248 × 10 0 ± 1.514062 × 10 0 5.984748 × 10 0 ± 1.259837 × 10 0 8.566069 × 10 0 ± 1.414316 × 10 0
f61.610282 × 10 8 ± 4.409940 × 10 8 2.683804 × 10 9 ± 1.897736 × 10 8 1.125136 × 10 7 ± 3.080402 × 10 7
f73.643252 × 10 1 ± 1.163031 × 10 0 3.635001 × 10 1 ± 1.050653 × 10 0 3.901514 × 10 1 ± 1.897382 × 10 0
f87.636827 × 10 0 ± 1.421298 × 10 0 5.947331 × 10 0 ± 1.286650 × 10 0 8.606066 × 10 0 ± 1.668634 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f101.541919 × 10 3 ± 2.100733 × 10 2 1.588640 × 10 3 ± 2.124062 × 10 2 1.633723 × 10 3 ± 2.570738 × 10 2
f111.480573 × 10 1 ± 2.154669 × 10 1 1.286354 × 10 1 ± 2.135075 × 10 1 6.873639 × 10 0 ± 1.395100 × 10 1
f129.665733 × 10 2 ± 3.204188 × 10 2 2.131759 × 10 2 ± 1.341268 × 10 2 2.351108 × 10 2 ± 1.392522 × 10 2
f131.507205 × 10 1 ± 7.899218 × 10 0 1.747101 × 10 1 ± 5.664911 × 10 0 1.616138 × 10 1 ± 5.997028 × 10 0
f142.133530 × 10 1 ± 9.644383 × 10 1 2.136757 × 10 1 ± 3.272063 × 10 0 2.194751 × 10 1 ± 3.183398 × 10 0
f152.465098 × 10 0 ± 1.482018 × 10 0 1.524179 × 10 0 ± 9.985145 × 10 1 1.212446 × 10 0 ± 8.005678 × 10 1
f166.228463 × 10 1 ± 7.399308 × 10 1 3.557114 × 10 1 ± 3.589944 × 10 1 4.984834 × 10 1 ± 6.435235 × 10 1
f173.486476 × 10 1 ± 5.016896 × 10 0 3.316091 × 10 1 ± 5.201363 × 10 0 3.505441 × 10 1 ± 6.723407 × 10 0
f182.146325 × 10 1 ± 3.044085 × 10 0 2.077898 × 10 1 ± 4.073286 × 10 1 2.081425 × 10 1 ± 4.281106 × 10 1
f195.078805 × 10 0 ± 1.737644 × 10 0 5.462701 × 10 0 ± 1.914274 × 10 0 5.006098 × 10 0 ± 1.699874 × 10 0
f203.322552 × 10 1 ± 5.277263 × 10 0 3.011395 × 10 1 ± 4.949041 × 10 0 2.957005 × 10 1 ± 7.071608 × 10 0
f212.068744 × 10 2 ± 1.548801 × 10 0 2.063994 × 10 2 ± 1.477586 × 10 0 2.087112 × 10 2 ± 1.873729 × 10 0
f221.000000 × 10 2 ± 0.01.000000 × 10 2 ± 0.01.000000 × 10 2 ± 0.0
f233.485940 × 10 2 ± 2.713869 × 10 0 3.483330 × 10 2 ± 2.099438 × 10 0 3.509930 × 10 2 ± 2.924522 × 10 0
f244.237707 × 10 2 ± 1.651448 × 10 0 4.251463 × 10 2 ± 1.305455 × 10 0 4.262602 × 10 2 ± 2.262393 × 10 0
f253.867463 × 10 2 ± 2.090470 × 10 2 3.867043 × 10 2 ± 1.523778 × 10 2 3.866993 × 10 2 ± 9.419400 × 10 3
f269.245550 × 10 2 ± 4.373209 × 10 1 9.254841 × 10 2 ± 3.466222 × 10 1 9.327218 × 10 2 ± 3.297828 × 10 1
f275.016845 × 10 2 ± 5.896865 × 10 0 4.926597 × 10 2 ± 7.915101 × 10 0 4.981355 × 10 2 ± 7.186093 × 10 0
f283.114948 × 10 2 ± 3.531211 × 10 1 3.064948 × 10 2 ± 2.600802 × 10 1 3.217193 × 10 2 ± 4.403163 × 10 1
f294.353955 × 10 2 ± 1.046230 × 10 1 4.392384 × 10 2 ± 1.303416 × 10 1 4.348371 × 10 2 ± 1.639610 × 10 1
f301.994910 × 10 3 ± 5.095366 × 10 1 1.969638 × 10 3 ± 1.586864 × 10 1 1.975658 × 10 3 ± 3.868023 × 10 1
FunctionL-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f20.0 ± 0.00.0 ± 0.00.0 ± 0.0
f30.0 ± 0.00.0 ± 0.00.0 ± 0.0
f45.856160 × 10 1 ± 0.05.867054 × 10 1 ± 7.703027 × 10 1 5.856160 × 10 1 ± 0.0
f57.138668 × 10 0 ± 1.934874 × 10 0 8.172221 × 10 0 ± 1.624578 × 10 0 7.300936 × 10 0 ± 1.628259 × 10 0
f65.733510 × 10 8 ± 1.749795 × 10 7 2.952184 × 10 8 ± 5.629595 × 10 8 5.360755 × 10 8 ± 1.577789 × 10 7
f73.729442 × 10 1 ± 1.342658 × 10 0 3.843567 × 10 1 ± 1.603208 × 10 0 3.783948 × 10 1 ± 1.356820 × 10 0
f87.340093 × 10 0 ± 1.966085 × 10 0 8.462537 × 10 0 ± 2.018838 × 10 0 7.322912 × 10 0 ± 2.011191 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f101.568173 × 10 3 ± 2.895061 × 10 2 1.524144 × 10 3 ± 2.684248 × 10 2 1.568425 × 10 3 ± 2.629532 × 10 2
f114.027302 × 10 0 ± 8.265100 × 10 0 4.306909 × 10 0 ± 8.200147 × 10 0 7.472353 × 10 0 ± 1.565693 × 10 1
f128.106433 × 10 1 ± 7.149354 × 10 1 9.996226 × 10 1 ± 6.182176 × 10 1 6.730995 × 10 1 ± 6.684616 × 10 1
f131.552083 × 10 1 ± 6.133034 × 10 0 1.610696 × 10 1 ± 5.854376 × 10 0 1.684447 × 10 1 ± 4.239543 × 10 0
f142.171615 × 10 1 ± 1.195599 × 10 0 2.195247 × 10 1 ± 1.241087 × 10 0 2.203156 × 10 1 ± 1.276744 × 10 0
f151.184498 × 10 0 ± 7.500999 × 10 1 1.166798 × 10 0 ± 7.381401 × 10 1 1.213596 × 10 0 ± 9.026724 × 10 1
f165.466729 × 10 1 ± 7.181852 × 10 1 3.527649 × 10 1 ± 4.731054 × 10 1 3.623068 × 10 1 ± 5.130314 × 10 1
f173.449224 × 10 1 ± 7.695939 × 10 0 3.388604 × 10 1 ± 6.728233 × 10 0 3.355475 × 10 1 ± 7.485698 × 10 0
f182.084219 × 10 1 ± 3.720811 × 10 1 2.006454 × 10 1 ± 3.965811 × 10 0 2.045489 × 10 1 ± 2.704982 × 10 0
f194.817716 × 10 0 ± 1.898779 × 10 0 4.076852 × 10 0 ± 1.356127 × 10 0 4.425025 × 10 0 ± 1.967937 × 10 0
f202.972537 × 10 1 ± 8.035391 × 10 0 2.704956 × 10 1 ± 6.845223 × 10 0 2.679755 × 10 1 ± 7.130091 × 10 0
f212.083516 × 10 2 ± 1.807313 × 10 0 2.085552 × 10 2 ± 1.968052 × 10 0 2.080673 × 10 2 ± 2.171051 × 10 0
f221.000000 × 10 2 ± 0.01.000000 × 10 2 ± 0.01.000000 × 10 2 ± 0.0
f233.505319 × 10 2 ± 3.690104 × 10 0 3.507215 × 10 2 ± 3.276329 × 10 0 3.498206 × 10 2 ± 3.709504 × 10 0
f244.263549 × 10 2 ± 2.213716 × 10 0 4.259219 × 10 2 ± 1.872762 × 10 0 4.262288 × 10 2 ± 2.572837 × 10 0
f253.866960 × 10 2 ± 7.127660 × 10 3 3.866957 × 10 2 ± 5.585311 × 10 3 3.866962 × 10 2 ± 6.091826 × 10 3
f269.561221 × 10 2 ± 4.264016 × 10 1 9.313090 × 10 2 ± 4.626237 × 10 1 9.371625 × 10 2 ± 4.567299 × 10 1
f274.932655 × 10 2 ± 7.774858 × 10 0 4.932667 × 10 2 ± 7.143169 × 10 0 4.899617 × 10 2 ± 6.927315 × 10 0
f283.042600 × 10 2 ± 2.111257 × 10 1 3.085200 × 10 2 ± 2.924355 × 10 1 3.000000 × 10 2 ± 0.0
f294.348186 × 10 2 ± 1.955936 × 10 1 4.340111 × 10 2 ± 1.408097 × 10 1 4.350100 × 10 2 ± 1.643612 × 10 1
f301.969240 × 10 3 ± 9.962199 × 10 0 1.972825 × 10 3 ± 2.132789 × 10 1 1.966696 × 10 3 ± 9.911583 × 10 0
Table 4. Results for CEC 2020 problems, 50 D .
Table 4. Results for CEC 2020 problems, 50 D .
FunctionL-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f20.0 ± 0.00.0 ± 0.00.0 ± 0.0
f30.0 ± 0.00.0 ± 0.00.0 ± 0.0
f45.235449 × 10 1 ± 4.437343 × 10 1 5.987294 × 10 1 ± 4.656380 × 10 1 5.533534 × 10 1 ± 4.638192 × 10 1
f51.858622 × 10 1 ± 2.721102 × 10 0 1.090016 × 10 1 ± 2.044454 × 10 0 1.643115 × 10 1 ± 2.349557 × 10 0
f61.826999 × 10 4 ± 9.247796 × 10 4 3.321219 × 10 7 ± 7.189487 × 10 7 1.159359 × 10 6 ± 1.645192 × 10 6
f76.474578 × 10 1 ± 1.859098 × 10 0 6.125818 × 10 1 ± 1.330959 × 10 0 6.533251 × 10 1 ± 2.386148 × 10 0
f81.786065 × 10 1 ± 2.476184 × 10 0 1.059279 × 10 1 ± 1.794339 × 10 0 1.642296 × 10 1 ± 2.506320 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f103.145981 × 10 3 ± 3.519306 × 10 2 3.289575 × 10 3 ± 3.466515 × 10 2 3.287218 × 10 3 ± 3.405172 × 10 2
f114.479648 × 10 1 ± 7.262966 × 10 0 2.847299 × 10 1 ± 3.609160 × 10 0 2.765099 × 10 1 ± 3.060072 × 10 0
f122.211379 × 10 3 ± 5.051756 × 10 2 1.791286 × 10 3 ± 4.034954 × 10 2 1.688626 × 10 3 ± 4.222375 × 10 2
f135.887959 × 10 1 ± 3.102544 × 10 1 5.325204 × 10 1 ± 2.649015 × 10 1 3.130917 × 10 1 ± 1.675635 × 10 1
f142.896163 × 10 1 ± 3.059013 × 10 0 2.516679 × 10 1 ± 1.954335 × 10 0 2.528063 × 10 1 ± 2.405140 × 10 0
f153.628762 × 10 1 ± 8.435876 × 10 0 2.536133 × 10 1 ± 3.286879 × 10 0 2.407813 × 10 1 ± 2.946462 × 10 0
f163.678356 × 10 2 ± 8.492899 × 10 1 3.514246 × 10 2 ± 1.238785 × 10 2 4.255240 × 10 2 ± 1.338552 × 10 2
f173.210085 × 10 2 ± 8.169067 × 10 1 2.288290 × 10 2 ± 7.119400 × 10 1 2.673055 × 10 2 ± 7.762843 × 10 1
f183.657321 × 10 1 ± 9.959862 × 10 0 2.552774 × 10 1 ± 2.227834 × 10 0 2.530384 × 10 1 ± 1.768854 × 10 0
f192.229348 × 10 1 ± 5.693372 × 10 0 1.598069 × 10 1 ± 3.193897 × 10 0 1.473478 × 10 1 ± 2.609984 × 10 0
f201.403643 × 10 2 ± 6.340770 × 10 1 1.545319 × 10 2 ± 5.105286 × 10 1 1.195946 × 10 2 ± 5.365881 × 10 1
f212.167895 × 10 2 ± 2.526317 × 10 0 2.112500 × 10 2 ± 2.551442 × 10 0 2.176929 × 10 2 ± 2.782886 × 10 0
f221.090866 × 10 3 ± 1.616684 × 10 3 2.492513 × 10 3 ± 1.798316 × 10 3 2.005652 × 10 3 ± 1.887346 × 10 3
f234.306553 × 10 2 ± 5.383210 × 10 0 4.292150 × 10 2 ± 5.785777 × 10 0 4.311355 × 10 2 ± 6.917497 × 10 0
f245.070785 × 10 2 ± 3.570835 × 10 0 5.065977 × 10 2 ± 3.750386 × 10 0 5.083627 × 10 2 ± 4.562561 × 10 0
f255.403542 × 10 2 ± 2.684350 × 10 1 4.823364 × 10 2 ± 6.109998 × 10 0 4.811605 × 10 2 ± 3.115865 × 10 0
f261.222280 × 10 3 ± 5.892037 × 10 1 1.137749 × 10 3 ± 5.618732 × 10 1 1.127770 × 10 3 ± 5.926815 × 10 1
f275.254951 × 10 2 ± 1.084758 × 10 1 5.098628 × 10 2 ± 9.343850 × 10 0 5.120206 × 10 2 ± 1.296607 × 10 1
f284.901310 × 10 2 ± 2.315774 × 10 1 4.625180 × 10 2 ± 1.261646 × 10 1 4.598068 × 10 2 ± 6.772419 × 10 0
f293.695921 × 10 2 ± 2.035753 × 10 1 3.693345 × 10 2 ± 1.245696 × 10 1 3.660005 × 10 2 ± 1.532970 × 10 1
f306.485235 × 10 5 ± 6.118087 × 10 4 6.093376 × 10 5 ± 3.282382 × 10 4 6.120136 × 10 5 ± 4.048972 × 10 4
FunctionL-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f20.0 ± 0.00.0 ± 0.00.0 ± 0.0
f30.0 ± 0.00.0 ± 0.00.0 ± 0.0
f46.309933 × 10 1 ± 4.861348 × 10 1 5.566785 × 10 1 ± 4.698233 × 10 1 4.758795 × 10 1 ± 4.251127 × 10 1
f51.151320 × 10 1 ± 2.639832 × 10 0 1.326464 × 10 1 ± 2.814235 × 10 0 1.142791 × 10 1 ± 3.021920 × 10 0
f68.654442 × 10 7 ± 9.973005 × 10 7 3.268284 × 10 7 ± 6.168995 × 10 7 1.133130 × 10 6 ± 1.432097 × 10 6
f76.261671 × 10 1 ± 2.090590 × 10 0 6.495951 × 10 1 ± 2.964682 × 10 0 6.340675 × 10 1 ± 2.379764 × 10 0
f81.129984 × 10 1 ± 2.991564 × 10 0 1.405517 × 10 1 ± 2.486846 × 10 0 1.080723 × 10 1 ± 2.867020 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f103.291896 × 10 3 ± 3.546988 × 10 2 3.251439 × 10 3 ± 3.159284 × 10 2 3.280655 × 10 3 ± 3.345557 × 10 2
f112.434579 × 10 1 ± 3.577627 × 10 0 2.366605 × 10 1 ± 3.226247 × 10 0 2.435336 × 10 1 ± 3.915678 × 10 0
f121.221989 × 10 3 ± 3.673815 × 10 2 1.231565 × 10 3 ± 3.631003 × 10 2 1.101802 × 10 3 ± 2.769593 × 10 2
f132.996118 × 10 1 ± 1.994031 × 10 1 2.509417 × 10 1 ± 2.325881 × 10 1 2.884285 × 10 1 ± 2.311451 × 10 1
f142.469884 × 10 1 ± 2.257149 × 10 0 2.434175 × 10 1 ± 2.010046 × 10 0 2.408265 × 10 1 ± 2.364445 × 10 0
f152.124024 × 10 1 ± 1.743464 × 10 0 2.098253 × 10 1 ± 1.794117 × 10 0 2.006068 × 10 1 ± 1.504636 × 10 0
f163.748420 × 10 2 ± 1.328029 × 10 2 4.129801 × 10 2 ± 1.369084 × 10 2 3.800765 × 10 2 ± 1.210071 × 10 2
f172.439444 × 10 2 ± 8.263835 × 10 1 2.795374 × 10 2 ± 1.127910 × 10 2 2.368689 × 10 2 ± 9.354595 × 10 1
f182.261803 × 10 1 ± 1.476114 × 10 0 2.290233 × 10 1 ± 1.433200 × 10 0 2.229708 × 10 1 ± 9.553674 × 10 1
f191.152990 × 10 1 ± 2.462608 × 10 0 1.167486 × 10 1 ± 3.407432 × 10 0 1.085855 × 10 1 ± 2.264789 × 10 0
f201.311696 × 10 2 ± 5.913158 × 10 1 1.463096 × 10 2 ± 8.239168 × 10 1 1.140987 × 10 2 ± 5.954660 × 10 1
f212.134822 × 10 2 ± 3.087330 × 10 0 2.150670 × 10 2 ± 2.682522 × 10 0 2.122163 × 10 2 ± 3.509845 × 10 0
f222.137994 × 10 3 ± 1.867652 × 10 3 1.807509 × 10 3 ± 1.762836 × 10 3 1.971311 × 10 3 ± 1.854554 × 10 3
f234.291962 × 10 2 ± 6.293038 × 10 0 4.299165 × 10 2 ± 5.907164 × 10 0 4.310311 × 10 2 ± 8.161084 × 10 0
f245.082245 × 10 2 ± 4.368923 × 10 0 5.068221 × 10 2 ± 3.243129 × 10 0 5.076325 × 10 2 ± 4.327285 × 10 0
f254.806366 × 10 2 ± 2.299212 × 10 0 4.809218 × 10 2 ± 2.728328 × 10 0 4.804641 × 10 2 ± 1.608251 × 10 0
f261.137797 × 10 3 ± 6.906183 × 10 1 1.140350 × 10 3 ± 5.343276 × 10 1 1.129101 × 10 3 ± 6.948572 × 10 1
f275.091280 × 10 2 ± 1.372703 × 10 1 5.062935 × 10 2 ± 8.969971 × 10 0 5.027209 × 10 2 ± 1.332971 × 10 1
f284.588490 × 10 2 ± 5.684342 × 10 14 4.588490 × 10 2 ± 5.684342 × 10 14 4.588490 × 10 2 ± 5.684342 × 10 14
f293.689658 × 10 2 ± 1.435579 × 10 1 3.645692 × 10 2 ± 1.356443 × 10 1 3.618460 × 10 2 ± 1.293478 × 10 1
f306.042557 × 10 5 ± 3.470244 × 10 4 6.052686 × 10 5 ± 2.884529 × 10 4 6.030053 × 10 5 ± 3.227074 × 10 4
Table 5. Results for CEC 2017 problems, 100 D .
Table 5. Results for CEC 2017 problems, 100 D .
FunctionL-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f22.115272 × 10 1 ± 7.398144 × 10 1 3.621475 × 10 5 ± 2.321462 × 10 6 2.681979 × 10 1 ± 1.323709 × 10 2
f32.318827 × 10 6 ± 3.171098 × 10 6 1.189379 × 10 5 ± 1.466348 × 10 5 1.798921 × 10 7 ± 2.335796 × 10 7
f41.058027 × 10 2 ± 6.168533 × 10 1 1.841503 × 10 2 ± 3.456924 × 10 1 1.985144 × 10 2 ± 1.507302 × 10 1
f56.018822 × 10 1 ± 7.650742 × 10 0 2.881188 × 10 1 ± 5.095198 × 10 0 3.903405 × 10 1 ± 4.302748 × 10 0
f62.768146 × 10 2 ± 2.094528 × 10 2 6.501323 × 10 4 ± 1.049251 × 10 3 1.157600 × 10 4 ± 3.452856 × 10 4
f71.548649 × 10 2 ± 6.128553 × 10 0 1.344337 × 10 2 ± 3.132118 × 10 0 1.395379 × 10 2 ± 5.342167 × 10 0
f85.553710 × 10 1 ± 7.466852 × 10 0 2.763175 × 10 1 ± 4.045294 × 10 0 3.991552 × 10 1 ± 3.695467 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f109.323340 × 10 3 ± 4.989688 × 10 2 9.672592 × 10 3 ± 5.809016 × 10 2 9.829538 × 10 3 ± 4.826719 × 10 2
f114.977651 × 10 2 ± 1.154684 × 10 2 1.398680 × 10 2 ± 4.061085 × 10 1 1.037060 × 10 2 ± 2.911234 × 10 1
f122.859917 × 10 4 ± 1.360122 × 10 4 1.881319 × 10 4 ± 9.745573 × 10 3 1.760383 × 10 4 ± 7.710254 × 10 3
f134.605786 × 10 2 ± 2.511310 × 10 2 1.473008 × 10 2 ± 3.623802 × 10 1 1.498773 × 10 2 ± 3.307054 × 10 1
f142.557099 × 10 2 ± 4.291528 × 10 1 7.044519 × 10 1 ± 1.206246 × 10 1 6.030425 × 10 1 ± 9.448970 × 10 0
f152.477823 × 10 2 ± 5.117302 × 10 1 2.053792 × 10 2 ± 4.548641 × 10 1 1.577945 × 10 2 ± 3.308829 × 10 1
f161.556113 × 10 3 ± 2.800367 × 10 2 1.608918 × 10 3 ± 2.842634 × 10 2 1.718316 × 10 3 ± 3.021229 × 10 2
f171.202785 × 10 3 ± 1.818120 × 10 2 1.096731 × 10 3 ± 2.224261 × 10 2 1.205546 × 10 3 ± 1.764183 × 10 2
f182.505036 × 10 2 ± 5.809274 × 10 1 2.154752 × 10 2 ± 5.280900 × 10 1 1.852876 × 10 2 ± 4.461551 × 10 1
f191.872595 × 10 2 ± 2.511073 × 10 1 1.325648 × 10 2 ± 2.787079 × 10 1 9.679019 × 10 1 ± 2.009190 × 10 1
f201.315492 × 10 3 ± 2.221504 × 10 2 1.415311 × 10 3 ± 1.801826 × 10 2 1.400301 × 10 3 ± 1.908523 × 10 2
f212.683218 × 10 2 ± 5.769608 × 10 0 2.521659 × 10 2 ± 4.632993 × 10 0 2.616113 × 10 2 ± 6.682406 × 10 0
f221.040296 × 10 4 ± 5.472091 × 10 2 1.082725 × 10 4 ± 3.703361 × 10 2 1.069680 × 10 4 ± 5.856186 × 10 2
f235.802729 × 10 2 ± 1.030188 × 10 1 5.690771 × 10 2 ± 1.305997 × 10 1 5.676531 × 10 2 ± 1.120739 × 10 1
f249.211472 × 10 2 ± 1.071240 × 10 1 9.011854 × 10 2 ± 8.832015 × 10 0 9.021128 × 10 2 ± 8.999005 × 10 0
f257.374379 × 10 2 ± 3.838504 × 10 1 7.328269 × 10 2 ± 3.355738 × 10 1 7.234544 × 10 2 ± 4.379540 × 10 1
f263.491742 × 10 3 ± 1.258458 × 10 2 3.205070 × 10 3 ± 9.033496 × 10 1 3.225848 × 10 3 ± 9.340949 × 10 1
f276.178230 × 10 2 ± 2.298992 × 10 1 5.971925 × 10 2 ± 2.012176 × 10 1 5.780971 × 10 2 ± 1.877897 × 10 1
f285.361838 × 10 2 ± 3.127036 × 10 1 5.351704 × 10 2 ± 2.750945 × 10 1 5.269297 × 10 2 ± 2.456277 × 10 1
f291.243196 × 10 3 ± 1.472281 × 10 2 1.155445 × 10 3 ± 1.799635 × 10 2 1.289836 × 10 3 ± 2.137914 × 10 2
f302.440230 × 10 3 ± 1.492281 × 10 2 2.361227 × 10 3 ± 1.419884 × 10 2 2.356999 × 10 3 ± 1.423812 × 10 2
FunctionL-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.05.038721 × 10 7 ± 1.209810 × 10 6
f21.673566 × 10 4 ± 8.143557 × 10 4 2.786236 × 10 3 ± 1.828659 × 10 4 1.619229 × 10 5 ± 9.865543 × 10 5
f32.598700 × 10 6 ± 3.046766 × 10 6 5.200551 × 10 6 ± 5.106623 × 10 6 1.867233 × 10 4 ± 2.098053 × 10 4
f42.001814 × 10 2 ± 1.061642 × 10 1 2.026510 × 10 2 ± 1.067640 × 10 1 1.986945 × 10 2 ± 9.981837 × 10 0
f52.282564 × 10 1 ± 5.374418 × 10 0 2.973385 × 10 1 ± 5.185268 × 10 0 2.126504 × 10 1 ± 4.943471 × 10 0
f64.324963 × 10 5 ± 3.097421 × 10 5 1.037825 × 10 5 ± 1.022042 × 10 5 2.023262 × 10 5 ± 1.675807 × 10 5
f71.258454 × 10 2 ± 4.037642 × 10 0 1.316895 × 10 2 ± 4.837579 × 10 0 1.250752 × 10 2 ± 5.169489 × 10 0
f82.262054 × 10 1 ± 5.062372 × 10 0 2.849150 × 10 1 ± 5.048100 × 10 0 2.191965 × 10 1 ± 4.776998 × 10 0
f90.0 ± 0.00.0 ± 0.00.0 ± 0.0
f101.003350 × 10 4 ± 5.181942 × 10 2 9.919141 × 10 3 ± 5.053399 × 10 2 9.883996 × 10 3 ± 5.542053 × 10 2
f115.681264 × 10 1 ± 2.615972 × 10 1 6.601018 × 10 1 ± 3.161370 × 10 1 6.454649 × 10 1 ± 3.386242 × 10 1
f121.084098 × 10 4 ± 4.835960 × 10 3 1.171083 × 10 4 ± 5.317012 × 10 3 1.629061 × 10 4 ± 7.410286 × 10 3
f131.252693 × 10 2 ± 2.955466 × 10 1 1.221712 × 10 2 ± 3.504948 × 10 1 1.239350 × 10 2 ± 3.319592 × 10 1
f144.026492 × 10 1 ± 4.654467 × 10 0 3.953127 × 10 1 ± 4.084129 × 10 0 3.520826 × 10 1 ± 3.548234 × 10 0
f159.450216 × 10 1 ± 3.081027 × 10 1 9.918180 × 10 1 ± 3.475378 × 10 1 9.316875 × 10 1 ± 3.463347 × 10 1
f161.762196 × 10 3 ± 2.895494 × 10 2 1.653987 × 10 3 ± 3.149891 × 10 2 1.610269 × 10 3 ± 3.171456 × 10 2
f171.177577 × 10 3 ± 2.796055 × 10 2 1.214972 × 10 3 ± 2.455955 × 10 2 1.195693 × 10 3 ± 2.879482 × 10 2
f189.865988 × 10 1 ± 2.318114 × 10 1 9.849790 × 10 1 ± 2.673075 × 10 1 8.753706 × 10 1 ± 1.577905 × 10 1
f195.375456 × 10 1 ± 6.426384 × 10 0 5.332317 × 10 1 ± 6.965122 × 10 0 4.849865 × 10 1 ± 6.528906 × 10 0
f201.361891 × 10 3 ± 2.361198 × 10 2 1.382355 × 10 3 ± 1.833415 × 10 2 1.440578 × 10 3 ± 2.375126 × 10 2
f212.498712 × 10 2 ± 4.793399 × 10 0 2.528137 × 10 2 ± 6.324763 × 10 0 2.474338 × 10 2 ± 4.624437 × 10 0
f221.089671 × 10 4 ± 5.882543 × 10 2 1.056397 × 10 4 ± 5.692429 × 10 2 1.064701 × 10 4 ± 6.777621 × 10 2
f235.685648 × 10 2 ± 1.051014 × 10 1 5.673360 × 10 2 ± 9.183117 × 10 0 5.682635 × 10 2 ± 9.777963 × 10 0
f248.976863 × 10 2 ± 8.162956 × 10 0 8.979166 × 10 2 ± 6.615433 × 10 0 8.970015 × 10 2 ± 7.936307 × 10 0
f256.981912 × 10 2 ± 4.440744 × 10 1 7.117357 × 10 2 ± 4.069772 × 10 1 6.780342 × 10 2 ± 4.988064 × 10 1
f263.120648 × 10 3 ± 8.549825 × 10 1 3.144199 × 10 3 ± 9.002875 × 10 1 3.076108 × 10 3 ± 9.972924 × 10 1
f275.746185 × 10 2 ± 2.402400 × 10 1 5.729895 × 10 2 ± 1.855867 × 10 1 5.600763 × 10 2 ± 2.094522 × 10 1
f285.161411 × 10 2 ± 1.977592 × 10 1 5.203338 × 10 2 ± 2.475601 × 10 1 5.205190 × 10 2 ± 2.297309 × 10 1
f291.268410 × 10 3 ± 2.138615 × 10 2 1.240134 × 10 3 ± 1.761199 × 10 2 1.265436 × 10 3 ± 1.745885 × 10 2
f302.301197 × 10 3 ± 1.302743 × 10 2 2.257047 × 10 3 ± 1.092659 × 10 2 2.254460 × 10 3 ± 1.133409 × 10 2
Table 6. Statistical comparison of algorithms with and without LBR, Mann–Whitney test, CEC 2020 problems.
Table 6. Statistical comparison of algorithms with and without LBR, Mann–Whitney test, CEC 2020 problems.
Algorithm 5 D 10 D 15 D 20 D
L-SHADE vs. L-SHADE-LBR0+/8 = /1-2+/7 = /1-2+/6 = /2-1+/8 = /1-
jSO vs. jSO-LBR0+/9 = /0-1+/9 = /0-3+/7 = /0-2+/8 = /0-
DISH vs. DISH-LBR0+/9 = /0-2+/8 = /0-2+/7 = /1-2+/8 = /0-
Table 7. Results for CEC 2020 problems, 5 D .
Table 7. Results for CEC 2020 problems, 5 D .
Func.L-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f22.10142 × 10 1 ± 9.06803 × 10 2 4.26469 × 10 1 ± 4.06753 × 10 1 3.67067 × 10 1 ± 2.07562 × 10 1
f35.28447 × 10 0 ± 1.11465 × 10 1 5.25455 × 10 0 ± 1.00384 × 10 0 5.26472 × 10 0 ± 2.30645 × 10 1
f41.12273 × 10 1 ± 2.46869 × 10 2 9.77072 × 10 2 ± 4.59113 × 10 2 1.09207 × 10 1 ± 4.67916 × 10 2
f50.0 ± 0.00.0 ± 0.06.24099 × 10 2 ± 1.87230 × 10 1
f60.0 ± 0.00.0 ± 0.00.0 ± 0.0
f80.0 ± 0.00.0 ± 0.00.0 ± 0.0
f91.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f103.42630 × 10 2 ± 1.42101 × 10 1 3.39473 × 10 2 ± 1.76526 × 10 1 3.39473 × 10 2 ± 1.76526 × 10 1
Func.L-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f24.26469 × 10 1 ± 4.06753 × 10 1 1.15746 × 10 0 ± 2.03217 × 10 0 1.35031 × 10 0 ± 2.35366 × 10 0
f35.25455 × 10 0 ± 1.00384 × 10 0 5.36259 × 10 0 ± 1.81877 × 10 1 5.01955 × 10 0 ± 1.19377 × 10 0
f49.77072 × 10 2 ± 4.59113 × 10 2 8.97824 × 10 2 ± 5.08256 × 10 2 9.31253 × 10 2 ± 5.43166 × 10 2
f50.0 ± 0.06.24099 × 10 2 ± 1.87230 × 10 1 0.0 ± 0.0
f60.0 ± 0.00.0 ± 0.00.0 ± 0.0
f80.0 ± 0.00.0 ± 0.00.0 ± 0.0
f91.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f103.39473 × 10 2 ± 1.76526 × 10 1 3.36315 × 10 2 ± 2.00340 × 10 1 3.19473 × 10 2 ± 8.65478 × 10 1
Table 8. Results for CEC 2020 problems, 10 D .
Table 8. Results for CEC 2020 problems, 10 D .
Func.L-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f21.17808 × 10 0 ± 1.45778 × 10 0 1.84593 × 10 0 ± 1.83951 × 10 0 1.84081 × 10 0 ± 1.86378 × 10 0
f31.09229 × 10 1 ± 2.57045 × 10 1 1.14946 × 10 1 ± 3.42269 × 10 1 1.14560 × 10 1 ± 3.15783 × 10 1
f42.66548 × 10 1 ± 3.96601 × 10 2 2.80188 × 10 1 ± 4.31393 × 10 2 2.68712 × 10 1 ± 4.43370 × 10 2
f51.29016 × 10 1 ± 1.31744 × 10 1 1.05751 × 10 1 ± 1.28970 × 10 1 1.53622 × 10 1 ± 1.78070 × 10 1
f67.61459 × 10 2 ± 4.86606 × 10 2 1.70091 × 10 2 ± 1.63884 × 10 2 6.68416 × 10 2 ± 3.85526 × 10 2
f77.99651 × 10 2 ± 1.71343 × 10 1 6.86388 × 10 2 ± 1.52324 × 10 1 1.40205 × 10 1 ± 2.24697 × 10 1
f81.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f92.34825 × 10 2 ± 1.07191 × 10 2 2.21541 × 10 2 ± 1.13696 × 10 2 2.02562 × 10 2 ± 1.02820 × 10 2
f104.16120 × 10 2 ± 2.22173 × 10 1 3.98009 × 10 2 ± 0.03.98009 × 10 2 ± 0.0
Func.L-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f21.84593 × 10 0 ± 1.83951 × 10 0 6.33319 × 10 1 ± 9.86976 × 10 1 1.86907 × 10 0 ± 1.69093 × 10 0
f31.14946 × 10 1 ± 3.42269 × 10 1 1.11115 × 10 1 ± 3.65883 × 10 1 1.13203 × 10 1 ± 3.67866 × 10 1
f42.80188 × 10 1 ± 4.31393 × 10 2 2.71553 × 10 1 ± 4.90826 × 10 2 2.70102 × 10 1 ± 4.86311 × 10 2
f51.05751 × 10 1 ± 1.28970 × 10 1 3.67719 × 10 1 ± 2.19512 × 10 1 4.29479 × 10 2 ± 8.28999 × 10 2
f61.70091 × 10 2 ± 1.63884 × 10 2 3.52478 × 10 1 ± 1.48118 × 10 1 1.31651 × 10 2 ± 1.05516 × 10 2
f76.86388 × 10 2 ± 1.52324 × 10 1 1.98363 × 10 1 ± 2.13177 × 10 1 1.12443 × 10 2 ± 5.56388 × 10 2
f81.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f92.21541 × 10 2 ± 1.13696 × 10 2 2.01274 × 10 2 ± 1.09669 × 10 2 2.21412 × 10 2 ± 1.13577 × 10 2
f103.98009 × 10 2 ± 0.03.98000 × 10 2 ± 4.77485 × 10 2 3.98000 × 10 2 ± 4.77485 × 10 2
Table 9. Results for CEC 2020 problems, 15 D .
Table 9. Results for CEC 2020 problems, 15 D .
Func.L-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f21.11444 × 10 0 ± 1.55539 × 10 0 1.02993 × 10 0 ± 1.03052 × 10 0 3.18849 × 10 0 ± 2.45075 × 10 0
f31.56548 × 10 1 ± 1.19675 × 10 1 1.59237 × 10 1 ± 2.36877 × 10 1 1.59270 × 10 1 ± 2.35040 × 10 1
f43.49454 × 10 1 ± 3.99602 × 10 2 3.62932 × 10 1 ± 5.18875 × 10 2 3.70914 × 10 1 ± 4.04489 × 10 2
f52.72357 × 10 0 ± 1.21873 × 10 0 8.36060 × 10 1 ± 2.85275 × 10 1 2.82653 × 10 0 ± 1.14969 × 10 0
f61.87866 × 10 0 ± 1.91004 × 10 0 7.52602 × 10 1 ± 7.72092 × 10 2 7.95332 × 10 1 ± 3.08505 × 10 1
f75.20809 × 10 1 ± 6.24344 × 10 2 5.13873 × 10 1 ± 5.19186 × 10 2 7.84397 × 10 1 ± 1.13748 × 10 1
f81.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f93.81793 × 10 2 ± 2.42286 × 10 1 3.88236 × 10 2 ± 2.36353 × 10 0 3.89188 × 10 2 ± 1.14562 × 10 0
f104.00000 × 10 2 ± 0.04.00000 × 10 2 ± 0.04.00000 × 10 2 ± 0.0
Func.L-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f21.02993 × 10 0 ± 1.03052 × 10 0 1.28117 × 10 0 ± 1.36090 × 10 0 9.91480 × 10 1 ± 1.03583 × 10 0
f31.59237 × 10 1 ± 2.36877 × 10 1 1.57505 × 10 1 ± 1.53824 × 10 1 1.59450 × 10 1 ± 2.06851 × 10 1
f43.62932 × 10 1 ± 5.18875 × 10 2 3.80445 × 10 1 ± 4.43296 × 10 2 3.72271 × 10 1 ± 4.30883 × 10 2
f58.36060 × 10 1 ± 2.85275 × 10 1 1.52051 × 10 0 ± 9.02712 × 10 1 8.68646 × 10 1 ± 4.30019 × 10 1
f67.52602 × 10 1 ± 7.72092 × 10 2 6.77266 × 10 1 ± 1.30617 × 10 1 7.71839 × 10 1 ± 1.49955 × 10 1
f75.13873 × 10 1 ± 5.19186 × 10 2 8.95471 × 10 1 ± 1.24307 × 10 1 5.06856 × 10 1 ± 3.73802 × 10 2
f81.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f93.88236 × 10 2 ± 2.36353 × 10 0 3.89321 × 10 2 ± 1.06397 × 10 0 3.89540 × 10 2 ± 6.90534 × 10 1
f104.00000 × 10 2 ± 0.04.00000 × 10 2 ± 0.04.00000 × 10 2 ± 0.0
Table 10. Results for CEC 2020 problems, 20 D .
Table 10. Results for CEC 2020 problems, 20 D .
Func.L-SHADEjSODISH
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f21.21924 × 10 1 ± 4.36906 × 10 2 2.92990 × 10 1 ± 4.85360 × 10 1 2.48343 × 10 1 ± 4.09772 × 10 1
f32.04116 × 10 1 ± 5.76301 × 10 2 2.05948 × 10 1 ± 1.81538 × 10 1 2.07393 × 10 1 ± 2.76334 × 10 1
f44.42784 × 10 1 ± 4.24157 × 10 2 4.64267 × 10 1 ± 4.97379 × 10 2 4.63768 × 10 1 ± 4.83692 × 10 2
f51.01366 × 10 2 ± 6.99335 × 10 1 2.01814 × 10 0 ± 1.12021 × 10 0 6.41469 × 10 0 ± 4.26252 × 10 0
f62.31513 × 10 1 ± 5.99600 × 10 2 2.66587 × 10 1 ± 6.28029 × 10 2 9.35968 × 10 1 ± 3.25339 × 10 1
f79.97161 × 10 1 ± 1.46926 × 10 0 6.74912 × 10 1 ± 1.05720 × 10 1 6.67590 × 10 1 ± 2.66941 × 10 1
f81.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f93.97749 × 10 2 ± 1.21733 × 10 0 3.97785 × 10 2 ± 1.10925 × 10 0 3.97699 × 10 2 ± 2.99098 × 10 0
f104.13658 × 10 2 ± 3.04777 × 10 3 4.13657 × 10 2 ± 1.79505 × 10 4 4.13657 × 10 2 ± 1.13687 × 10 13
Func.L-SHADE-LBRjSO-LBRDISH-LBR
f10.0 ± 0.00.0 ± 0.00.0 ± 0.0
f22.92990 × 10 1 ± 4.85360 × 10 1 2.58541 × 10 1 ± 4.19616 × 10 1 2.51160 × 10 1 ± 4.09793 × 10 1
f32.05948 × 10 1 ± 1.81538 × 10 1 2.05692 × 10 1 ± 2.22046 × 10 1 2.05423 × 10 1 ± 1.49670 × 10 1
f44.64267 × 10 1 ± 4.97379 × 10 2 4.64425 × 10 1 ± 3.98595 × 10 2 4.66623 × 10 1 ± 4.65800 × 10 2
f52.01814 × 10 0 ± 1.12021 × 10 0 3.72641 × 10 0 ± 1.37385 × 10 0 2.11216 × 10 0 ± 1.23587 × 10 0
f62.66587 × 10 1 ± 6.28029 × 10 2 3.30280 × 10 1 ± 2.01361 × 10 1 2.39600 × 10 1 ± 5.25111 × 10 2
f76.74912 × 10 1 ± 1.05720 × 10 1 8.89405 × 10 1 ± 1.51139 × 10 1 6.72319 × 10 1 ± 1.19042 × 10 1
f81.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.01.00000 × 10 2 ± 0.0
f93.97785 × 10 2 ± 1.10925 × 10 0 3.97576 × 10 2 ± 1.65870 × 10 0 3.98260 × 10 2 ± 1.24851 × 10 0
f104.13657 × 10 2 ± 1.79505 × 10 4 4.13657 × 10 2 ± 1.79505 × 10 4 4.13657 × 10 2 ± 1.13687 × 10 13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stanovov, V.; Akhmedova, S.; Semenkin, E. Differential Evolution with Linear Bias Reduction in Parameter Adaptation. Algorithms 2020, 13, 283. https://doi.org/10.3390/a13110283

AMA Style

Stanovov V, Akhmedova S, Semenkin E. Differential Evolution with Linear Bias Reduction in Parameter Adaptation. Algorithms. 2020; 13(11):283. https://doi.org/10.3390/a13110283

Chicago/Turabian Style

Stanovov, Vladimir, Shakhnaz Akhmedova, and Eugene Semenkin. 2020. "Differential Evolution with Linear Bias Reduction in Parameter Adaptation" Algorithms 13, no. 11: 283. https://doi.org/10.3390/a13110283

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop