A Mutual-Evaluation Genetic Algorithm for Numerical and Routing Optimization

Many real-world problems can be formulated as numerical optimization with certain objective functions. However, these objective functions often contain numerous local optima, which could trap an algorithm frommoving toward the desired global solution. To improve the search efficiency of traditional genetic algorithms, this paper presents a mutual-evaluation genetic algorithm (MEGA). A novel mutual-evaluation approach is employed so that the merit of selected genes in a chromosome can be determined by comparing the fitness changes before and after interchanging with those in the mating chromosome. According to the determined genome merit, a therapy crossover can generate effective schemata to explore the solution space efficiently. The computational experiments for twelve numerical problems show that the MEGA can find near optimal solutions in all test benchmarks and achieve solutions with higher accuracy than those obtained by eight existing algorithms. This study also uses the MEGA to find optimal flow-allocation strategies for multipath-routing problems. Experiments on quality-of-service routing scenarios show that the MEGA can deal with these constrained routing problems effectively and efficiently. Therefore, the MEGA not only can reduce the effort of function analysis but also can deal with a wide spectrum of real-world problems.


Introduction
Many engineering optimization issues can be formulated as global optimization problems with numerical functions. When solving a complex problem, the particular challenge is that algorithms may be trapped in local optima and fail in finding global optima. Recently, genetic algorithms (GAs) have received considerable attention for solving complex and unstructured problems [1,2]. However, traditional genetic algorithm (TGA) often suffers from the drawbacks of premature convergence and weak exploitation capabilities [3].
To overcome the deficiencies of the TGA, this paper proposes a mutual-evaluation approach to incorporate with the TGA as a mutual-evaluation genetic algorithm (MEGA). The proposed therapy crossover can implicitly generate effective schemata to efficiently exploit the given search space without explicitly analyzing the solution space. The performance of the proposed MEGA is experimented on 12 well-known numerical functions and compared with four well-known evolutionary algorithms (EAs) and four existing modified GAs. The experimental results show that the proposed MEGA is able to increase the accuracy by several orders of magnitude in almost all the cases. That is, MEGA can effectively approach the global optimum without being trapped in many local optima.
Because of the simplification property of the mutual-evaluation approach, the MEGA is suitable to deal with a wide spectrum of real world problems. In this paper, the proposed MEGA is also realized to deal with multipath routing problems in a multicommodity network, where more than one routes will be connected for each origin-destination (OD) pair. The goal is to find an optimal flow-allocation strategy to minimize the total transmission cost and satisfy all qualityof-service (QoS) requirements at the same time. The performance of the proposed MEGA is experimented by optimizing several multipath routing problems. The experiment results show that the MEGA outperforms the TGA dramatically. Furthermore, the search ability of the MEGA is robust in obtaining consistent results.
The rest of the paper is organized as follows. Section 2 describes the proposed mutual-evaluation approach. The main operations of the MEGA are introduced in Section 3.

Journal of Applied Mathematics
In Section 4, the experimental studies on 12 numerical optimization problems are described. The results are compared with eight existing optimization algorithms. In Section 5, the MEGA is realized to deal with a real-world problem: multipath routing problem. A novel representation of chromosomes for routing problems is also proposed. Finally, conclusions and contributions are offered in Section 6.

Mutual Evaluation Approach
To overcome the deficiencies of the TGA, modified GAs should keep evolutional population as diverse as possible to improve algorithms' exploration capability for discovering new solution area. In regard to evaluation approaches, existing GAs can be classified into three categories: (1) chromosome-oriented, (2) gene-oriented, and (3) schemaoriented.

Chromosome-Oriented Approaches.
Chromosome-oriented GAs concentrate on the chromosome fitness in the evolutionary process. The TGA adopts chromosome fitness to evaluate the quality of whole individual chromosome [4,5]. A number of researches recently focus on incorporating mathematical methods with the chromosome-oriented GAs to alleviate the deficiency of premature convergence [4]. For example, Yao et al. proposed fast evolutionary programming (FEP) with Cauchy mutation to solve the premature convergence deficiency [6]. To improve slow finishing deficiency, Tu and Lu developed a stochastic genetic algorithm (StGA) with mathematical methods, which not only can find global optima, but also can reduce the computational effort [7]. By merging niche techniques and Nelder-Mead's simplex method, Wei and Zhao proposed the niche hybrid genetic algorithm (NHGA) to alleviate premature convergence and weak exploitation deficiencies of TGA [3].

Gene-Oriented Approaches.
In gene-oriented GAs, the compact genetic algorithm (cGA) is a common representative [8]. The cGA represents the population as a probability vector over the set of individuals to mimic the order-one behavior of TGA. The cGA manipulates the gene distribution and essentially evolves each gene individually [8]. To enhance cGA's performance, Ahn and Ramakrishna proposed a strong elitism version of cGA [9] and then Rimeharoen et al. introduced a moving average technique to update the probability vector [10]. Although the cGA reduces memory requirement and offers many advantages, its limitation is the assumption of the independency between individual genes.

Schema-Oriented
Approaches. The schema-oriented GAs explore the exact schemata by borrowing from the schema theorem proposed by Holland [11]. A schema is a pattern within a chromosome defined by fixing the values of specific chromosome loci. The increase of effective schemata enables the efficient search within a solution space and guides the evolution of the population in approaching the global optimal solution [12]. Yen and Shyu proposed a statistical gene evaluation method that uses simple statistical quantities to investigate the individual gene influence which suggests better choices for a gene evolution [13]. Kubota et al. proposed the virus-evolutionary genetic algorithm (VEGA) [14] that simulates coevolution of a virus population and a host population. VEGA applies horizontal propagation and vertical inheritance in a population with virus infection operators and genetic operators [14]. In this paper, the proposed MEGA is a new schema-oriented GA in which an innovative mutual-evaluation approach is used to achieve highly efficient evolution with necessary robustness. The MEGA does not only evaluate a chromosome by its fitness but also analyze genome's merit to improve the population's quality.

The Proposed Mutual-Evaluation Approach.
According to the biological concept of the genetic engineering, a gene splicing is a process that manipulates genes outside the traditional random reproductive process. The proposed MEGA integrates a mutual-evaluation approach in a novel therapy crossover operation to produce offspring by introducing isolation, manipulation, and reintroduction of gene splicing techniques to improve the chromosome's fitness. The algorithm randomly selects two parents from a mating pool of generation . Let the parent with superior fitness be named as the good parent ⃗ good ( ) and the other one is the bad parent ⃗ bad ( ). The MEGA generates a therapy mask to indicate which gene loci in a chromosome are chosen for crossover points. For each bit in a chromosome, we uniformly generate a random number in interval [0, 1] and compare the number with a predefined therapy rate . If the random number is less than , its mask bit of the corresponding locus is set to value 1, which means that this gene locus belongs to the therapy genome. Otherwise, the mask bit is 0 that means the gene in this locus will not change during crossover operation.
Step 1. According to Darwin's evolution theory, two crossover parents are combined to produce new offspring in the hope that the fitness of next generation may improve gradually. In this paper, the therapy crossover wants to preserve parents' advantage and enhance population's diversity at the same time. Thus, offspring inherits the majority of parents' properties from the good parent than the bad one. Each mating parent has a different therapy rate according to its fitness; for example, good parent ⃗ good ( ) has a lower therapy rate (e.g., good = 0.45 in this paper) than that of the bad one ⃗ bad ( ) (e.g., bad = 0.9). On average, 45% of genes in the good parent will be merged with those genes in the bad parent. Example 1. A minimization function ( ) = ⃗ × ⃗ is adopted here to illustrate the rationale of the therapy crossover. Let two parents with five genes at generation be ⃗ , where ⃗ good ( ) is the good parent with better fitness ( ⃗ good ( )) = 5 and ⃗ bad ( ) is the bad one with worst fitness ( ⃗ bad ( )) = 20. Two therapy masks for ⃗ good ( ) and ⃗ bad ( ) are randomly generated by comparing five random numbers with therapy rates good = 0.45 and bad = 0.9, respectively. Without loss of generality, we assume that these two masks are ⃗ good = [0 1 0 1 0] and ⃗ bad = [1 1 1 1 0] . That is, the second and fourth genes of ⃗ good ( ) should be merged, and most of the genes in ⃗ bad ( ) should be merged except for the fifth one.
Step 2. For comparison purpose, two auxiliary chromosomes ⃗ good and ⃗ bad are generated for ⃗ good ( ) and ⃗ bad ( ), respectively. The auxiliary chromosome clones genes from the corresponding parent and replaces the selected therapy loci with those genes in the other parent. Thus, in (1), ⃗ good copies all genes from ⃗ good ( ) and then replaces the genes (its locus marked by ⃗ good ) by those genes in ⃗ bad ( ). On the other hand, (2) describes the value of each gene in ⃗ according to its therapy mask ⃗ bad : Notation (¬) denotes logic negation (usually expressed by "NOT") which operates on one Boolean value and returns its complement as result.
Example 2. Thus, (1) guides ⃗ good to copy the majority of genes (i.e., 1st, 3rd, and 5th) from ⃗ good ( ) and the minority of genes (i.e., 2nd and 4th) from ⃗ bad ( ), that is, ⃗ Step 3. Because it is difficult to determine the merit of a set of genes in a chromosome, a simple but effective method proposed in this paper is the mutual-evaluation approach, which measures the merit of two selected genomes by comparing the fitness changes before and after interchanging the genome with the other mating chromosome. Comparing the fitness change before and after the genome replacement (i.e., ( ⃗ good ( )) versus ( ⃗ good ) and ( ⃗ bad ( )) versus ( ⃗ bad )) can realize the substitution effect and can be used to represent the relative merit of these genomes.

Mutual-Evaluation Genetic Algorithm (MEGA)
The main operations of the MEGA are initialization, mutual evaluation, selection, crossover, mutation, and replacement. We depict the flowchart of the MEGA in Figure 1 and describe their functionality in this section.

Encoding and Initialization.
For illustration, the following minimization problem with fixed boundaries is considered:

Selection Operation.
Fitness of each chromosome represents the objective function value of this solution, denoted as for the th chromosome. The MEGA employs a traditional roulette selection method as a discriminator of the solution quality that uses chromosome's fitness to create a selective pressure towards global optimal solution. Those chromosomes with higher fitness should have a greater selection chance, thus creating a selective pressure towards high-fitness solutions [1].
3.2.1. Therapy Crossover. The proposed MEGA incorporates the mutual-evaluation approach with a therapy crossover to enhance the exploitation ability and speed up the convergence rate. According to the determined relative merit of genome, the parents linearly combine their genomes to generate a new genome for their offspring like (4) for ⃗ good ( + 1) and (5) for ⃗ bad ( + 1), respectively. If a mask bit in its therapy mask is 0 (e.g., 1 = 0), the gene at the locus of the good parent does not change (i.e., 1 = 1 ). Otherwise, if a mask bit is 1 (e.g., 2 = 1), the locus in the crossover child (e.g., ⃗ good ( + 1)) inherits from both genetic materials from two parents (i.e., 1 × coef + 1 × (1 − coef)). The principle of this linear combination is to retain the favorable schemata in the evolution process. Therefore, each gene of crossover child can be reproduced by the following equations: Coefficient coef = 0.8+(rand/2) is a random value in interval [0.8, 1.3] to reduce the deficiency of premature convergence. The rand function returns a uniformly distributed pseudorandom number between 0 and 1.

Example 4.
In this example, we assume the random coeffi- The child's fitness ( ⃗ bad ( + 1)) = 6.56 is better than ( ⃗ bad ( )) = 20. Therefore, both children have better fitness values than that of their parents.
The exclusive features of the therapy crossover include that (1) the merit of each genome is evaluated individually; and (2) the gene merit facilitates the MEGA to perform an efficient search by adaptively shifting emphasis on significant genome without explicit functional analysis [15]. That is, the therapy crossover can avoid frequently throwing away potential schemata in inferior chromosomes and inherit the genetic advantages of superior chromosomes without loss of genetic diversity.

Partial-Gaussian Mutation.
In this paper, a partial-Gaussian mutation used in the MEGA can increase population diversity to enhance its exploration ability. The original Gaussian mutation proposed by Hinterding in 1995 can converge to near-optimal solutions of some multimodal optimization problems [16]. Our partial-Gaussian mutation concentrates on exploiting potential optimal areas and speeds up the convergent effect. At the beginning of evolution, a high mutation rate is assigned to sample the search space extensively. And then, the mutation rate decreases with time to fine tune solutions and concentrates on exploiting potential optimal areas. Equation (6) calculates the mutation rate as where maximal generation is 3000. A random number in interval [0, 1] is generated for each gene and then compared with the mutation rate . If the mutation rate is greater than or equal to the random number, this gene value will be flipped by adding a unit Gaussian distributed random value to the chosen gene. Otherwise, no mutation occurs at this gene. This mutation operation can only be used for integer and float genes.

Reproduction Operation.
The MEGA adopts a replacement-with-elitism method to monotonously enhance the solution quality. A number of the elite parents can survive into next generation for preventing good solutions from being lost through a nondeterministic selection operation. Successive population consists of three evolutionary sources: (1) the elite 10% chromosomes can survive to next generation; (2) 80% offsprings are produced by the crossover operation; and, (3) the rest 10% chromosomes are produced by the mutation operation.

Numerical Test Functions.
Numerical experiments are conducted to demonstrate the robustness and reliability of the proposed MEGA. This work selects 12 well-known benchmark functions, which cover broad range functionality characteristics with two categories: unimodal functions (Functions 1 -6 ) and multimodal functions (Function 7 -12 ).  ( 5 ), and the noisy quadratic function ( 6 ). Unimodal functions are relatively easy to solve but the difficulty increases as the problem dimension goes high. The high-dimension multimodal functions are a generalized Schwefel's function 7 ( 7 ), the Rastrigin's function 6 ( 8 ),     a modified Ackley's Path Function 10 ( 9 ), the Griewank's function 8 ( 10 ), a generalized Penalized Function 1 ( 11 ), and a generalized Penalized Function 2 ( 12 ). Multimodal functions represent the most difficult class of problems, which possess many local optima and could trap an algorithm into one of its local optimal solutions.

Algorithm Implementation and Parameter
Settings. In all cases, the population size is 150, in which the number of elite individuals is 15; the therapy crossover produces 120 individuals, and the mutation produces 15 ones. The therapy rates for good parent and bad parent are good = 0.45 and bad = (1 − good ) = 0.55, respectively. For each test function, 50 independent trials with different seeds are performed using the MATLAB environment.
For complexity analysis, the mean number of function evaluations serves as a measure of required computational effort for an algorithm. Different from the crossover in traditional GAs, the exclusive feature of MEGA is the usage of the mutual-evaluation approach for genome therapy. Because  Table 2 summarizes the experimental results of the MEGA in 50 trials, which include (1) the global optimum ( min ), (2) the mean best function value, The first thing we can observe from Table 2 is that the obtained results are equal or really close to the "known" optimal solutions. Particulary, this paper uses a computational precision of 60 digits after point. Thus, the results "0" on 8 and 10 in Table 2 mean that they are less than 10 −60 . The standard deviation with respect to the functions 8 and 10 is equal to zero; that is, the results of all 50 runs reach the optimum. All the obtained results approach the "known" optimal values with small differences. Secondly, the small standard deviations for all test functions also indicate that the MEGA consistently converges to the near-optimal solutions in all 50 trails. Finally, all of the mean numbers of function evaluations are relatively small. Therefore, the proposed MEGA can address a variety of numerical optimization functions effectively. To further analyze the solution capability of the proposed MEGA, the following sections describe the comparisons between the MEGA and two groups of optimization algorithms for the 12 benchmark functions.

Comparison with Other Evolutionary Algorithms.
The performance of the MEGA is compared with four state-ofthe-art EAs: particle swarm optimization (PSO) [17], evolutionary optimization (EO) [17], conventional evolutionary programming (CEP) [6], and FEP [6]. The comparison of the experimental results for 12 test functions is shown in Table 3.
For all the unimodal functions ( 1 -6 ) in Table 3, we can observe that the MEGA can achieve dramatically the highest accuracy than others, while other four algorithms experience premature convergence on functions 3 , 4 , and 5 . For all the multimodal functions ( 7 -10 ), the results shown in Table 3 clearly indicate that the MEGA can identify the actual optima of these functions with the highest accuracy. The MEGA can achieve better solution accuracy than other four EAs for all 12 benchmark functions.
The computational efforts required for the algorithms are measured by their mean numbers of function evaluations and depicted in Table 4. Obviously, the comparison demonstrates that the MEGA outperformed all the four algorithms for all the 12 functions with respect to the convergent ability. Therefore, the comparison indicates that the MEGA is both efficient and effective in solving the unimodal and multimodal benchmark functions.

Comparison with Other Genetic Algorithms.
The performance of the MEGA is compared with those of four well-known GAs: cluster-based adaptive mutation genetic algorithm (CMGA) [18], orthogonal genetic algorithm with quantization (OGA/Q) [19], hybrid taguchi genetic algorithm (HTGA) [20], and StGA [7]. The OGA/Q is a quantizedversion of the OGA, which incorporates with the Taguchi method to minimize the effect of chromosome variation without eliminating the population diversity [21]. The HTGA enhanced the Taguchi method as a new operation to adapt a dynamically extended precision method from a lowprecision solution space to a high-precision one [20]. The above algorithms have been executed to solve the test functions and the results were reported in the literature. We will use these existing results for a direct comparison in Tables 5  and 6.
As the termination criteria used in these four algorithms are different, to make a fair comparison basis, we let the solution qualities obtained by our MEGA be slightly better than those of the four algorithms (in Table 5), and then, compared the mean computational effort at the given accuracy (in Table 6). For the unimodal functions ( 1 -6 ), the convergence rate of an algorithm is a more important issue   9 . This finding implies that the proposed MEGA is robust and effective in solving the multimodal functions and can perform as good as the four state-of-the-art GAs.

Multipath Routing Problems.
Multipath traffic engineering becomes more attractive for ubiquitous networks to satisfy the required QoS of mobile network applications [22]. For a given network topology with available link capacities, it is required to determine the optimal distribution of traffic requests on multipath routing subject to the constraints imposed by QoS specifications. The service network is modeled as a connected weighted, directed graph = ( , ), where the = {V 1 , V 2 , . . . , V } is the vertex set of , and the = { 1 , 2 , . . . , } is a finite set of edges. We consider a link-state routing environment, where each node has knowledge of its neighbor links and each source node knows (1) traffic load from node to node ( ̸ = ), (2) QoS requirements/constraints (e.g., delay constraint), (3) residual bandwidth of each link, (4) link cost for the traffic, and (5) time delay for the traffic to pass through each link at a certain time period.
The link = ( , ) ∈ connects a source node ∈ to a destination node ∈ with positive cost , capacity , and a delay function (i.e., : R + → R + ). The cost function represents the traffic-related delay of each edge and limit to positive cost. The delay requirement specifies the upper bound of delay tolerance, denoted as Δ .
In multipath routing problems, messages are routed for a set of OD pairs , where denotes the traffic demand volume (or bandwidth) for OD pair . The set of all permissible paths for the th OD pair is denoted as . Let indicator function be a nonnegative continuous variable denoting the traffic flow allocated on path of request . Indicator function is 1 if link belongs to path for request , and 0 otherwise; that is, if ∈ otherwise, where = 1, 2, . . . , . The aggregate flow on link is denoted as , which must satisfy the capacity feasibility constraint; that is, = ∑ ∈ ∑ ∈ ≤ , for all ∈ . The multipath routing problem is to find multiple paths to transmit the traffic demand for all OD pairs such that the total routing cost is minimal and the QoS constraints should be satisfied. The multipath routing can be modeled as a combinatorial linear programming problem in the following: objective function, The objective function (in (7)) is to minimize the total routing cost in the multipath routing problem. The first constraint (in (8)) calculates the aggregate traffic on link . The second constraint assures the capacity constraint on each link by (9). The third constraint (in (10)) ensures the total delay of each OD pair to satisfy a prespecified path-delay bound Δ . The constraint in (11) enforces that traffic demand of each OD pair should be satisfied. Equation (12) lets an auxiliary binary variable be 1 if the traffic along path is larger than zero, and 0 otherwise. Equation (13) sums up the used routes and assigns the number to the th gene .

Algorithm Implementation for
Multipath Routing Problems 5.2.1. Encoding. The encoding method is the most important steps towards solving real world problems using EAs. In this paper, the encoding method maps all multipath OD pairs into a chromosome based on route aspect. Since there are so many candidate paths between two nodes in the network graph, traditional GAs may consume considerable computational effort in searching infeasible solutions because genetic operations do not always preserve feasibility. Therefore, to reduce the search space, this work uses the shortest path routing algorithm to precalculate the first shortest paths and record in a routing table.
The MEGA maintains a population of chromosomes to optimize a given objective function for the multipath routing problem. Each chromosome can be represented by a twodimensional array of integers. The first-dimension genes ⃗ = [ 1 2 ⋅ ⋅ ⋅ ] represent the number of used routes for each OD pair. The second-dimension genes record the route number of each used path for realizing the demand of each OD pair. Thus, if the gene value of the th gene is = , the vector [ 1 , 2 , . . . , ] represents route numbers of subflows for OD pair . Gene is an integer in interval [1, ] to represent a route number for the OD pair in its routing table. We use an example of network topology to illustrate the relationship between a chromosome, 2D genes, and routing tables in Figure 2. Figure 2(a) depicts an example of a network graph, link parameters, and two multipath routings. Parameters along links are triple (cost, delay, and bandwidth). The first OD pair, that is, 1 = (1, 3), has two routing paths 12 = {(1, 4), (4, 5), (5, 3)} and 11 = {(1, 2), (2, 3)} to transmit traffic from node 1 to node 3. Thus, the first-dimension gene 1 = 2. 1-4-5-3 · · · · · · R 1-4-7-5-3 Routing table for OD pair (6,8) Route no. Route path  The route numbers in its routing table (in Figure 2(c)) corresponding to paths 12 and 11 are 2 and 1; thus, the second-dimension genes for 1 (in Figure 2(b)) record 11 = 2 and 12 = 1, respectively. In Figure 2(a), the second OD pair, that is, 2 = (6, 8) routes its traffic along three subflows along paths 21 , 23 , and 24 . Therefore, the first-dimension gene is 2 = 3 and the corresponding second-dimension genes are 21 = 1, 22 = 3, and 23 = 4, respectively.

Therapy Crossover for Multipath Routing.
Each time the selection operation chooses two crossover parents from the population. The proposed mutual-evaluation approach calculates the merit of two selected genomes by comparing the changes in the chromosome fitness before and after interchanging the first-and second-dimension genomes with the other mating chromosome. The first-dimension genes can linearly combine with those in the other chromosome by using the proposed therapy crossover (in Section 3.2.1 (4) and (5)). The obtained results should be transformed into integer type, and therefore, the second-dimension genes should be modified. If the obtained result of the first-dimension gene is larger than before, we randomly select a suitable number of genes from the other parent to add into the seconddimension genes in this chromosome. If the result equals to the original value, we randomly select a small number of genes from the other parent to replace the original genes. Otherwise, a suitable number of genes should be randomly selected to remove from this chromosome.

Mutation for Multipath
Routing. Mutation operation performs on an individual chromosome to flip one or more genes with a small probability (typically 0.001) and ensures that no point in the search space has a zero probability of being searched. According to a mutation probability, the mutation randomly selects a subset of genes and chooses new paths from its routing table. Thus, the route numbers of these new paths replace the original values of selected genes. The resulting chromosome is a new multipath routing plan that can increase population diversity.

Test Platform and Performance Metrics.
In this paper, we use the well-known network generation tool [23] to create an asynchronous network based on the Waxman's techniques [24]. The network in Figure 3 illustrates a random generated  graph with 14 nodes and the average degree of each node is four. Figure 3 only shows the cost/delay information along one direction link (from a smaller-ID node to a larger-ID one) to reduce the complexity of the graph representation. All links are assumed to have 1.5 Mbps of bandwidth capacity. For each test scenario, the OD pairs are randomly generated five times for each scenario to decrease the selection bias. Two kinds of transmission rates are assigned for the OD pairs: 0.25 Mbps (lowrate) and 0.5 Mbps (highrate). The performance of the proposed MEGA is evaluated based on the following ways.
(1) The total routing cost: This cost reflects the algorithm's ability to construct multiple paths for all OD pairs by using low-cost and lightly utilized links.
(2) The maximum end-to-end delay for OD pairs: It indicates the algorithm's ability to satisfy the delay bound imposed by the service level agreement of applications.
An algorithm's effectiveness in allocating network resources can be judged by monitoring how frequently that algorithm fails to construct a set of acceptable OD pairs. There are two kinds of failure. One is that the created OD pair does not satisfy its delay bound. The other one is that the algorithm cannot find unsaturated links to create a path for OD pairs.

Computational Experiments for Multipath
Routing Problems. The performance of the proposed MEGA is compared with the TGA for two kinds of multipath routing scenarios with respect to different transmission rates and different number of total OD pairs in the network. In the first test scenario, the experiments are performed for 100 independent runs with 10 generations in the MATLAB environment. The evolution curves of total routing costs with respect to lowrate and highrate are compared in Figures 4(a) and 4(b), respectively. The evolution curves in Figure 4(a) show that the efficiency of the MEGA is better than that of the TGA in both two transmission rates. Particulary, even though the transmit rate increases two times than the low-rate case, the growing ratio of the total routing cost in Figure 4(b) is still less than that obtained by the TGA. The second scenario simulates the stress test to measure the robustness of these two algorithms. We increase the number of OD pairs from 3 to 10 and conduct 100 independent trials for each test case. The OD pairs are randomly generated five times for each scenario to simulate the great diversity of OD-pair selection for performance evaluation.
Tables 7 and 8 are the "mean" and "standard deviation" experimental results obtained by the proposed MEGA and TGA for five cases in the low-rate scenario (with 0.25 Mbps transmission rate). The mean experimental results obtained by the proposed MEGA are better than the results of the TGA on all the test cases. That is, the MEGA can find better routing paths to serve all multipath transmission requirements than TGA (in Table 7). For the standard deviations in 100 independent runs, the proposed algorithm achieved superior results compared with the TGA in all the test cases (in Table 8). This finding implies that the proposed algorithm is robust in solving multipath routing problems and can perform better than the TGA.
In the high-rate scenario, the test network cannot afford the data flows if the number of OD pairs is larger than 6. Thus, the following experiments only increase the number of OD pairs from 3 to 6 and the transmission rate is 0.5 Mbps for the high-rate cases. We simulate 100 independent trials for each test case and depict the "mean" and "standard deviation" of the experimental results obtained by the MEGA and TGA in Tables 9 and 10, respectively. Experimental results indicate that the proposed MEGA outperforms the TGA with respect to the "mean" and "standard deviation" of routing costs for all high-rate scenarios. Furthermore, the search ability of the MEGA is robust in obtaining consistent results and performs better than the TGA.

Conclusions and Future Works
To the best of our knowledge, the proposed MEGA is the first mutual-evaluation approach, which calculates each genome merit by interchange-compare-replace method. The genome evaluation facilitates the MEGA to perform an efficient search by dynamically shifting emphasis to significant genomes in the feasible space without abdicating any portion of the candidate schemata. The therapy crossover is also proposed to preserve better-performance schema patterns. Simpler than other modified approaches, the proposed MEGA can preserve high quality genomes during evolution period without using extra analyzing techniques.
The performance of the proposed algorithm was measured using 12 benchmark functions. The performance was compared with four existing EAs and four well-known GAs. The experiment results show that the MEGA is able to find near-optimal solutions, even though other algorithms experience difficulties in approaching the global optima on some functions. The behavior of the algorithm is also consistent as indicated by a small standard deviation among the 50 trials for each test function. Furthermore, the MEGA can increase the accuracy by several orders of magnitude than other algorithms in almost all test functions. That is, the MEGA can outperform the existing global optimization algorithms with a dramatic improvement in terms of effectiveness. Furthermore, because of the simplification property of the mutual-evaluation approach, this study slightly modified the encoding method of the MEGA to solve the multipath routing problems in multicommodity networks. The bandwidth-delay constraints are introduced to enforce the service quality of multimedia applications. The experimental results show that the MEGA not only can solve QoS constrained multipath routing problems, but also can outperform the TGA. That is, the proposed MEGA not only can reduce the effort of explicit function analysis but also can deal with a wide spectrum of real world problems.
The contributions of the paper are as follows.
(1) We develop a novel mutual-evaluation approach which incorporates with a GA that has never been jointly considered in the literature. (2) We introduce a novel therapy crossover, which not only can evolve superior genomes but also can achieve global optima. (3) We introduce a novel chromosome representation for the MEGA to address multipath routing problems in a multicommodity network. (4) Experimental results show that the proposed MEGA can achieve significant performance gain over several well-known algorithms under the considered scenarios. The proposed MEGA has high exploration and exploitation abilities as a robust, statistically sound, and quickly convergent algorithm.
We have observed that there are many researches for routing problems. In the future, we will further compare with more state-of-the-art methods on QoS multipath routing. Particulary, the simplicity property of the MEGA can help to route dynamic services across heterogeneous environments. Therefore, developing a distributed MEGA to enhance its scalability for highly dynamic environments is also our future work.