Cloud Particles Evolution Algorithm

Many evolutionary algorithms have been paid attention to by the researchers and have been applied to solve optimization problems. This paper presents a new optimization method called cloud particles evolution algorithm (CPEA) to solve optimization problems based on cloud formation process and phase transformation of natural substance. The cloud is assumed to have three states in the proposed algorithm. Gaseous state represents the global exploration. Liquid state represents the intermediate process from the global exploration to the local exploitation. Solid state represents the local exploitation. The cloud is composed of descript and independent particles in this algorithm.The cloud particles use phase transformation of three states to realize the global exploration and the local exploitation in the optimization process. Moreover, the cloud particles not only realize the survival of the fittest through competition mechanism but also ensure the diversity of the cloud particles by reciprocity mechanism.The effectiveness of the algorithm is validated upon different benchmark problems. The proposed algorithm is compared with a number of other wellknown optimization algorithms, and the experimental results show that cloud particles evolution algorithm has a higher efficiency than some other algorithms.


Introduction
Many real-world problems which are classified as global optimization problems are very complex and are quite difficult to solve.Various optimization algorithms are developed to solve these problems.Most of these algorithms learn and imitate a variety of intelligent behavior in nature.These algorithms exploit a set of potential solutions and detect the optimal solution through cooperation and competition among the individuals of the population [1].For example, Genetic algorithm [2] proposed by Professor Holland is a computational model of biological evolution which imitates natural selection and genetic mechanism.Ant colony algorithm [3] is a heuristic bionic optimization algorithm which imitates the foraging behavior of the ant colony.In the ant colony algorithm, an individual is called an artificial ant.The artificial ants search probabilistically in the solution space to create candidate solutions.The evaluating and the updating of these candidate solutions are based on the pheromone associated with each one of them.The candidate solution with the maximum amount of pheromone is considered to be the optimal solution of the problem.Particle swarm optimization (PSO) [4] is a bionic intelligent computing method proposed by Kennedy and Eberhart which imitates the flying and the foraging behavior of birds.In PSO, the optimal solution is called a leader.The leader guides each particle to search.Besides, the particles cooperate with each other to realize exploration.These behaviors are directed by a combination of swarm's previous best (gBest) and their own previous best (pBest) [5].The artificial bee colony algorithm (ABC) [6,7] is a nonderivative optimization method based on self-organizational models of bee and swarm intelligence.In ABC, bees are divided into three groups which are called "employed" bees, "onlooker" bees, and "scouts" bees."Onlookers" bees look for food sources.Food source is the search goal of bees."Employed" bees work on food source.In addition, "scouts" bees search for the new food source near the hive.Artificial immune system (AIS) [8] is an adaptive system inspired by the immunology which imitates human immunological functions, principles, and models to solve complex problems.In AIS, the external organisms are called antigen and the internal organisms are called antibodies.Immune clonal operation provides multistrategy conditions for recombination and mutation [9].The clonal selection selects the excellent antibodies from the subpopulation generated by clonal proliferation.

Cloud Particles Evolution Algorithm
2.1.Cloud Formation.The idea of the proposed CPEA is inspired from cloud formation.The state of the cloud particles' changing process [10][11][12] is shown in Figure 1.
After the wet air rises, it will go through a set of changes because of some certain reasons temperature and pressure.Then, the water vapor is formed.Water vapor plays an important role in this process.When the water vapor is saturated, it will adsorb on the cloud condensation nuclei (CNN) and the initial cloud particles are formed.Cloud particles are condensed and sublimated continuously by absorbing water vapor around the cloud.In the process of the condensation, the cloud particles are getting closer and colliding.Then the larger cloud particles are formed.Ice phase is generated when the temperature is below 0 ∘ C at this time.The emergence of ice crystals damages the stable status of cloud phase structure.Ice crystals grow up quickly to be the snow crystals because of the condensation of the water vapor.On the one hand, water vapor transfers to ice crystals and ice crystals begin to grow up as supercooled water droplets evaporate.On the other hand, the supercooled water droplets will be frozen after the colliding to make the snow particles grow bigger.If the supercooled water vapor which has been collided is too many, the snow crystals will be transferred to be solid substances.Then, the solid substances such as snow and hail will be dropped.If the solid substances descent in the area which has a higher temperature over 0 ∘ C, they will melt into raindrop.The water in our planet will be irradiated by the sunshine, and it will become the water vapor again in the air.The cloud is formed and changed like this.

Phase and Phase Transformation.
There are various kinds of substance in nature.Substances commonly exist in three states: solid, liquid, and gaseous.Figure 2 is a simple diagram to show how water transforms among the three states.Evaporation, condensation, and melting are the three ways for water transforming.In order to describe the different forms of the substance, the phase can be used as a symbol for states of substances.In thermodynamics, phase is defined as some parts which have the same quality [13].A phase transformation is the transition which is from one phase to another in thermodynamics.
Phase transformation is relevant to Gibbs free energy.In thermodynamics, the Gibbs free energy [14] is a thermodynamic potential that measures the "useful" or processinitiating work obtainable from a thermodynamic system at a constant temperature and pressure.

Population Evolution Mechanism.
Darwin pointed out that mutation and natural selection are the two important factors for the evolution, which explained many natural phenomena successfully [15].Nowak from Harvard University ranked the cooperation as the third important factor of the species evolution.Evolution is a unity of the cooperation and the opposition [16].Every gene, cell, and organism of the individual should strengthen its own evolution process at the expense of beating its rivals.Therefore, individuals often conflict with each other because of limited resources [17].However, the complex life organizational structure in different levels is produced in the process of evolution, which indicates that the cooperation exists in the evolution.Therefore, Nowak summarizes five mechanisms for the cooperation evolution such as kin selection, direct reciprocity, indirect reciprocity, network reciprocity, and group selection (shown in Figure 3).Haldane believes that the altruism exists in the evolution [18].According to the study of the social insect, Hamilton [19] used mathematic model for quantitative description and defined the Hamilton rules.The rule states natural selection favors cooperation if the donor and the recipient have kin relationships.Direct reciprocity requires several repeated encounters between two individuals which can offer help with each other.The direct reciprocity theory framework is repeated Prisoner's Dilemma.The classical game strategy is tit for tat, TFT [20].The interactions among individuals in nature are often asymmetrical and fleeting.Therefore, indirect reciprocity is a more prevalent reciprocity form.Indirect reciprocity refers to establishing a good reputation by helping others, which will be rewarded by others later.In Indirect reciprocity, one individual acts as donor while the other acts as recipient.Indirect reciprocity will favor the cooperation evolution, if ; the probability of one individual's reputation exceeds the cost-to-benefit ratio of the altruistic act, namely,  > / [21].Network reciprocity is a new form of reciprocity.Cooperators can form network clusters in order to help each other.The cooperation mechanism of this spatial reciprocity is called network reciprocity.Group selection is a minimalist stochastic model [22].The population is divided into several groups.Cooperators help each other in their own groups while defectors do not help.The pure cooperator groups grow faster than the pure defector groups while the defector groups reproduce faster than the cooperator groups in any mixed group.Therefore, defectors are often selected in the lower level (within groups) while cooperators are often selected in the higher (between groups) level.The proposed algorithm which is similar to other metaheuristic algorithms begins with an initial population called the cloud particles.At the beginning, let the initial state be cloud gaseous phase which is composed of many cloud particles (shown in Figure 4).The best individual (best cloud particle) is chosen as a seed or nucleus.Then, the cloud particles are going to have a condensational growth or collisions and coalescence growth (shown in Figure 5).If the population evolves successfully, the cloud particles are going to change from cloud gaseous state to liquid state by condensation.Cloud changes from gaseous phase to cloud liquid phase.This process is called phase transformation.When cloud particles are in liquid state, cloud particles transform from global exploration to local exploitation gradually.In this process, according to the idea of the indirect reciprocity mentioned before, if  > / the algorithm will favor the cooperation evolution or it will change from the liquid phase to solid phase by phase transformation.When cloud particles are in solid state (shown in Figure 6), the cloud particles condense into ice crystals at the beginning, and then ice crystals curdle into hailstones.The hailstone is the best solution.
It can be seen from the biological evolutionary behavior that the individuals compete to survive.However, competitive ability is always asymmetry because the different individual has different fitness.Usually, it is easier for disadvantaged individuals to compete with the individuals in the same situation while they seldom compete with the advantaged individuals [23].Therefore, competition and cooperation which are two different survival strategies for living organisms are formed in order to adapt to the environment.In the CPEA model, according to the natural phenomena that the population develops cooperatively when they evolve, these two strategies can be shown as follows.The populations carry out the competitive evolution in cloud gaseous phase.In cloud liquid phase, the populations employ competitive mechanism for improving local exploitation ability and employ reciprocal mechanism to ensure the diversity of the population.Cloud solid phase indicates the algorithm has found the optimal solution area; therefore the algorithm realizes fast convergence by competitive evolution.The optimization process of CPEA model is shown in Figure 7.
In CPEA, the initialization, phase transformation driving force, condensation operation, reciprocity operation, and solidification operation are described as follows.

Initialization.
When the metaheuristic methods are used to solve the optimization problem, one population is formed at first.Then, the variables which are involved in the optimization problem can be represented as an array.In this algorithm, one single solution is called one "cloud particle." The cloud particles have three parameters [24].They are Expectation (  ), Entropy (  ), and Hyperentropy (  ).  is the expectation values of the distribution for all cloud particles in the domain.  is the range of domain which can be accepted by linguistic values (qualitative concept).In another word,   is ambiguity.  is the dispersion degree of entropy (  ).That is to say,   is the entropy of entropy.In addition,   can be also defined as the uncertainty measure for entropy.In  dimensional optimization problems, a cloud particle is an array of 1 × .This array is defined as follows: To start the optimization algorithm, a matrix of cloud particles which size is × is generated by cloud generator (i.e., population of cloud particles).Cloud particles generation is described as follows: Φ(, ) is the normal random variable which has an expectation  and a variance ;  is the number of cloud particle;  is the cloud particle generated by cloud generator.Consider The matrix  generated by cloud generator is given as (rows and column are the number of population and the number of design variables, resp.) Each of the decision variable values ( 1 ,  2 , . . .,   ) can represent real values for continuous or discrete problems.The fitness of a cloud particle is calculated by the evaluation of objective function () given as In the function above,  and  are the number of cloud particles and the number of decision variables.The cloud particles which number is  are generated by the cloud generator.The cloud particle which has the optimal value is selected as the nucleus.The total population consists of   subpopulations.Each subpopulation of cloud particles selects a cloud particle as the nucleus.Therefore, there are   nuclei.

Phase Transformation Driving
Force.The phase transformation driving force (PT) can reflect the evolution extent of the population.In the proposed algorithm, phase transformation driving force is the factor of deciding phase transform.During the evolution, the energy possessed by a cloud particle is called its evolution potential.Phase transformation driving force is the evolution potential energy difference between the cloud particles of new generation and the cloud particles of old generation.The phase transformation driving force is defined as follows: The cloud evolution system realizes phase transformation if PT ≤ .EP represents the evolution potential of cloud particles.gen represents the maximum generation of the algorithm.bestindividual represents the best solution of the current generation.genindex is the number of the current evolution.

Condensation Operation.
The population realizes the global exploration in cloud gaseous phase by condensation operation.The population takes the best individual as the nucleus in the process of the condensation operation.In the cloud gaseous phase, the condensation growth space of cloud particle is calculated by the following equations.cd is the condensation factor.generation is the number of current evolution.   is the entropy of the th subpopulation.   is hyperentropy of the th subpopulation.Consider Generally, the optima value will act as the nucleus.However, the nucleus of the last subpopulation is the mean of the optima values coming from the rest of the subpopulations in reciprocity operation.The purpose of the reciprocity mechanism is to ensure the diversity of the cloud particles.For example, circle represents the location of the global optimum and pentagram represents the location of the local optima in Figure 8.

Solidification Operation.
The algorithm realizes convergence operation in cloud solid phase by solidification operation.sf is the solidification factor: where  = 1, 2, . . .,   and  = 1, 2, . . ., .   represents the number of the subpopulations. represents the number of dimensions.
The pseudocode of CPEA is illustrated in Algorithms 1, 2, 3, and 4.   represents the number of the subpopulation.Childnum represents the number of individuals of every subpopulation.The search space is [−, ].
Ali et al. [25] and Rahnamayan et al. [26].They are unimodal functions, multimodal functions, and rotate functions.All the algorithms are implemented on the same machine with a T2450 2.4 GHz CPU, 1.5 GB memory, and windows XP operating system with Matlab R2009b.

Parameter Study.
In this section, the setting basis for initial values of some parameters involved in the cloud particles evolution algorithm model is described.Firstly, we hope that all the cloud particles which are generated by the cloud generator can fall into the search space as much as possible.
(1) Initialize  (number of dimensions),   = 5; Childnum = 30;   =2/3r;   =   /1000; (2) for  = 1:   (3) for  = 1:  (4) if  is even (5)    = /2 (6) else ( 7) end if (9) end for (10) end for (11) Use cloud generator to produce cloud particles (12) Evaluate Subpopulation ( 13) Selectnucleus ( 14) while requirements are not satisfied do (15) Use cloud generator to produce cloud particles (16) Evaluate Subpopulation (17) Selectnucleus (18) Evolutionary algorithm of the cloud gaseous state (19) Evolutionary algorithm of the cloud liquid state (20) Evolutionary algorithm of the cloud solid state (21) end while Algorithm 4: Cloud particles evolution algorithm: CPEA.Secondly, the diversity of the population can be ensured effectively if the cloud particles are distributed in the search space uniformly.Therefore, some experiments are designed to analyze the feasibility of the selected scheme in order to illustrate the rationality of the selected parameters.It can be seen from Figure 9 that the number of the valid cloud particles reduces as the   increases.The average number of valid cloud particles is basically the same when   is (1/3), (1/2), and (2/3), respectively.However, the number of valid cloud particles is relatively fewer when   is (3/4), (4/5), and .It can be seen from Figure 10 that all the cloud particles fall into the research space when   is (1/3), (1/2), and (2/3), respectively.However, the cloud particles concentrate in a certain area of the search space when   is (1/3) or (1/2).By contrast, the cloud particles are widely distributed throughout all the search space [−2, 2] when   is (2/3).Therefore, set   = (2/3).The dashed rectangle in Figure 10 represents the decision variables space.

Definition of 𝐻 𝑒 . Ackley Function (multimode) and
Sphere Function (unimode) are selected to calculate the number of fitness evaluations (NFES) of the algorithm when the solution reaches  stop with different   .Let the search space  be [−32, 32] (Ackley) and [−100, 100] (Sphere),  stop = 10 −3 . 1 =   ,  2 = 0.1  ,  3 = 0.01  ,  4 = 0.001  ,  5 = 0.0001  ,  6 = 0.00001  , and  7 = 0.000001  .The reported values are the average of the results for 50 independent runs.As shown in Figure 11, for Ackley Function, the number of fitness evaluations is least when   is equal to 0.001  .For Sphere Function, the number of fitness evaluations is less when   is equal to 0.001  .Therefore,   is equal to 0.001  in the cloud particles evolution algorithm.
3.1.4.Definition of /.Many things in nature have the mathematical proportion relationship between their parts.If the substance is divided into two parts, the ratio of the whole and the larger part is equal to the ratio of the larger part to the smaller part; this ratio is defined as golden ratio, expressed algebraically in which  represents the golden ration.Its value is The golden section has a wide application not only in the fields of painting, sculpture, music, architecture, and management but also in nature.Studies have shown that plant leaves, branches, or petals are distributed in accordance with the golden ratio in order to survive in the wind, frost, rain, and snow.It is the evolutionary result or the best solution which adapts its own growth after hundreds of millions of years of long-term evolutionary process.
In this model, competition and reciprocity of the population depend on the fitness and are a dynamic relationship which constantly adjust and change.With competition and reciprocity mechanism, the convergence speed and the population diversity are improved.Premature convergence is avoided, and large search space and complex computational problems are solved.Indirect reciprocity is selected in this paper.Here, set / = (1 − 0.618) ≈ 0.38.

Comparison Strategies.
Five comparison strategies come from the literature [33] and literature [26] to evaluate the algorithms.The comparison strategies are described as follows.
Successful Run.A run during which the algorithm achieves the fixed accuracy level within the max number of function evaluations for the particular dimension is shown as follows: Error.The error of a solution  is defined as () − ( * ), where  * is the global minimum of the function.The minimum error is written when the max number of function evaluations is reached in 30 runs.Both the average and standard deviation of the error values are calculated.
Convergence Graphs.The convergence graphs show the median performance of the total runs with termination by either the max number of function evaluations or the termination error value.

Number of Function Evaluations (NFES).
The number of function evaluations is recorded when the value-to-reach is reached.The average and standard deviation of the number of function evaluations are calculated.

Acceleration Rate (AR).
The acceleration rate is defined as follows, based on the NFES for the two algorithms: in which AR > 1 means Algorithm 2 (2) is faster than Algorithm 1 (1).
The average acceleration rate (AR ave ) and the average success rate (SR ave ) over  test functions are calculated as follows: while it obtains relatively bad results in other functions.The reason for this phenomenon is that the changing rule of   is set according to the priori information and background knowledge instead of adaptive mechanism in CPEA.Therefore, the results of some functions are relatively worse.Besides,   which can be optimized in order to get better results is an important controlling parameter in the algorithm.

Comparison of JADE, DE, PSO, ABC, CEBA, and CPEA.
This section compares JADE, DE, PSO, ABC, CEBA, and CPEA in terms of SR and NFES for most functions at  = 30 except  03 ( = 2).The termination criterion is to find a value smaller than 10 −8 .The best result of different algorithms for each function is highlighted in boldface.
As shown in Table 3, CPEA can reach the specified precision successfully except for function  14 which successful run is 96.In addition, CPEA can reach the specified precision faster than other algorithms in function  01 ,  02 ,  04 ,  11 ,  14 , and  17 .For other functions, CPEA spends relatively long time to reach the precision.The reason for this phenomenon is that, for some functions, the changing rule of   given in the algorithm makes CPEA to spend more time to get rid of the local optima and converge to the specified precision.
As shown in Table 4, the solution obtained by CPEA is closer to the best solution compared with other algorithms except for functions  13 and  20 .The reason is that the phase transformation mechanism introduced in the algorithm can improve the convergence of the population greatly.Besides, the reciprocity mechanism of the algorithm can improve the adaptation and the diversity of the population in case the algorithm falls into the local optimum.Therefore, the algorithm can obtain better results in these functions.
As shown in Table 5, CPEA can find the optimal solution for the test functions.However, the stability of CPEA is worse than JADE and DE ( 25 ,  26 ,  27 ,  28 and  31 ).This result shows that the setting of   needs to be further improved.

Comparison of CMAES, JADE, DE, PSO, ABC, CEBA, and CPEA.
The convergence curves of different algorithms on some selected functions are plotted in Figures 13 and 14.The solutions obtained by CPEA are better than other algorithms with respect to functions  01 ,  11 ,  16 ,  32 , and  33 .However, the solution obtained by CPEA is worse than CMA-ES with respect to function  04 .The solution obtained by CPEA is worse than JADE for high dimensional problems such as function  05 and function  07 .The reason is that the setting of initial population of CPEA is insufficient for high dimensional problems.

Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and
CPEA.The algorithms DE, ODE, RDE, PSO, ABC, CBEA, and CPEA are compared in terms of convergence speed.The termination criterion is to find error values smaller than 10 −8 .The results of solving the benchmark functions are given in Tables 6 and 7. Data of all functions in DE, ODE, and RDE are from literature [26].The best result of NEFS and the SR for each function are highlighted in boldface.The average success rates and the average acceleration rate on test functions are shown in the last row of the tables.

Test on the CEC2013 Benchmark Problems.
To further verify the performance of CPEA, a set of recently proposed CEC2013 rotated benchmark problems are used.In the experiment, -test [37] has been carried out to show the differences between CPEA and the other algorithm."General merit over contender" displays the difference between the number of better results and the number of worse results, which is used to give an overall comparison between the two algorithms.Table 8, we can see that PSO, ABC, PSO-FDR, PSO-cf-Local, and CPEA work well for  01 , where PSO-FDR, PSO-cf-Local, and CPEA can always achieve the optimal solution in each run.CPEA outperforms other algorithms on  01 ,  02 ,  03 ,  07 ,  08 ,  10 ,  12 ,  15 ,  16 , and  18 ; the cause of the outstanding performance may be because of the phase transformation mechanism, which leads to faster convergence in the late stage of evolution.

Conclusions
This paper presents a new optimization method called cloud particles evolution algorithm (CPEA).The fundamental concepts and ideas come from the formation of the cloud and the phase transformation of the substances in nature.In this paper, CPEA solves 40 benchmark optimization problems which have difficulties in discontinuity, nonconvexity, multimode, and high-dimensionality.The evolutionary process of the cloud particles is divided into three states: the cloud gaseous state, the cloud liquid state, and the cloud solid state.The algorithm realizes the change of the state of the cloud with the phase transformation driving force according to the evolutionary process.The algorithm reduces the population number appropriately and increases the individuals of subpopulations in each phase transformation.Therefore, CPEA can improve the whole exploitation ability of the population greatly and enhance the convergence.There is a larger adjustment for   to improve the exploration ability of the population in melting or gasification operation.When the algorithm falls into the local optimum, it improves the diversity of the population through the indirect reciprocity mechanism to solve this problem.
As we can see from experiments, it is helpful to improve the performance of the algorithm with appropriate parameter settings and change strategy in the evolutionary process.As we all know, the benefit is proportional to the risk for solving the practical engineering application problems.Generally, the better the benefit is, the higher the risk is.We always hope that the benefit is increased and the risk is reduced or controlled at the same time.In general, the adaptive system adapts to the dynamic characteristics change of the object and disturbance by correcting its own characteristics.Therefore, the benefit is increased.However, the adaptive mechanism is achieved at the expense of computational efficiency.Consequently, the settings of parameter should be based on the actual situation.The settings of parameters should be based on the prior information and background knowledge in some certain situations, or based on adaptive mechanism for solving different problems in other situations.The adaptive mechanism will not be selected if the static settings can ensure to generate more benefit and control the risk.Inversely, if the adaptive mechanism can only ensure to generate more benefit while it cannot control the risk, we will not choose it.The adaptive mechanism should be selected only if it can not only increase the benefit but also reduce the risk.In CPEA,   and   are statically set according to the priori information and background knowledge.As we can see from experiments, the results from some functions are better while the results from the others are relatively worse among the 40 tested functions.
As previously mentioned, parameter setting is an important focus.Then, the next step is the adaptive settings of   and   .Through adaptive behavior, these parameters can automatically modify the calculation direction according to the current situation while the calculation is running.Moreover, possible directions for future work include improving the computational efficiency and the quality of solutions of CPEA.Finally, CPEA will be applied to solve large scale optimization problems and engineering design problems."+, " "−, " and "≈" denote that the performance of CPEA ( = 10) is better than, worse than, and similar to that of other algorithms.
Figure 8(a) displays the position of nucleus before reciprocity operator.Figure 8(b) shows the individuals generated by cloud generator based on Figure 8(a).
Figure 8(c) shows the position of nucleus with reciprocity operator based on Figure 8(a).

Figure 8 (
d) shows the individuals generated by cloud generator based on Figure 8(c).By comparing Figure 8(b) with Figure 8(d), it can be seen that the distribution of the individuals in Figure 8(d) is more widely than those in Figure 8(b).

Figure 9 :
Figure 9: Valid cloud particles of different   .

3. 1 . 1 .
Definition of   .Let the search range  be [−, ] and  is the search radius.Set  = 2,  1 = −2/3,  2 = 2/3,   = 0.001  .Five hundred cloud particles are generated in every experiment.The experiment is independently simulated 30 times.The cloud particle is called a valid cloud particle if it falls into the search space.
Figure 9 is a graph which shows the number of the valid cloud particles generated with different   , and Figure 10 is a graph which shows the distribution of valid cloud particles generated with different   .

Figure 10 :Figure 11 :
Figure 10: Distribution of cloud particles in different   .
The cost-to-benefit ratio of the altruistic act is (1−0.618)≈ 0.38 and 0.618 is the golden ratio.  represents the number of the subpopulations.  is the nucleus of the th subpopulation.Consider

Table 1 :
[35] section compares CMA-ES, DE, RES, LOS, ABC, CEBA[34], and CPEA in terms of mean number of function evaluations for successful runs.Results of CMA-ES, DE, RES, and LOS are taken from literature[35].Each run is stopped and regarded as successful if the function value is smaller than  stop.As shown in Table 1, CPEA can reach  stop faster than other algorithms with respect to  07 and  17 .However, CPEA converges slower than CMA-ES with respect to  18 .When the Average number of function evaluations to reach  stop of CPEA versus CMA-ES, DE, RES, LOS, ABC, and CEBA on  07 ,  17 ,  18 , and  33 function.The bold entities indicate the best results obtained by the algorithm. represents dimension.dimension is 30, CPEA can reach  stop relatively fast with respect to  33 .Nonetheless, CPEA converges slower while the dimension increases.
[36]Comparison of CMA-ES, DE, PSO, ABC, CEBA, andCPEA.This section compares the error of a solution of CMA-ES, DE, PSO, ABC[36], CEBA, and CPEA under a given NFES.The best solution of every single function in different algorithms is highlighted with boldface.As shown in

Table 2 :
Comparison on the ERROR values of CPEA versus CMA-ES, DE, PSO, ABC, and CEBA on function  01 ∼ 13 .

Table 8
gives the information about the average error, standard deviation, and -value of 30 runs of 7 algorithms over 300,000 FEs on 20 test functions with  = 10.The best results among the algorithms are shown in bold.Table 8 consists of unimodal problems ( 01 - 05 ) and basic multimodal problems ( 06 - 20 ).From the results of

Table 3 :
Comparison on the SR and the NFES of CPEA versus JADE, DE, PSO, ABC, and CEBA for each function.

Table 6 :
Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA.The best result for each case is highlighted in boldface.

Table 7 :
Comparison of DE, ODE, RDE, PSO, ABC, CEBA, and CPEA.The best result for each case is highlighted in boldface.

Table 8 :
Results for CEC2013 benchmark functions over 30 independent runs.