A Novel Tournament Selection Based Differential Evolution Variant for Continuous Optimization Problems

Differential evolution (DE) is a powerful global optimization algorithm which has been studied intensively by many researchers in the recent years. A number of variants have been established for the algorithm that makes DE more applicable. However, most of the variants are suffering from the problems of convergence speed and local optima. A novel tournament based parent selection variant of DE algorithm is proposed in this research.The proposed variant enhances searching capability and improves convergence speed of DE algorithm. This paper also presents a novel statistical comparison of existing DE mutation variants which categorizes these variants in terms of their overall performance. Experimental results show that the proposed DE variant has significance performance over other DE mutation variants.


Introduction
Storn and Price [1] have proposed DE in 1995 which is a stochastic population based evolutionary algorithm to find global optima for any given function [2].DE has been applied to variety of real life problems such as electrical power systems [3], control systems [4], electromagnetism [5], image processing [6], bioinformatics [7], signal processing [8], chemical engineering [9], and many more.Different state-of-the-art versions of DE algorithm like jDE [10], SaDE [11], JADE [12], EPSDE [13], and so forth are based on parameter selection, parameter adaption, strategy selection, and/or strategy adaption mechanisms.DE state-of-the-art proves to be very powerful and uses conventional DE strategies as a strategy selection pool/strategy adaption and/or uses parameter selection pool/parameter adaption.During the process of optimization, DE algorithm evolves population of potential solutions by exploring the entire search space overtime to locate the optima.
The DE algorithm mutation strategies are formed by the linear combination of existing population members.The trial vector and target vector forms the mutant vector in DE.DE mutation strategies can be formed by any combinations of current vector, random vector(s), better vector, and best vector.In any mutation strategy, the order, number, and name of vector(s) are very important.The behavior of DE algorithm is influenced by the selection of mutation strategy and crossover scheme along with their control parameters: mutation probability "" and crossover rate "CR." As stated, there are many mutation strategies associated with DE algorithm but current literature does not show any concrete study where all of the popular mutation strategies are analyzed.It has been observed by the authors that some of the DE mutation strategies are more useful than the current dominating ones.Therefore, in this paper, a comparative analysis is also carried out to identify best performer DE mutation variants.Furthermore, this research also proposes a novel DE mutation variant.The proposed DE variant is based on the well known tournament selection criteria.The proposed DE mutation strategy is abbreviated as TSDE in this paper.To generate a mutant vector in DE algorithm, some parents will be selected through tournament selection criteria.The TSDE variant has the ability to increase the convergence speed since it uses the combination of best vectors and random vectors from the population.Tournament selection will be helpful to incorporate some diversity in TSDE since global best vector has less chance of selection as a parent.To generate a mutant vector, the proposed TSDE utilizes two best individuals instead of global best from the whole population which makes it more convenient to escape from the local optima problem.
The most state-of-the-art DE algorithms focuses on basic DE mutation strategies and the incorporation of this proposed TSDE mutation variant will be helpful in improving the performance of DE algorithm.
The rest of the paper is organized into the following sections.Section 2 presents literature survey.In Section 3, the DE algorithm operators are discussed.In Section 4, the proposed TSDE variant is presented.Section 5 presents test functions and parameter setting details.In Section 6, experimental results along with discussion are given.Section 7 contains statistical analysis.Section 8 compares proposed TSDE with other well known heuristic algorithms.Finally, conclusion and future work are presented in Section 9.

Literature Survey
As stated above, many state-of-the-art versions of the DE algorithms are developed, for example, jDE, SaDE, JADE, and DEGL.Variations in the mutation strategies of DE algorithm have also been the result of intensive research in this area.Qu and Suganthan [29] have introduced the concept of ensemble restricted tournament selection in DE algorithm.They have selected a set of  population members and calculated the Euclidean distance of  population members from the current individual.The member having smallest distance from current individual is selected as a tournament winner and it became the current individual if its fitness is also better.Epitropakis et al. [30] have introduced two new mutation strategies that incorporate spatial information of neighbors to form new population individual.Their new mutation strategies "DE/nrand/1" and "DE/nrand/2" corresponding to two well known mutation strategies "DE/rand/1" and "DE/rand/2" are important part of the current literature.Qiu et al. [31] have introduced the concept of tournament selection criteria in DE algorithm.They have selected all individuals used in their mutation strategy by using the concept of tournament selection by taking a tournament of size 3. Their research work is applied on "DE/rand/1," "DE/rand/2," "DE/current to best/2," and "DE/rand to best/2" DE mutation strategies which proved to be competitive to original DE algorithm and also enhanced jDE algorithm.Neighbourhood based mutation strategy for DE algorithm was introduced by Das et al. [32].This strategy implements the concepts of local model that works on small neighbourhood of each individual while global model mutates each individual using the neighbourhood concept of entire population.The concept of neighbourhood is implemented on a commonly used DE mutation strategy "DE/target-tobest/1/bin." Sutton et al. [33] have introduced the concept of focus search by imposing a selective pressure for DE mutation strategies.The selective pressure is implemented by using the concept of rank based selection of vectors used in generating donor vector in DE mutation.A bias is imposed on selecting a parent using rank based selection criteria that is based on the fitness of population members.Epitropakis et al. [34] have introduced the concept of proximity based mutation operator in DE algorithm.Their new framework utilizes the information of neighbours in evolving the population.Each population member is assigned a selection probability that is inversely proportional to the distance from mutated vector.Proximity based mutation is applied to "DE/rand/1, " "DE/best/1," "DE/rand/2," "DE/best/2," "DE/rand-to-best/1," and "DE/rand-to-best/2" mutation strategies.Research result shows the significance of proximity based mutation operator in DE algorithm.
DE algorithm has a number of conventional mutation strategies where the performances of these strategies may differ from each other.Zaharie [35] has analyzed the influence of crossover variant on the choice of control parameters of the DE algorithm.He has used "DE/rand/1" for different parameter settings of population size, mutation probability, and the crossover rate by binomial and exponential schemes.Babu and Munawar [36] have applied various DE strategies for the optimal design of shell and tube heat problem which performed better than "DE/rand" for the same problem.They have used different DE variants for the selection of optimal design using various control parameter settings.They conclude that "DE/best" strategies have better performance than "DE/rand" for shell and tube heat problem.Jeyakumar and Velayutham [24] present empirical performance analysis of DE variants for unconstrained global optimization problems where they have shown that binomial variants have superior performance compared to the exponential variants.Mezura-Montes et al. [37] have performed comparative study of DE mutation variants for continuous function optimization application.They found through the empirical study that "DE/best/1/bin" has produced good results.Ali et al. [38] have compared Mixed Strategies Differential Evolution (MSDE) with the five different DE variants for function optimization.They have shown that performance of MSDE is better than five variants of DE algorithm.Mutation operators are also used in other well known algorithms like PSO [39].
Qin et al. [11] introduced both control parameter adaption and strategy adaption mechanisms in DE algorithm (SaDE).In strategy adaption scheme, they used a strategy candidate pool of four conventional strategies: "DE/rand/1/bin, " "DE/rand/2/bin," "DE/rand-to-best/2/bin," and "DE/currentto-rand/1." Each target vector generates a trial vector based on the learning period (LP) over previous generations based on the success rate.In parameter adaption, SaDE adjusts control parameter CR based on the median value of CR that is calculated based on previous CR values that have successfully generated trial vector.Zhang and Sanderson [12] introduced new parameter adaption method and have used "DE/current to pbest/1" conventional strategy that is based on "DE/current to best/1."Mallipeddi et al. [13] introduced ensemble based crossover and mutation DE strategy (EPSDE) and their corresponding control parameter scheme.They have used a pool of different crossover and mutation strategies and a pool of values for each associated control parameter.Each target vector generates a trial vector based on the assigned strategy and the parameter values.Successful combination of the mutation and crossover strategy and associated parameters values are stored in the pool.EPSDE uses "DE/current-to-rand/1/bin" and JADE mutation strategy along with binomial and exponential crossover.Islam et al. [40] have introduced modified mutation strategy, Mathematical Problems in Engineering 3

Variant name
Equation a modification in conventional binomial crossover scheme, and a parameter adaption scheme in DE algorithm.They have used "DE/current-to-gr best/1" and call it a less greedy version of "DE/current-to-best/1."Their mutation strategy "DE/current-to-gr best/1" uses best vector from % population of individuals to generate the trial vector for each target vector.In pbest crossover mutant vector can swap p-top ranked individuals of current generation instead of current parent using binomial crossover scheme.In parameter adaption scheme scale factor adaption is based on the Cauchy distribution and crossover probability adaption is based on the Gaussian distribution.Wang et al. [41] have introduced composite DE (CoDE) variant.In this scheme, they have used a pool of three trial vector generation strategies and a pool of three parameter setting combinations.The trial vector strategies used are "rand/1/bin," "rand/2/bin," and "currentto-rand/1" while the parameter setting combinations are [ = 1.0,Cr = 0.1], [ = 1.0,Cr = 0.9], and [ = 0.  1 shows different binomial and exponential mutation strategies used with algorithm.Current literature shows variation in naming in these strategies, which is already indicated in Table 1; for example, Four different names are used in the literature for mutation strategy  7 .Throughout this paper,   denotes the target vector (or current vector),   represents the trial vector, and V  as a mutant vector.In this paper,     states the   th random vector for th generation, V   will refer to th component of donor vector at th generation,  best  states the best vector at th generation, and    refers to current vector at th generation.

DE Algorithm
DE algorithm has three different parameters: a population of size NP, crossover control parameter CR, and difference vector amplification parameter .Each population member in DE is represented as a -dimensional parameter vector.In DE algorithm, population is initialized randomly which is supposed to cover the entire search space.Each vector in the DE is represented by  , , where  = 1, 2, 3, . . ., NP and  is generation number.New offspring in DE algorithm are generated by mutation, crossover, and selection operators.Three different terminologies of vectors donor vector, trial vector, and target vector are used in DE algorithm.Donor vector is a vector that is created in the mutation operation, trial vector is created in the crossover operation, and target vector is the current vector of population.The detail of these operators with index  = 1, 2, 3, . . ., NP,  = 1, 2, 3, . . .,  is as follows.
Mutation.In mutation operation, mutant vector also called donor vector is created.Donor vector V ,+1 of th population member is calculated by adding the weighted difference of two vectors to third vector: where random indices =  3 , and  is mutation probability parameter.
Crossover.DE crossover strategies control the number of inherited components from the mutant vector to form a target vector.Binomial and exponential are main crossover schemes.The DE crossover rate parameter (CR) influences the size of perturbation of the base (target) vector to ensure the population diversity [43].Following are the binomial and exponential crossover schemes.
Selection.In DE algorithm, new population members are formed using selection operation.Selection operator uses the greedy approach by comparing fitness of trial vector  ,+1 with the fitness of target vector  , ; the vector having best fitness is selected as a member of new population.For selection, the following equation is used: where fitness ( ) function calculates the fitness value of objective function.
There are several DE algorithm mutation variants/strategies that are formed by the linear combination of existing population members.The trial vector and target vector form the mutant vector in DE.Throughout this paper,   denotes the target vector (or current vector),  is the running index,   represents the trial vector, and V  is a mutant vector.In DE algorithm, different mutation schemes are used to create the trial vector by using any combination of current, best, and random vectors.The behavior of DE algorithm is influenced by the selection of mutation strategy and crossover scheme along with their control parameters: mutation probability "" and crossover rate "CR" [44].
Difference vector (DV) is the difference of two mutating vectors that is used to form offspring in the population [45].To form the mutant vector in DE, some researcher uses a random value  ∈ (0, 1) as a coefficient multiplier with the first difference vector [1] and mutation probability "" as a coefficient multiplier with the other difference vector(s) [17].Some researchers have used only "" as a coefficient multiplier with the difference vector(s) to form the mutant vector [44].To reduce the number of control parameters of DE algorithm, we use  =  [28].A list of binomial and exponential mutation strategies that are used in this research is given in Table 1.

Proposed Tournament Selection Based DE (TSDE) Variant
Tournament selection is one of the famous selection approaches used in genetic algorithms [46].The parameter associated with this scheme is the size of tournament that selects number of individuals to participate in the competition of tournament.To decrease the risk of premature convergence, the loss of diversity should be kept as low as possible in the population [47].To enhance the performance of population based algorithms, various researches have used tournament selection in their research work [48,49].
The DE algorithm sometimes faces the problem of slow and/or premature convergence [32].Exploration and exploitation are very important aspects that are helpful in improving the convergence acceleration and solution quality of evolutionary algorithms [50].The mutation strategies "DE/rand/1" and "DE/rand/2" are helpful for exploration [11] and "DE/best/1," "DE/rand-to-best," and "DE/best/2" are more helpful for exploitation than exploration due to less population diversity [14].To balance the exploration and exploitation ability of proposed TSDE variant, both random and best vectors are used.
A novel DE mutation strategy based on the selection of parent is introduced in this research work.The new mutation strategy utilizes the knowledge of some best performing vectors and some random vectors in creating the mutant vector.Like mutation strategies "DE/rand/1" and "DE/rand/2," the proposed mutation strategy uses random vector as a base vector along with two difference vectors.To generate a mutant vector in DE, each difference vector utilizes the knowledge of one distinct best vector that will be helpful in incorporating diversity and improving convergence in the proposed mutation strategy.The two difference vectors and two distinct best vectors selected through tournament selection criteria will be helpful in balancing the exploration and exploitation.
To minimize the loss diversity, a tournament of small size is selected so that less fit individuals may have a chance of selection as a parent.TSDE generates two best vectors using tournament selection mechanism.The selection of parents is based on the tournament selection criteria by taking a tournament of size 3 of randomly selected population members.The equation of the proposed variant TSDE is as follows: The

𝑔
is the second best individual from the second tournament.The vectors used in each tournament are selected randomly from the current population.
The proposed TSDE variant has the ability to increase the convergence speed since it uses the combination of best vectors and random vectors from the population.Tournament selection may incorporate some diversity in TSDE since global best vector has less chance of selection as a parent because tournament of small size is used instead of whole population.Selection of two best individuals instead of global best from the whole population makes it more opportune to escape from the local optima problem that accordingly increases the solution quality and convergence speed of DE algorithm.The most state-of-the-art DE algorithm focuses on basic DE mutation strategies and the incorporation of this proposed mutation strategy will be helpful in improving the performance of DE algorithm.
The proposed mutation strategy has the following features.
(1) A random number will be used as a base vector since mutation strategies "DE/rand/1" and "DE/rand/2" use random base vectors that maintain better exploration ability but have slow convergence speed [11].
(2) Although DE mutation strategies "DE/best/1," "DE/ rand-to-best," and "DE/best/2" have fast convergence [14], these are biased towards global best value that may have a chance to stick in local optima problem [11].So to improve the convergence speed of proposed Vector Difference vector Figure 1: Graphical view of proposed TSDE.
TSDE, two best vectors (best 1 , best 2 ) are selected using tournament selection criteria.(3) The equation of proposed mutation strategy is where  is scaling factor of each difference vector.The motivation of this mutation strategy is an existing DE variant "DE/rand to best/1" [23,51] that places perturbation of random vector towards global best vector and contains two difference vectors with equation Two best vectors (best 1 , best 2 ) in TSDE will be helpful in escaping local optima problem and will be helpful in improving the solution quality.(4) Two difference vectors are used where each difference vector utilizes one best vector selected through tournament criteria, since mutation strategies having two difference vectors (DV) produce better perturbation mode than one difference vector mutation strategies [41].
The members of tournament are selected randomly which almost makes it equal chance of selection of higher or lower performing individuals from the population.This variant maintains randomness since members of the tournament are selected randomly and we select the best performing amongst the tournament to ignore the poor performing individuals.The proposed variation is simple since it uses two best individuals selected using the tournament selection technique and other random vector selected randomly from the current population.The performance of TSDE is compared with other DE variants in Section 6.
A graphical representation of proposed mutation strategy is shown in Figure 1, where mutant vector V   is generated Algorithm 1 shows the pseudocode of the proposed DE algorithm (TSDE).The proposed method is also implemented through computer simulation which will be discussed in a later part in the paper.

Test Functions and Parameter Setting
In order to evaluative the performance of the proposed variant and existing variants of DE algorithm, a comprehensive set of 37 -dimensional benchmark functions is used.These benchmark functions are commonly used for multidimensional global optimization problems.These functions are given in Table 2 along with the equation, search space, dimension, and optimum value of each function.-dimensional functions have been used for extensive comparison of DE algorithm variants for various dimensions.
The control parameters of the DE algorithm are Np (population size),  (dimension),  (scaling factor), and CR (crossover control parameter).Different researchers have used different  and CR parameter values for function optimization of -dimensional functions.Like Noman and Iba [52] have used  = 0.9, CR = 0.9, Brest et al. [10] have used  = 0.5, CR = 0.9; Ali and Torn [53] have used CR = 0.5,  ∈ [0.4,1]; Piotrowski and Napiorkowski [54] have used CR = 0.5 and  = 0.5; Dong et al. [55] have used  = 0.7 and CR = 0.5 in their research work.In this research work, we have used control parameters  = 0.7, CR = 0.5 and population size Np = 30.The dimensions size used is 10, 20, and 30 and the numbers of training iterations for average fitness value are 5000, 10000, and 15000 for 10D, 20D, and 30D, respectively.The DE algorithm iterates over 30 independent runs and the average fitness value of these 30 runs is used in the results.Same parameter values are used for all functions and all mutation strategies.Experimental results of number of function calls (NFC) are generated for maximum NFC 10 4 * DIM [56].To find out the number of function calls, VTR is set to 0.0001 and Max-NFC for 10D, 20D, and 30D is, respectively, 100,000, 200,000, and 300,000 for all functions and all mutation strategies.

Experimental Results and Discussion
This section presents simulation results and its comparative analysis for DE mutation strategies given in Table 1 using Step function (separable, unimodal) where Function "15" (separable, multimodal) Mishra 2 global optimization problem (nonseparable, multimodal) MultiModal global optimization problem (separable, multimodal) Quintic global optimization problem (separable, multimodal) Stochastic global optimization problem (separable, multimodal)  parameter setting given in Section 5.In order to perform a fair comparison, all results are averaged over 30 independent runs.Both binomial and exponential crossover schemes are implemented for all benchmark functions given in Table 2 and DE mutation strategies given in Table 1 to obtain the results on a set of benchmark functions given in Table 2 for the DE mutation strategies given in Table 1.Both average fitness value and number of function calls (NFC) are used for performance evaluation of DE mutation strategies as used by most of the evolutionary computing algorithms [14].First, average fitness values of the functions are considered over 30 independent runs for the performances evaluation of DE mutation strategies.For easy observation, best values of mutation strategies are given as bold face for every case.Experimental results of average fitness for all DE mutation strategies and benchmark functions are reported in Table 2.
Because of space limitation, the results of binomial and exponential strategies of proposed TSDE, five most commonly used strategies [14], and one other best performing strategy are reported in Tables 3 and 4. Figures 2-9 show average fitness and NFC Dendrograms for the proposed DE mutation variant and existing DE mutation strategies.Figure 10 shows the convergence graph of binomial and exponential versions of the proposed TSDE and other selected mutation strategies.NFC experimental results for DE mutation strategies and benchmark functions are reported in

Statistical Analysis
The It is important to discuss that a number of DE mutation strategies are considered in this research work and all these strategies perform different from each other.Dendrograms are also considered to be an effective tool to determine the clusters number.For DE strategies, more number of clusters in the Dendrogram indicates that DE mutation strategies have varying performance from one another, which shows the diverse nature of these strategies in DE algorithm.The diverse nature of DE variants will be helpful to make this algorithm more broadly applicable in various optimization problems of distinct characteristics.
The convergence graphs of selected benchmark functions for 10 dimensions (10D) are given in Figure 10.Logarithmic convergence graphs of average fitness show number of iterations horizontally and performance vertically.Horizontal iterations are shown as a multiple of 100. Figure 10 contains the convergence graphs of five commonly used DE mutation strategies ( 1 ,  2 ,  3 ,  4 , and  7 ), the proposed mutation variants ( 41 and  42 ), and one other better performing variant  10 .Convergence graphs of few selected functions are given in the paper.Convergence graph shows that the proposed "TSDE( 41 )" mutation variant has better performance in most of the cases and proposed that "TSDE( 42 )" has better performance in few cases.The mutation strategy  10 ("DE/rand to best/1/bin") proves itself to be another better performing DE strategy; it is following "TSDE ( 41 )" in most cases of convergence graphs.It is important to mention here that  10 ("DE/rand to best/1/bin") strategy has never been declared to be one of the best performing variants or among top performing variants of DE algorithm.The convergence graph clearly shows the dominating performance of the proposed DE mutation variant "TSDE( 41 )." Most of the researchers have concentrated on very few basic DE strategies like "DE/rand/1/bin ( 1 )," "DE/best/1/bin ( 2 ), " "DE/rand/2/bin ( 3 )," "DE/best/2/bin ( 4 )," and so forth, and some key performing mutation strategies remained unnoticed in the DE research.This study will prove to be fruitful for DE researchers having choice of different variants for scientific and real time applications.

Comparison of TSDE with PSO and GA
PSO and GA are two well known metaheuristics that are applied on variety of real world problems [57].The simulation results of proposed "TSDE( 41 )" is compared with two other well known metaheuristics: particle swarm optimization and genetic algorithm.In order to compare TSDE with PSO and GA, the parameter setting for GA is crossover rate CR = 0.5 and mutation rate  = 0.05 [58] while for PSO inertia weight  ∈ (0.4, 0.7), 1 = 1.49618, and 2 = 1.49618 are used [59].For fair comparison between PSO, GA, and TSDE, the population size, dimensions, number of iterations, and standard benchmark functions are kept the same.Results of average fitness value and NFC are presented in Tables 5  and 6, respectively.Proposed strategy TSDE has dominating performance than evolutionary metaheuristics PSO and GA in most of the cases of 10D, 20D, and 30D for average fitness value and NFC.TSDE proves itself to be one of the powerful variants of DE algorithm that can be a significant addition in DE research.

Conclusion
In this paper, a novel technique TSDE mutation variant based on the selection of parent vectors to generate the new population is proposed.The proposed variants "TSDE( 41 )" and "TSDE( 42 )" are compared with the existing DE mutation strategies.TSDE( 41 ) mutation variant has dominating performance among the DE mutation strategies.This new mutation variant TSDE will prove to be a valuable addition in DE literature.A comprehensive set of well known dimensional benchmark functions have been used to evaluative the performance of the proposed TSDE as well as existing mutation strategies of DE algorithm.In literature, such an extensive comparison performed with many number of mutation strategies is lacking.This comparison unveils some very beneficial facts related to the DE research improvements.The other important aspect of this research is the statistical comparison of the DE mutation strategies; such a comparison is also not available in the literature.Experimental results consist of average fitness value and number of function call parameters of DE algorithm that are reported in Tables 3 and  4. Since total strategies are 42 and it is not possible to report all the mutation strategies in few tables, so proposed TSDE mutation variant, five commonly used mutation strategies, and one better performing mutation strategy are reported in the results.It can be concluded that there is a deviation in the performance of DE algorithm variants and the selection of DE variant affects the performance result of DE algorithm.Better results can be obtained by choosing the better performing strategies.The proposed binomial "TSDE( 41 )" has the leading performance among the DE mutation strategies for both fitness value and number of function call performance parameters, the proposed exponential "TSDE( 42 )," mutation strategy "DE/rand to best/1/bin ( 10 )," " 7 " and "DE/best/1/bin ( 2 )" are the other better performance mutation strategies.This research work also reveals one of the best performing mutation strategies "DE/rand to best/1/bin ( 10 )" that has been rarely brought into play.This research work will prove to be a significant addition in DE literature.The results of TSDE are also compared with two well known heuristics, GA and PSO, that show dominating performance of TSDE.The further studies are still required to explore the limitations, advantages, and flaws in this direction.The main focus in this work is not the overwhelming and existing DE variants but new tracks to work on the DE variants in optimization.One possible direction of future work of this paper is to develop a mathematical model for the DE mutation strategies.

Table 1 :
List of DE mutation strategies.
equation of this variant uses two best vectors  best 1

Table 4 .
The results of average fitness value and NFC are used to evaluate the performance of DE mutation strategies.Commonly used mutation strategies, proposed mutation  2 , and so forth, convergence graphs show that  41 has better convergence speed than these mutation strategies.Research result shows that proposed mutation strategies perform better for the functions having various characteristics like unimodal/multimodal and separable/nonseparable.The proposed mutation strategy ( 41 ) has better average fitness performance for separable functions  1 ,  2 ,  12 ,  13 ,  19 (10D, 20D),  16 (10D),  22 (10D, 20D),  18 ,  19 ,  20 ,  21 ,  22 (10D),  26 , and  28 ; nonseparable functions  3 ,  7 ,  10 ,  11 ,  18 ,  24 , and  25 ; unimodal functions  2 ,  3 ,  11 ,  12 ,  13 ,  16 (10D),  26 , and  21 ; and multimodal functions  1 ,  7 ,  10 ,  18 ,  19 (10D, 20D),  22 (10D, 20D),  18 ,  19 ,  22 (10D),  24 ,  25 ,  26 , and  28 .One of the rarely used DE mutation strategies  mutation strategy ( 42 ) has better average fitness performance for separable functions  5 (10D, 30D),  12 ,  14 ,  22 ,  22 (10D), and  26 ; nonseparable functions  6 (30D),  24 ,  25 , and  30 ; unimodal functions  12 and  14 ; and multimodal functions  5 (10D, 30D),  6 (30D),  22 ,  22 (10D),  24 ,  25 ,  26 , and  30 .Average fitness result shows that  41 has dominating performance among all strategies and although it shares its best performance with other mutation strategies like  2 ,  7 , and  10 in most of the cases, the convergence graphs show that the proposed TSDE has fast convergence than other mutation strategies like  2 ,  7 , and  10 and vice versa.Results of number of function calls show that TSDE( 41 ) has clearly leading performance among all other DE mutation strategies in most of the cases.The binomial TSDE mutation variant  41 has best performance in most of the functions having various characteristics like separable, nonseparable, unimodal, and multimodal. 41 outperforms other stratefunctions.The mutation strategies that fail to reach VTR for any function with specified parameter contain dashed values in Table 4 for the NFC results.It can be noted from the results that the proposed TSDE performs better than all others including  10 .Proposed TSDE can be one of the powerful mutation strategies in DE algorithm.Experimental result can be summarized that the proposed binomial "TSDE( 41 )" has the dominating performance among all DE mutation strategies for both average fitness values and number of function call values.The strategy 41and  42 ) but reported for selected mutation strategies only due to space issue.Experimental results are analyzed on the basis of commonly used performance parameters average fitness and NFC.The benchmark functions are -dimensional functions having varying nature.From average fitness results, it can be observed that proposed binomial TSDE( 41 ) version has dominating performance among all strategies; the mutation strategy  10 that has never announced to be one of the best performing mutation strategies in DE has the second best performance; proposed exponential TSDE( 42 ) and two other well known DE mutation strategies  7 and  2 have subsequent best performance.Although in most of the cases proposed  41 shares its best performance with other mutation strategies like  10 ,  7 , 10 has second best average fitness performance for separable functions  1 ,  2 ,  12 (10D),  14 ,  19 ,  16 (10D),  18 ,  24 ,  20 ,  21 , and  26 ; nonseparable functions  7 ,  18 , and  21 (10D, 30D); unimodal functions  2 ,  12 (10D),  14 ,  16 (10D),  20 , and  21 ; and multimodal functions  1 ,  7 ,  18 ,  19 ,  21 (10D, 30D),  18 ,  24 , and  26 .The proposed exponential 10 ("DE/rand repeated to best/1/bin") is the second bestFigure 10: 10D: logarithmic convergence graphs of average fitness for some functions showing number of iterations horizontally and average fitness vertically.

Table 4 :
Number of function calls (mean ±

Table 4 :
Continued. in average fitness value and is runner-up in most of the cases of number of function call values.One interesting aspect of this variant ( 10 ) is that the performance of this strategy is never announced to be one the best strategies of the DE algorithm.The mutation strategies  7 and  2 ("DE/best/1/bin") that are already known to be best performing DE variant have subsequent best performance along with exponential "TSDE( 42 )" and vice versa.It can also be summarized that the worst performing mutation strategies among the considered strategies are  6 ("DE/rand repeating & current to rand/1/bin"),  8 ,  11 (DE/rand to best/1),  14 ("DE/rand to best/2")  16 ("DE/rand to current/2"),  19 ("DE/rand/3"), and  20 ("DE/best/3") and vice versa.
results of DE mutation strategies are analyzed statistically by calculating the rank correlation and by generated Dendrograms.To find the rank correlation, rank of each DE mutation strategy is calculated to find the correlation between DE strategies.The performance of mutation strategies will be observed during the analysis for various dimensions with the help of Dendrogram.Similarity/dissimilarity matrix is used to generate the Dendrogram.Similarity matrixes are obtained using "Statistical Package for the Social Sciences" (SPSS).After similarity matrix, pairwise distances between DE strategies are calculated.Pairwise distance is used to create hierarchical cluster tree linkage that is used to create the Dendrogram.Dendrogram actually generates the clusters from the dissimilarity matrix using agglomerative clustering.Dendrograms are considered to be an effective way of cluster analysis that specifies the level of similarity joining of any two clusters.The leaf of clusters in Dendrograms contains the DE mutation strategy's number and -axis contains the dissimilarity scale that shows the joining distance between clusters.The height of each branch in Dendrograms shows that which DE strategies are similar or different from each other.The greater height of DE strategies branches means that strategies are more different from each other.Consider the Dendrograms generated from average fitness values and number of function call values given in Figures2-9, respectively.A similarity measure among the DE strategies is explored by generating Dendrograms from the average fitness values and number of function call performance parameters.The best performing mutation strategies are contained in the similar performing cluster shown in transparent circles and most worst performing strategies are contained in the dotted circle in each Dendrogram.The Dendrograms are generated for 10D, 20D, 30D and the average of all (10D, 20D, and 30D) is generated for both performance metrics of fitness value and number of function calls.The best performing and worst performing strategies among all mutation strategies are contained in the corresponding groups in most of the cases.The best performing strategies "TSDE( 41 )," "DE/rand repeated to best/1/bin ( 10 )," "DE/best/1/bin ( 2 )," and "( 7 )" are contained in the same group in most of the cases that can be considered as best performing mutation strategies that are group shown in the transparent circle in Figures2-9.The worst performing group of mutation strategies is encircled dashed line and it contains "DE/rand repeating & current to rand/1/bin ( 6 )," "DE/rand & current to best/2/bin ( 14 )," "DE/rand & current to rand/1/bin ( 16 )," "DE/rand/3/bin ( 19 )," and "DE/best/3 ( 20 )" and vice versa.

Table 6 :
Number of function calls (mean ±