A Hybrid Mutation Chemical Reaction Optimization Algorithm for Global Numerical Optimization

This paper proposes a hybrid metaheuristic approach that improves global numerical optimization by increasing optimal quality and accelerating convergence. This algorithm involves a recently developed process for chemical reaction optimization and two adjustment operators (turningand mutation operators). Threetypes of mutation operators (uniform,nonuniform, and polynomial) werecombinedwithchemicalreactionoptimizationandturningoperatortofindthemostappropriateframework.Thebestsolution amongthesethreeoptionswasselectedtobeahybridmutationchemicalreactionoptimizationalgorithmforglobalnumerical optimization.Theoptimalquality,convergencespeed,andstatisticalhypothesistestingofouralgorithmaresuperiortothose previoushighperformancealgorithmssuchasRCCRO,HP-CRO2,andOCRO.


Introduction
Various optimization approaches have been used to solve existing global optimization problems in science, engineering, and related fields.Evolutionary algorithms (EAs), inspired by the biological evolution and social behavior of organisms, are the most powerful among optimization frameworks.Several EAs have been successfully used to solve global optimization problems.Genetic algorithms (GAs) [1,2] were developed based on the natural process of evolution.Particle swarm optimization (PSO) [3] was inspired by the social behavior of bird flocking or fish schooling.Memetic algorithms (MAs) [4] are similar to GAs except that MAs initiate a local search before undergoing the evolutionary process [5].Ant colony optimization [6] resembles PSO based on the capability of ants to find the shortest route between their nest and a food source.Artificial bee colony [7,8] was inspired by the echolocation behavior of bees.Harmony search [9] was based on the natural musical performance of a musician searching for an optimal state of harmony.
A recent study has proposed a new chemical reactioninspired EA named "chemical reaction optimization" [10,11].CRO is a low-latency algorithm that mimics molecular interactions in chemical reactions to reach a low-energy stable state.An effective version of CRO is the real-coded chemical reaction optimization (RCCRO) [12], which can process a large set of continuous problems [13][14][15][16][17]. HP-CRO, a hybrid algorithm based on PSO and CRO, produces more distinguished results compared with RCCRO [18].However, orthogonal chemical reaction optimization (OCRO) algorithms for global numerical optimization problems [19] are the best when compared with RCCRO and HP-CRO.OCRO is efficient in solving high-dimension functions but not in solving low-dimensional functions.
The current study proposes a hybrid mutation chemical reaction optimization algorithm for global numerical optimization (MCRO).This strategy can be used to solve high-dimension and low-dimension functions better and faster.The proposed framework combines CRO with a turning operator and a mutation operator.Our algorithm was determined by choosing the best solution that estimates the performance of the combination with three types of mutation operators (uniform, nonuniform, and polynomial).
This paper is organized as follows.Section 2 provides a brief review of preliminary studies on CRO and mutation 2 Mathematical Problems in Engineering operator.Section 3 presents the proposed approach: MCRO.Section 4 describes the evaluation methods and simulation results and discusses the comparisons of algorithms.Section 5 concludes this paper and suggests proposals for future work.

Preliminary
2.1.Optimization.In science, mathematics, engineering, and other related fields, optimization is the selection of an optimal solution from a set of feasible alternatives.In general, an optimization problem includes minimizing or maximizing a function by systematically selecting input values from a given feasible set [20].The problem function () is scalar, where a variable  demonstrates a particular solution and  is usually a vector of  components.The  components specify the dimensions of .An optimization problem can be subjected to a vector of constraints  = { 1 ,  2 , . . .,   } that limits the feasible region, where  corresponds to the total number of constraints [12].The aim of optimization is to generate the optimal solution and manage the convergence speed of a problem function.Our research focuses on function minimization; the goal of minimum optimization is to find the minimum solution ś ∈  and ( ś ) ≤ (), ∀ ∈ .

CRO.
In principle, chemical reactions are governed by two rules of thermodynamics.First, energy cannot be created or destroyed but can be transformed from one entity to another; second, the entropy of a system tends to increase, where entropy is the measure of the degree of disorder [21].CRO was inspired by the nature of chemical reactions and mimics the interactions of molecules in a chemical reaction to reach a low-energy stable state.Potential energy is the energy stored in a molecule with respect to its molecular configuration; the system becomes disordered when potential energy is converted to other forms [11].Molecules stored in a container are vital to the manipulation of agents.Each molecule contains a profile that includes several attributes, such as molecular structure (), current potential energy (PE), and current kinetic energy (KE).In CRO reaction, the initial reactants in high-energy states incur a sequence of collisions.Molecules collide either with other molecules or with the walls of the container, pass through energy barriers, and become the final products in low-energy stable states.Four types of elementary CRO reactions capture the transition of high-energy molecules to stable states [10].
On-wall ineffective collision is the reaction created when a molecule hits the wall of a container and then bounces back.This reaction only slightly changes the molecular structure  when PE  + KE  ≥ PE ώ occurs.The modification of molecular structure ώ is processed by neighborhood search operator as ώ \in neighborhood ().The central energy buffer is updated by extracting and storing a certain portion of its KE.The profile of the molecule is updated as KE ώ = (PE  − PE ώ + KE  ) × , where  is a random number that  ∈ [, 1].The new PE is calculated using the new  or ώ .
Intermolecular ineffective collision refers to two or more molecules that collide with each other and then separate.The profiles of the molecules and the central energy buffer are updated when PE The number of molecules is unaltered after the collision, but the molecular structures are produced from their own neighborhoods by neighborhood search operator.Similar to on-wall ineffective collision, this reaction also slightly changes the molecular structure.
Decomposition occurs when a molecule hits the wall of a container and then splits into two or more molecules.This elementary reaction is applied to finish local search and to explore other regions.The profiles of the molecules and the central energy buffer are updated when PE  + KE  ≥ PE ώ 1 + PE ώ 2 and when the energy buffer is sufficient.This reaction significantly alters the molecular structures of the resultant molecules.
Synthesis represents the situation when two or more molecules collide and mingle to a single molecule.The profiles of the molecules and the central energy buffer are updated when This reaction strongly and significantly alters the resultant molecular structure.
The algorithm of CRO is shown as Algorithm 1 [12].On-wall and intermolecular ineffective collisions are local searches, whereas decomposition and synthesis are global searches.

Mutation Operator.
Mutation is a natural process that changes a DNA sequence to a new variation endowing individual survival advantages in biology.Mutation operator, which has been applied to GAs, simulates the natural biological evolution of species by exchanging information and producing offspring chromosomes [1,2].In GAs, mutation operator is designed to work on an individual chosen; this individual is randomly selected or chosen by defined condition.The new gene value is modified based on userdefinable probability or other conditions such as a distribution.GAs containing mutation can widely maintain offspring chromosomes, instead of being limited to the genes available in the initial population.Therefor GAs with mutation can generate better solution and prevent premature convergence compared to GAs without mutation.
The concept of MCRO is initialed by supposing that advantage of mutation operator would empower global numerical optimization while combining with CRO.Then we reformed the mutation operator algorithm in accordance with CRO infrastructure, such as replacing an individual in GAs by a molecule (), assigning probability and distribution similar to user-definable probability or condition in GAs, and changing gene value of GAs with molecular structure ().In CRO, a molecular structure contains amount number of moles and a pair of moles is presented as a set of  = { 1 ,  2 , . . .,  V }, where V is the number of moles.Mutation operator algorithm generates the new molecular structure ( ώ ) based on mutation operator formula relating to mutation type.There are many types of mutation and these types depend on the representation itself.Several mutation Input: Objective function , constraints, and the dimensions of the problem (1) / * initialization phase * / (2) Assign parameter value to PopSize, InitialKE, StepSize, buffer, KELossRate, On-wallColl, (3) DecThres, SynThres, The set of molecules in this Container are molecules 1, 2, . . ., PopSize (4) for each of the molecule do (5) Assign a random solution to the molecular structure .Calculate the PE by () and evaluate.(6) Assign the KE with InitialKE (7) end for (8) / * Iterations phase * / (9) while (the stopping criteria not met) do (10) Get  randomly in interval [0, 1] (11) if ( > On-wallColl) then (12) select a molecule  from Container Randomly (13) if (Decomposition criterion met) then (14) Decomposition (15) if Success then Add a new molecule   to Container (16) else (17) On-wallIneffectiveCollision (18) end if (19) else (20) select two molecule  1 and  2 from Container Randomly (21) if (Synthesis criterion met) then (22) Synthesis (23) if Success then Remove molecule  2 from Container (24) else (25) Inter-molecularIneffectiveCollision (26) end if (27) end if (28) evaluate and keep any new minimum solution (29) end while (30) / * The Output phase * / (31) Output the best solution and its function value.

Proposed Approach: MCRO
Mutation chemical reaction optimization (MCRO), which is the improved hybrid of CRO and two adjustment operators, generates various molecular structures to produce highquality results and accelerate convergence.The turning operator and the mutation operator will be explained in this section.

Turning Operator.
Turning operator is a new operator that is merged into a subalgorithm named neighborhood search operator or (⋅).Neighborhood search operator is used in three types of elementary reactions (on-wall ineffective collision, intermolecular ineffective collision, and decomposition) to transform the molecular structure from the neighborhood of the operand [12].The new molecular structures of the solution are calculated as Formula (1), and turning operator is generated once for an objective function (  ) by using Formula (2).This new operator can highly improve the optimal quality and reliability of the algorithm: Turning operator of   = { random number interval [1, 1 − 300]} , (1 ≤  ≤ 23) . (2)

Mutation
Operator.Mutation operator is migrated into MCRO to improve solution diversity and accelerate convergence by generating various molecular structures.These structures are generated by using the processing mutation formula.This action increases the probability of finding the optimal solution and avoids being trapped into a local optimal solution; therefore, it accelerates convergence.In our research, the best mutation operator that matches our algorithm was selected from three types of mutation operators (uniform, nonuniform, and polynomial).
Uniform and nonuniform mutations are presented in basic GAs and their extending development of algorithm [22,25].Polynomial mutation is a popular one that was first introduced in Nondominated Sorting Genetic Algorithm (NSGA) and NSGA-II [23,24].Uniform, nonuniform, and polynomial mutations are calculated as Formulas (3), (4), and (5), respectively.As mutation operator algorithm mentioned in Section 2.3, input data is a set of molecular structures presented as  = { 1 ,  2 , . . .,  V }, where V is the number of moles,   is a member of , 1 ⩽  ⩽ V. Likewise output data is a set of molecular structures presented as ώ .According to opportuneness, the variable names such as  and ώ maybe present differently when they appear in other algorithms such as ,  ώ ,  1 , and  2 : Generally, the optimization method consists of three phases: initialization, iteration, and output, which are all included in our proposed approach.We applied mutation operator to the initialization and iteration phases of the algorithm.
A mutation operator was applied to each new molecule in the initialization phase.In the iteration phase, we processed mutation operator in all types of elementary reaction, namely, on-wall ineffective collision, intermolecular ineffective collision, decomposition, and synthesis.Mutation operator was incorporated to change the molecular structure.We compared the results before and after the performance of mutation operator and selected the better result.We believe mutation is one potential component to improve the performance of MCRO, since mutation operator can spread the search space by randomly sampling new points.Furthermore this mutation operator algorithm increases the chance of generating more powerful result not less than twice of original CRO for every elementary reaction.Following are the four elementary reactions' explanations.
In on-wall ineffective collision of MCRO, a molecule  is chosen.The first new molecular structure ( ώ ) of  is generated by neighborhood search operator which includes turning operator as mentioned in Section 3.1 and then the first fitness or PE ώ is calculated by ( ώ ).The new KE  is generated and buffer is updated when the condition PE  + KE  ≥ PE ώ occurs.Next the second new molecular structure () is produced by computing mutation operator algorithm and then the second fitness or tempPE is calculated Input: a molecule  and buffer (1) Obtain ώ = ℎ() with turning operator Calculate the tempPE by () (9) if tempPE better than PE ώ then Replace ώ with , PE ώ with tempPE (10) Update the profile of  by  = ώ , PE  = PE ώ and KE  = KE ώ (11) end if Output  and  Algorithm 3: On-wall ineffective collision of MCRO.by ().The profile of molecule M is updated as the best solution by comparing the first fitness (PE ώ) and the second fitness (tempPE).Therefore this action contains two probabilities and achieves the best solution for each time while there is only one choice for original CRO.The on-wall ineffective collision of MCRO is shown in Algorithm 3.
Intermolecular ineffective collision of MCRO algorithm starts by selecting two molecules  1 and  2 .Subsequently, neighborhood search operator with turning operator transforms the molecular structures of  1 and  2 from  1 and  2 to ώ 1 and ώ 2 , respectively, and evaluates PE ώ 1 by ( ώ 1 ) and PE ώ 2 by ( ώ 2 ).The new KE are generated as KE ώ Furthermore, the algorithm executes mutation operator to build two new molecular structures,  1 and  2 , and processes fitness of two molecules tempPE 1 by ( 1 ) and tempPE 2 by ( 2 ).Meanwhile, the best fitness between PE ώ 1 and tempPE 1 is selected to be the optimal solution of  1 , similar to the optimal solution of  2 is the best fitness between PE ώ 2 and tempPE 2 , and then the profiles of molecules  1 and  2 are updated.Prove that there are two choices to estimate the most outstanding result for this reaction that acted to a molecule; therefor this elementary reaction of MCRO consists of four chances to meet the optimum solution for a time, while the same reaction of original CRO is just giving two choices (one choice for a molecule).Algorithm 4 illustrates intermolecular ineffective collision of MCRO algorithm.
Decomposition of MCRO algorithm begins by selecting a molecules  from container.Copy a molecular structure () of  to two new variables  1 and  2 .Next the molecular structures  1 and  2 are transformed to ώ 1 and ώ 2 , respectively, by computing neighborhood search operator (with turning operator), and furthermore evaluate the fitness PE ώ 1 as ( ώ 1 ) and the fitness PE ώ 2 as ( ώ 2 ).The profiles of molecules are updated when Create new molecules   (10) else if  1 +  ≥ 0 then (11) Get  1 ,  2 ,  3 , and Create new molecules   (16) end if (17)  (Success = TRUE), the two new molecular structures  1 and  2 are established by operated mutation operator algorithm and then the fitness of tempPE 1 and tempPE 2 are evaluated by f (temp 1 ) and f (temp 2 ), respectively.The first minimum fitness chooses the best one from PE ώ 1 or tempPE 1 , and the second minimum fitness selects the best one between PE ώ 2 and tempPE 2 .Furthermore profiles of  and   are updated.This concludes that there are four opportunities to reach the best result in decomposition of MCRO similar to intermolecular ineffective collision of MCRO; in other words, they are twice the same reaction in original CRO.Decomposition of MCRO is represented as Algorithm 5.
In synthesis of MCRO algorithm, two molecules  1 and  2 are chosen.Secondly, neighborhood search operator (with turning operator) generated the two molecular structures of  1 and  2 from  1 and  2 to ώ 1 and ώ 2 , respectively, and the first fitness or PE ώ 1 is calculated by ( ώ 1 ).Subsequently, the new KE ώ 1 is computed when Afterward mutation operator algorithm generates ώ 1 to  ώ 1 , and the fitness generated by mutation operator algorithm or tempPE is calculated by ( ώ 1 ).Finally the profile of molecule  1 is updated by choosing the best solution between the fitness (PE ώ 1 ) and the fitness generated by mutation operator (tempPE).There are two probabilities included in this elementary reaction while synthesis of original CRO provides only one.Synthesis of MCRO algorithm is presented as Algorithm 6.

Simulation Results and Discussion
4.1.Benchmarks and Parameters.The proposed MCRO was simulated to solve the 23 objective functions based on RCCRO [12].The benchmark objective functions are classified into three categories, as shown in Table 1.

Category I (High-Dimensional Unimodal Functions).
Category I is composed of high-dimensional functions  1 - 7 and includes 30 molecules in each container.These functions are easy to solve because each objective function has only one global minimum.

Category II (High-Dimensional Multimodal Functions).
Category II provides high-dimensional functions  8 - 13 and includes 30 molecules in each container.These functions consist of many local minimums and can solve the most difficult problem.

Category III (Low-Dimensional Multimodal Functions).
Category III includes low-dimensional functions  14 - 23 Goldstein-Price function Hartman's family Shekel's family do Mutation of ώ 1 to  ώ 1 (7) Calculate the tempPE by ( ώ The mutation operator parameters in our experiment are probability and DistributionIndex.We set probability = 0.6 and DistributionIndex = 0.9 for every objective function.

MCRO Solution.
To obtain the most appropriate MCRO solution, we considered three types of mutation operators: polynomial, uniform, and nonuniform.Our approach estimates start with alloyed CRO with turning operator, consequently integrating mutation operator to the algorithm.The first, second, and third options are the algorithms with polynomial, uniform, and nonuniform mutations, respectively.We generated the results of these three solutions by running each objective function 25 times.The optimal quality of the solution of each function is the mean of the best global minimum at 25 repeats.If the means of the results are equal, the second key for comparison is the standard deviation.All simulations included in Section 4.3 are performed on the same personal computer and same computing environment with computer Intel Core i7, CPU 3.10 GHz, and Ram 4 GB.
Table 3 compares the results of Category I functions for the three solutions.All of the solutions yield good results for  1 ,  2 ,  3 ,  4 , and  6 .The algorithm with polynomial mutation is the leader for  5 , and the algorithm with uniform mutation generates the best result for  7 .The average result of Category I functions is the same for the three solutions.Table 4 compares the results of Category II functions for the three solutions.The results of the three solutions are as good as those of the two other solutions in  9 ,  10 , and  11 .The solution with polynomial mutation performs best for  12 and  13 .The solution with nonuniform mutation generates the best result for  8 .The average result of highdimensional multimodal functions is led by the solution combined with polynomial mutation, nonuniform mutation, and then uniform mutation.
Table 5 compares the results of Category III functions.The solution with polynomial mutation computes the best result for  14 ,  16 ,  17 ,  18 , and  19 .Meanwhile, the solution with uniform mutation is the leader for  15 ,  20 ,  21 ,  22 , and  23 , and the solution with nonuniform mutation performs best for 15 .The polynomial solution does not generate the best result for the remaining functions, but it produces close results to those of the leader for the first five functions.The average ranking of Category III functions for the three solutions in order of solution is polynomial mutation, uniform mutation, and nonuniform mutation.
The overall ranking is shown in the bottom of Table 5.The solution with polynomial mutation derives the most excellent results, followed successively by uniform mutation and nonuniform mutation.Hence, we prefer polynomial mutation operator for MCRO.

Experimental Method and Results
. The aim of the present experiment is to evaluate the performance of MCRO by comparing it with three previously published algorithms: RCCRO4 [12], the best version of RCCRO; HP-CRO2 [18], the best version of HP-CRO and a hybrid of CRO and PSO; and OCRO [19], which represents the best results among the three algorithms.We generated the result for MCRO, OCRO, HP-CRO2, and RCCRO4 by running each objective function 25 times.We provided three experiment methods to evaluate the capability of comparison algorithms.These three methods were evaluated for optimal solution quality, convergence speed, and statistical hypothesis testing.
(1) Optimal Solution Quality Evaluation.For each objective function, the optimal solution or the best global minimum is the most important result in each runtime.In a computing time, each round of iteration contains the local minimum as a result.Consequently, if the current global minimum is worse than the local minimum, then the current global minimum is replaced by the local minimum.The best global minimum at each computing time is generated when the program meets the stop condition of MCRO and other optimization algorithms.The mean of the best global minimum at 25 computing times is the optimal solution of each function.If the optimal solutions of the competitors are equal, then the second key for comparison is the standard deviation.The following optimal solution quality evaluations are demonstrated.Table 6 presents the optimal solution quality evaluation of MCRO, OCRO, HP-CRO2, and RCCRO4 for highdimensional unimodal functions or Category I ( 1 - 7 ).MCRO performs best for  1 ,  2 ,  3 ,  4 ,  6 , and  7 but not for  5 .MCRO and OCRO are the best for  1 ,  2 ,  3 , and  6 ; HP-CRO2 and RCCRO4 derive the same results as MCRO and OCRO for  6 .RCCRO4 generates the best result for  5 .The ranking of optimal solution quality for Category I functions in descending order is as follows: MCRO, OCRO, HP-CRO2, and RCCRO4.
Table 7 compares the optimal solution quality for Category II functions.MCRO performs best for  9 ,  10 ,  11 , and  12 ; OCRO generates the same results as MCRO for  9 ,  10 , and  11 ; HP-CRO2 is the best for  13 ; RCCRO4 has the best result for  8 .The ranking of optimal solution quality for Category II functions in descending order is the same as that for Category I.
Table 8 compares the results of Category III functions.MCRO is the best in this category because it generates the best results for 8/10 functions:  14 ,  16 ,  17 ,  18 ,  19 ,  21 ,  22 , and  23 .MCRO is followed by OCRO, which is the best for  15 .HP-CRO2, which generates the best result for  20 , is in third rank, followed by RCCRO4 in fourth rank.The ranking of optimal solution quality for Category III functions in descending order is the same as that for Categories I and II.
The overall ranking presented at the bottom of Table 8 shows that MCRO performs best in the optimal solution quality evaluation.MCRO is followed successively by OCRO, HP-CRO2, and RCCRO4.
(2) Convergence Speed Evaluation.Convergence speed is an essential issue that indicates the performance of algorithms in optimization.In our experiment, the convergence speed of the algorithm was calculated by counting the number of iterations (FEs) before the algorithm converges with the accepted result.The accepted result is the results given there are the average FEs needed to reach the threshold expressed as acceptable solutions for each objective function such as acceptable solution of  1 is 1 − 60, the robust algorithm has to generate the result at lease not worse than accepted result.The algorithm with fewer FEs is more powerful than that with greater FEs.After evaluating the convergence speed on the basis of iteration number, we evaluated the same parameter in MCRO, OCRO, HP-CRO2, and RCCRO4 by drawing a convergence curve for specific functions ( 4 ,  7 ,  9 ,  10 ,  11 , and  23 ) in a particular run.MCRO remains the most advantageous among the four algorithms.
(3) Statistical Hypothesis Testing.Friedman test [26], a wellknown nonparametric statistical test, was used to assess the statistical hypothesis.Friedman test evaluates the differences or equality of results among competitor algorithms.We calculated the Friedman statistic according to  2 distribution with  − 1 degrees of freedom.Then, we estimated the  value according to normal approximations.The Friedman statistic   is presented as Formula (6).Parameter  is the number of objective functions, which is 23 in our test. represents the number of comparison algorithms.In this case, we compared our proposed algorithm MCRO with three previous algorithms (OCRO, HP-CRO2, and RCCRO4); hence  = 4. Friedman rank or   is searching by transforming the quality results of MCRO, OCRO, HP-CRO2, and RCCRO4 to ranks for each objective function ().In cases of equal ranks, the average ranks are assigned: Table 10 demonstrates the results of Friedman rank test in terms of the comparison of optimal solution quality and convergence speed.MCRO achieves the best rank in the statistic test comparing both the optimal solution quality and the convergence speed on the basis of the corresponding Friedman ranks 1.59 and 1.22, the statistic   values 21.76398 and 35.96694, and the  values 7.30387 − 05 and 7.60987 − 08.In addition, the results explicitly display significant differences across the competing algorithms.In multiple comparisons, the post hoc test was utilized to compare a set of algorithms with a control algorithm (best algorithm).The post hoc test starts with finding the test statistic  of the th and th algorithms, which is processed using Formula (7).
where   and   are the Friedman ranks of the compared algorithms.
The unadjusted  value was calculated using the  value from the table of normal distribution (0, 1).The  value was then compared with an appropriate level of significance .
In the present study, the following three computing procedures were applied to obtain the adjusted  values (APVs): Bonferroni-Dunn procedure [27] presented as Formula (8), Holm procedure [28] assigned as Formula (9), and Hochberg procedure [29] expressed as Formula (10).Subscripts  and  are indexes that each corresponds to a concrete hypothesis, and   is the -value obtained for the th hypothesis.
The Bonferroni-Dunn procedure adjusts the value  through a single step:   The results of the APVs for the Friedman test are shown in Table 11, where the control algorithm is the proposed algorithm MCRO.These results comprise the  value, unadjusted  values, and APVs that were extracted using the post hoc procedures and computed using the three procedures.We compared the results for both optimal solution quality and convergence speed.MCRO is superior to RCCRO4, HP-CRO2, and OCRO in all post hoc procedures considered.
(4) Computation Time.Normally, the numbers of computation times are often proportional to capability of computer.Unlike the other results such as optimal solution quality and convergence speed, although the simulations perform on various personal computers, these results are very approximate.Hence our experiments disregard using the computation time as an evaluation method.However, the computation times appearing in Table 12 indicate that, in our simulated situation, the computer times for these four comparison algorithms are closed.
4.4.Discussion.Section 4.2 recommends that the best MCRO solution is the combination of CRO, turning operator, and polynomial mutation operator.The solution with polynomial mutation obtains the most appropriate results over uniform and nonuniform mutations.
The experimental results in Section 4.3 substantiate that MCRO is more powerful than OCRO, HP-CRO2, and RCCRO4.MCRO can generate the best result on most of the 23 objective functions (see Figure 1).We determined several issues, such as optimal solution quality, convergence speed, and statistical hypothesis testing, to evaluate algorithm performance.Table 6 shows that MCRO can generate the best optimal solution for all high-dimensional unimodal functions.We can obtain the optimal solution 0 (equal  min ) and StdDev 0 for five functions ( 1 ,  2 ,  3 ,  4 , and  6 ).MCRO achieves the best average rank in solving high-dimensional multimodal functions (Table 7) and performs the optimal  ).For low-dimensional functions (Table 8), the results of MCRO are close to the true optimal, indicating that this algorithm is outstanding and consistent.Moreover, some comparison algorithms can generate the same results as MCRO in solving several objective functions, but MCRO converges faster than the other algorithms (Table 10).

Conclusions and Future Work
This paper presents a hybrid optimization algorithm MCRO that solves global numerical optimization problems by combining CRO with turning and mutation operators.The most excellent MCRO framework is the utilization of the solution with polynomial mutation.Therefore, we provide polynomial mutation operator into our proposed approach.Our algorithm is advantageous because the nature of CRO is suitable for global numerical optimization, and turning operator is a perfect method to improve optimal solution quality and algorithm reliability.Furthermore, polynomial mutation operator accelerates algorithm convergence.To analyze the evaluation results, we compared the results of MCRO and other previously published algorithms (OCRO, HP-CRO2, and RCCRO4) in terms of optimal solution quality, convergence speed, and statistical hypothesis testing.There results indicate that MCRO is superior to other algorithms in solving global numerical optimization problems.Several potential parameters such as Popsize, InitialKE, buffer, KELossRate, On-wallColl, DecThres, SynThre, StepSize, probability, and DistributionIndex were also discussed.MCRO is useful in solving global numerical optimization problems with high-dimensional unimodal, highdimensional multimodal, and low-dimensional functions.However, whether or not MCRO will also be excellent for solving other problem types is unclear.Hence, the applicability of MCRO to other problem types requires further research and consideration.
Our future work will focus on investigating the applicability of MCRO to other numerical optimization problems and practical engineering optimization problems.We also intend to develop a new hybrid metaheuristic approach to solve optimization problems.

Table 3 :
Comparisons of the result of the three solutions for Category I ( 1 - 7 ) (the best results are marked in bold).

Table 4 :
Comparisons of the result of the three solutions for Category II ( 8 - 13 ) (the best results are marked in bold).

Table 5 :
Comparisons of the result of the three solutions for Category III ( 14 - 23 ) (the best results are marked in bold).
Table 9 compares the convergence speeds of the algorithms.MCRO demonstrates the fastest convergence speed for all functions, except for  8 , where HP-CRO2 is the best, and for  19 and  21 , where OCRO is the best.In addition, the distance between MCRO and the second best of each function is significant.For  1 ,  3 ,  4 ,  10 ,  11 , and  14 , MCRO converges the accepted results with 691, 819, 338, 896, 1056, and 314 FEs, respectively, whereas OCRO, respectively, has 23253, 25053, 31245, 65200, 64947, and 2818 FEs for the same functions.For  22 , MCRO and RCCRO4 undergo 28 and 4721 FEs, respectively.The convergence speed ranking of comparison algorithms for all objective functions from fastest to slowest is as follows: MCRO, OCRO, HP-CRO2, and RCCRO4.

Table 6 :
Optimal solution quality for Category I ( 1 - 7 ) (the best results are marked in bold).

Table 7 :
Optimal solution quality for Category II ( 8 - 13 ) (the best results are marked in bold).

Table 8 :
Optimal solution quality for Category II ( 14 - 23 ) (the best results are marked in bold).

Table 9 :
Convergence speed (the best results are marked in bold).

Table 10 :
Result of Friedman test.

Table 11 :
Adjusted  values for the Friedman test (control algorithm: MCRO).

Table 12 :
Computation times. min ) and StdDev 0 for two functions ( 9 and  11