Decision variable contribution based adaptive mechanism for evolutionary multi-objective cloud workflow scheduling

Workflow scheduling is vital to simultaneously minimize execution cost and makespan for cloud platforms since data dependencies among large-scale workflow tasks and cloud workflow scheduling problem involve large-scale interactive decision variables. So far, the cooperative coevolution approach poses competitive superiority in resolving large-scale problems by transforming the original problems into a series of small-scale subproblems. However, the static transformation mechanisms cannot separate interactive decision variables, whereas the random transformation mechanisms encounter low efficiency. To tackle these issues, this paper suggests a decision-variable-contribution-based adaptive evolutionary cloud workflow scheduling approach (VCAES for short). To be specific, the VCAES includes a new estimation method to quantify the contribution of each decision variable to the population advancement in terms of both convergence and diversity, and dynamically classifies the decision variables according to their contributions during the previous iterations. Moreover, the VCAES includes a mechanism to adaptively allocate evolution opportunities to each constructed group of decision variables. Thus, the decision variables with a strong impact on population advancement are assigned more evolution opportunities to accelerate population to approximate the Pareto-optimal fronts. To verify the effectiveness of the proposed VCAES, we carry out extensive numerical experiments on real-world workflows and cloud platforms to compare it with four representative algorithms. The numerical results demonstrate the superiority of the VCAES in resolving cloud workflow scheduling problems.


Introduction
Workflows have been commonly used to describe data processing applications from diversified fields, such as the Internet of Things and bio-informatics [1][2][3].These workflows often comprise of large-scale data-dependent tasks, which are computation and data intensive.Then, executing various workflow applications calls for powerful highperformance infrastructures.With substantial advantages, such as economies of scale, on demand supply of resources, high elasticity, and reliability, cloud computing is attracting more and more enterprises or individuals to deploy their big data processing workflows [4,5].
Workflow scheduling in cloud computing is a key technology for achieving the reduction of both execution cost and makespan to gain more profits for cloud providers and ensure the quality of service for cloud consumers [6].Workflow scheduling problem involves determining the mappings from tasks to resources and the task order on each resource, and is a classic NP-complete [7,8].Also, the execution cost and makespan of workflow scheduling are two conflicting optimization objectives [9].So far, multi-objective evolutionary algorithms have become popular to search a set of compromise solutions within an acceptable time [10][11][12].To solve the workflow scheduling problem in cloud computing, some studies design new evolution and selection operators to improve the classical multi-objective evolutionary algorithms.
Over the past decade or so, designing efficient evolution operators to reproduce new solutions for multi-objective workflow scheduling problem has attracted considerable research interest [13].First of all, the popular list-based workflow scheduling methods were embedded into the multiobjective evolutionary optimization framework as evolution operators [14][15][16].Secondly, bio-inspired evolution operators, such as artificial neural network [17][18][19], ant colony optimization [20], firefly algorithm [21], particle swarm optimization [22][23][24], and grey wolf optimization [25], were modified as evolution operators to solve the multi-objective workflow scheduling problem.Thirdly, integrating heuristic rules and bio-inspired optimization techniques to reproduce offspring populations has become a popular technological path.For instance, Choudhary et al. [26] combined the gravitational search method and the list-based workflow scheduling method to solve bi-objective workflow scheduling problem in cloud computing.Hosseini et al. [27] merged the simulated annealing and a task duplication strategy to optimize the makespan and execution cost of workflows.Mohammadzadeh et al. [28] integrated the antlion and grasshopper optimization algorithms to balance throughput, makespan, cost, and energy consumption of executing workflows in cloud platforms.Zhang et al. [29] enhanced the list-based workflow scheduling method with a local search mechanism to balance the makespan and energy consumption of workflow execution.
At the same time, some studies went into designing selection operators to balance multiple conflicting objectives of workflow execution in cloud computing.For example, Zhou et al. [30] merged a fuzzy-dominance-based environmental selection and a list-based workflow scheduling method to minimize execution cost and makespan of workflows in cloud computing.Kumar et al. [31] integrated the entropy weight mechanism into a multi-criteria decision-making framework to balance makespan, execution cost, reliability, and energy consumption.Ye et al. [32] improved a knee point driven evolutionary method to balance makespan, reliability, execution cost, and the mean durations of all workflow tasks.Pham et al. [33] focused on the volatility of spot cloud resources and improved the multi-objective evolutionary algorithm to make a trade-off between makespan and execution cost for workflows in cloud computing.
In the evolutionary optimization community, a multiobjective optimization problem is generally considered large-scale if it has at least one hundred decision variables [34].The multi-objective cloud workflow scheduling involves hundreds or even thousands of decision variables, and is a typical large-scale multi-objective optimization problem.However, the existing relevant studies evolve all decision variables as a whole and allocate evolution opportunities to each variable equally.This results in the low efficiency of these existing studies.
The recent research results in the evolutionary computation community demonstrate that cooperative coevolution [34,35] has become a crucial and effective way to solve large-scale multi-objective optimization problems.In cooperative coevolution approaches, all the decision variables are classified into multiple groups, and decision variables in different groups are evolved in a round-robin manner [36,37].These static classification techniques work well when problems' decision variables are fully or partially separable.However, this is not the case for multi-objective cloud workflow scheduling with nonseparable decision variables caused by data dependencies among tasks.
Besides, multi-objective cloud workflow scheduling poses an imbalance feature among decision variables regarding their contributions to optimization objectives.For instance, delaying the completion of a workflow task on the critical paths [38] often successively delays the completion of many tasks, including its successor tasks and other tasks being executed after the delayed tasks and their successor tasks.Whereas slightly delaying other tasks on the non-critical paths may not cause this chain reaction.The imbalance feature means that we should equip different decision variables with different evolution opportunities.This motivates us to design a decision variable contribution based adaptive mechanism to dynamically adjust the variable grouping and allocate evolution opportunities during the evolution process.Our main contributions in this paper are as follows.
• We define the contribution of a decision variable as the fitness improvement of the solution generated by perturbing this decision variable.Then, we try to dynamically measure the contribution of each decision variable and classify them according to their contributions.• We design an adaptive mechanism to dynamically allocate more evolution opportunities for the variable groups with more contributions to generate offspring solutions efficiently.• In the context of fifteen real-world workflows and the Amazon Elastic Compute Cloud, we compare the proposal with four state-of-the-art multi-objective cloud workflow scheduling algorithms.The results demonstrate the competitive performance of the proposal in simultaneously optimizing execution cost and makespan.This paper is organized as follows.In the second section formulates the multi-objective workflow scheduling problem.In the third section designs the proposed VCAES, followed by experimental verifications in the fourth section.In the final section concludes this paper.

Problem formulation
This section first describes the models of workflows and cloud resources, and then formulates the model for multiobjective workflow scheduling in cloud computing.

Workflow model
Without loss of generality, a workflow application is described by a Directed Acyclic Graph (DAG), whose vertices and directed edges represent workflow tasks and data dependencies, respectively.Formally, we construct the directed acyclic graph for a workflow application as Ψ = {T , D}, where T = {t 1 , t 2 , . . ., t n } denotes the vertex set corresponding to task set, D ⊆ T × T denotes the edge set corresponding to data dependencies among tasks.The existence of an edge d i, j ∈ D means that t j 's start demands t i 's output results.Generally, task t i is referred to as an direct predecessor of task t j , and t j is referred to as an direct successor of t i .Regarding a task t i , all its direct predecessors is expressed as a set P(t i ), and all its direct successors is expressed as a set S(t i ).
Figure 1 gives a visual example of a directed acyclic graph for a workflow with seven tasks, i.e., T = {t 1 , t 2 , . . ., t 7 }.An edge d 1,2 denotes the data dependency from t 1 to t 2 , meaning that t 2 's start have to wait for t 1 's output results.In Fig. 1, regarding task t 6 , the set of its direct predecessors is P(t 6 ) = {t 3 , t 4 }, and the set of its direct successors is S(t 6 ) = {t 7 }.

Cloud resource model
This paper targets the popular cloud paradigm, i.e., Infrastructure as a Service (IaaS).In this paradigm, cloud providers offer multiple types of cloud resources on demand [39,40].
The differences between different types of cloud resources mainly lie in their charging prices and performance configurations, such as number of CPU cores, memory size, and network bandwidth.Assuming that cloud platforms offer m types of resources, then we model them as Γ = {1, 2, . . ., m}, where τ ∈ Γ denotes the τ -th resource type.Regarding a type τ , we employ pr(τ ) and con(τ ) to represent its price and configurations.Then, a cloud resource of type τ is modeled as r τ k = {k, pr(τ ), con(τ )}, where k denotes the index of resource r τ k .Refer to well-known cloud providers (e.g., Amazon EC21 and Alibaba Cloud ECS2 ), this study follows the resource charging basis of pay-as-you-use.Under this rule, any consumer can rent any number of resources on demand and is charged according to the real usage time.In general, cloud resources are charged based on the number of billing periods, and the partial period will be rounded up to one more.In case that the period length is 60.0 min, the number of billing periods for 60.01 min is two.

Multi-objective scheduling cloud workflows
Since cloud resources are available on demand, we build a resource pool based on the maximum resource requirements for running a workflow.Assuming the maximum parallelism of the workflow is p, the resource pool includes p resources of each type.Then, we describe the resource pool as: R = r 1 1 , r 1 2 , . . ., r 1 p , r 2 p+1 , r 2 p+1 , . . ., r 2 2• p , . . ., r m m• p .The decision vector x = {x 1 , x 2 , . . ., x n } is used to represent the mappings from workflow tasks to cloud resources, where the value of decision variable x i is decoded as the index of the cloud resource mapped to the i-th task.It is worth noting that the value range of each decision variable is an integer from 1 to m • p.
Given a decision vector, assume that the task t i is mapped to resource r τ k .The start time st i,k of task t i refers to the maximum time of receiving all the input data and the available time of the mapping resource.
On resource r τ k , we assume the set of tasks being executed before task t i as follows: where I (t p ) indicates t p 's order number on resource r τ k .Then, the start time st i,k of task t i on cloud resource r τ k can be described as follows: where ft b,k indicates t b 's finish time on resource r τ k , ft p, * indicates the finish time of task t p , and dt p,i indicates the data transfer time from t p to t i .Before scheduling, task t i 's execution time et i,k on cloud resource r τ k can be predicted by the computation length of task t i and the CPU frequency of the mapped resource r τ k .The relationship among st i,k , et i,k , and ft i,k is described as follows: ( Given a decision vector, the set of all tasks mapped to cloud resource r τ k can be formulated as: With the task set T k , the start time ut k and end time nt k of renting resource r τ k can be computed as follows: ( Based on the above analysis, the first optimization objective, i.e., minimizing the execution cost, can be formulated as follows: where C indicates the length of a billing period for cloud resources.
The second optimization objective of this paper is to minimize the workflow's makespan, which corresponds to the maximum finish time of all the tasks in this workflow.The second optimization objective can be formulated as follows: Thus, the model for multi-objective workflow scheduling problem in cloud computing can be summarised as follows: Pareto-dominance has been widely employed to compare solutions in the multi-objective optimization field.
Pareto-dominance: Assuming x 1 and x 2 are two feasible solutions.x 1 is regarded to dominate x 2 (denoted as x 1 ≺ x 2 ) if and only if the two objectives of x 1 are not inferior to that of x 2 (i.e., f j (x 1 ) ≤ f j (x 2 ), ∀ j ∈ {1, 2}) and x 1 is better than x 2 on at least one objective (i.e., f j (x 1 ) < f j (x 2 ), ∃ j ∈ {1, 2}).
Pareto Set/Front: All the Pareto-optimal solutions are defined as Pareto-Set (PS) in the decision space and Pareto-Front (PF) in the objective space.

Algorithm design
Given a workflow scheduling solution, the importance of each workflow task varies greatly.For instance, adjusting the mapping from a critical task to a resource often successively affects the execution of many tasks, including its successors and other tasks being executed after these tasks and their successors.Whereas adjusting the mapping from a non-critical task to a resource may have no impact on the execution cost and makespan of the workflow.Also, the importance of each workflow task varies from solution to solution.Then, decision variables corresponding to different workflow tasks pose an imbalance feature to optimization objectives.To deal with the large-scale decision variables in cloud workflow scheduling, the VCAES incorporates a novel cooperative coevolution (CC) mechanism to dynamically measure the contributions of decision variables and adaptively allocate evolution opportunities for each group of decision variables based on their contributions.The proposed VCAES follows the framework of traditional multi-objective evolutionary optimization, including initialization, reproduction operator, and selection operator, as shown in Algorithm 1.

Algorithm 1 The overall framework of VCAES
Input: The problem in (8); population size N ; memory length l; group size of decision variables s Output: A final population P 1: P ← Randomly generate a population 2: M ← O l×n 3: V ← Generate a set of reference vectors 4: g ← Initialize the number of generations for CC 5: while does not reach termination condition do 6: [P, K ] ← AdaptiveCoEvolution(P, M, V , g) 7: for j = 1 → l − 1 do 8: end for 11: end for 12: end for 15: Q ← Generate a new population 16: P ← EnvironmentalSelection(P Q, N ) 17: end while As illustrated in Algorithm 1, the inputs of the proposed VCAES are the multi-objective cloud workflow scheduling problem, the population size, memory length for recording the variable contributions, the number of decision variables in one group.Once the VCAES reaches the termination condition, it will output an up-to-date population.
In the initialization stage, one population is generated randomly (Line 1).Next, an l × n matrix M is initialized to collect the contribution of each variable over the past l iterations (Line 2).The element in row j and column i represents the contribution of the i-th variable in the previous j-th iteration.Also, a set of uniformly distributed reference vectors are initialized to assist in calculating variable contributions (Line 3).In addition, the number of iterations for cooperative co-evolution during each generation is initialized (Line 4).These iterations will be allocated to each group of variables in proportion to the overall contribution of variables in the corresponding group.
After the initialization stage, the VCAES enters the main loop.During each generation, the proposed adaptive cooperative coevolution mechanism is triggered to distribute the decision variables into many groups and allocate evolution opportunities to each group according to variable contributions (Line 6).Next, the memory matrix of variable contributions is updated (Lines 7-14).It is worth noting that the decision variables in the cloud workflow scheduling are related to each other.The VCAES generates a new population by evolving all variables in each generation (Line 15).After that, the non-dominated sorting and elitist-preserving method in NSGA-II [41] is employed to select an offspring population P from the combined population P Q (Line 16).
Before introducing the proposed adaptive cooperative coevolution mechanism, we define and illustrate the variable contributions.
Suppose Q is a population that generated by evolving the decision variables in group G(i) and fixing other decision variables, the contribution of each variable in G(i) is defined as: where FI(q) denotes the fitness improvement of solution q.
For a set of reference vector V , the one associated with solution q is defined as v * = arg v∈V min f (q), v , where f (q), v represents the acute angle between two vectors.Suppose p is the associated solution of the reference vector v * , then the fitness improvement of solution q can be calculated as follows: = f q =(0.5, 0.5) v Fig. 2 An example workflow of seven tasks where Fit(q, v) denotes the fitness value of solution q with respect to the reference vector v, which can be calculated as (11) where f = f (q)− z denotes the translated objective vector, and z denotes the ideal point.
An intuitive example on calculating contribution is given in Fig. 2. Suppose previous solution p and the new solution q are associated to the reference vector v = (0.5, 0.5), and their objective vectors are f ( p) = (1.3,1.0) and f (q) = (0.7, 0.8), respectively.The ideal point is supposed to be z = (0, 0).Based on the definition in (11), we can obtain the fitness of these two solutions as Fit( p, v) = 1.8385 and Fit(q, v) = 1.1314.Then, the fitness improvement of solution q is FI(q) = max{0, 1.8385 − 1.1314} = 0.7071.If the solution q is generated by evolving variables {1, 3, 5}, then the contribution of these variables are defined as Algorithm 2 gives the pseudo-code of the proposed adaptive cooperative coevolution mechanism, which dynamically distributes the variables with higher contributions into the same groups and assigns more evolution opportunities to accelerate the population convergence.
As illustrated in Algorithm 2, the inputs of the function AdaptiveCoEvolution() are the current population, the contribution of each variable during the previous l generations, the set of reference vectors, the population size, and the number of generations for cooperative coevolution search.Then, the outputs of this function are the updated population and the contribution of each variable.
At first, function AdaptiveCoEvolution() calculates the total contribution of each variable during the previous generations (Lines 1-4), where the notation H (i) represents the total contribution of the i-th variable during the previous l generations.Then 10: Q ← Generate an offspring population by evolving variables in G( j) 12: P ← EnvironmentalSelection(P Q, N ) 13: end for 14: c ← Calculate the contributions of variables in G( j) as ( 9) 15: for i = 1 → n do 16: end for 18: end for order according to their total contributions (Line 5) and distributes the variables into a series of groups (Line 6).In this way, the variables with similar contributions are distributed into the same group, which is helpful to distribute variables with larger contributions into the same groups to gain more evolution opportunities.
After that, function AdaptiveCoEvolution() successively evolves each group of variables, and measure their contributions.For a group of variables, the number of evolutionary iterations is calculated based on the sum of their contributions (Line 9).During each iteration, this function performs the reproduction operator on the corresponding variables to generate a new population (Line 11), and performs the selection operator to select the offspring population (Line 12).After evolving a group of variables, this function measures their contributions (Line 14).Also, it updates the contribution of the corresponding variables (Lines 15-17), where notation K (i) represents the contribution of the i-th variable.
Regarding Function AdaptiveCoEvolution(), it costs O(n •l) to calculate the contribution of each decision variable (Lines 1-4, Algorithm 2).The time complexity of sorting decision variables is O(n log n) (Line 5, Algorithm 2).This function takes O(N • n) to reproduce an offspring population (Line 11, Algorithm 2).According to the analysis [41], the time complexity of the environmental selection is O(N 2 ) (Line 12, Algorithm 2).Then, it takes O(n to adaptively each group of decision variables (Lines 8-18, Algorithm 2).Then, the time complexity of Function AdaptiveCoEvolution() is Regarding the algorithm VCAES, the time complexity of updating the memory matrix is O(n • l) (Lines 7-11, Algorithm 1).The time complexities of reproducing an offspring population and selecting elitist solutions are O(N • n) and O(N 2 ) (Lines 15-16, Algorithm 1), respectively.Thus, the time complexity of the VCAES during one generation is

Numerical experiments
To investigate the performance of the proposed VCAES, we compare it with four representative multi-objective cloud workflow scheduling algorithms: MOELS [42], EMS-C [9], WOF [43], LSMOF [44].The MOELS and EMS-C follow the framework of NSGA-II [41] and incorporate new reproduction operators to generate offspring populations by evolving all the variables.The WOF and LSMOF are representative large-scale multi-objective evolutionary algorithms based on problem transformation.

Experimental setting
The eight types of workflows from different application domains, i.e., Inspiral (Gravitational physics), CyberShake (Earthquake), EpiGenomics (Biology), Montage (Astronomy), Sipht (Bioinformatics), BLAST (bioinformatics), Cycles (agroecosystems), and Seismology (seismology), have been widely used in evaluating cloud workflow scheduling algorithms.We select multiple task sizes from each workflow in the experiments.Besides, the DAG diagrams of some workflow instances with around 30 tasks are illustrated in Fig. 3.It is clear that these workflow instances cover various complicated structures, including in-tree, out-tree, fork-join, pipeline, and mixture.For more details on these workflows, please refer to the Pegasus repository. 3ive types of resources instances, i.e., t3.nano, t3.micro, t3.small, t3.medium, t3.medium, and t3.large, from Amazon EC2 4 are employed to simulate the cloud environments.We summarize the parameters of the five resource types in Table 1.Besides, we set the length of a billing period to 60 s and the bandwidth among resource instances to 5.0 Gbps.
Hypervolume [45] metric is designed to measure the quality of a population concerning both convergence and diversity, and has been frequently used to evaluate the performance of multi-objective evolutionary algorithms.Assume r = {r 1 , r 2 } is a reference point.The hypervolume value of a population P, corresponding to the volume between the reference point and the objective vectors of the solutions in  P, can be calculated as follows.
where L( ) represents the Lebesgue measure.
Referring to the settings in MOELS and EMS-C, we set the population size of the five algorithms as 100.The maximum number of fitness evaluations is set as the stop condition for all the five algorithms and is set as 500,000.
The five algorithms are independently repeated 31 times on each workflow instance to mitigate random effects.We run all the experiments on a PC with two Cores i7-6500U CPU @ 2.50 GHz 2.59 GHz, 8.00 GB RAM, Windows 10.

Comparison result
Tables 2 and 3 report the average and standard deviation (in brackets) of the hypervolume values for the algorithms MOELS, EMS-C, WOF, LSMOF, and VCAES in scheduling the 27 workflow instances to cloud resources.For each workflow instance, the largest hypervolume value among the five algorithms is highlighted in bold.Besides, the Wilcoxon's rank-sum test with α = 0.05 is used to examine the significant differences between the comparison algorithm and the proposed VCAES in terms of hypervolume metric.The marks −, +, and ≈ represent that the comparison algorithm is significantly worse than, better than, and similar to VCAES.
The comparison results in Tables 2 and 3 show that the proposed VCAES performs significantly better than the four baseline algorithms.Specifically, the proposal generates larger hypervolume values than MOELS, EMS-C, WOF, and LSMOF on 20, 22, 26, and 25 out of the 27 workflow instances respectively.
An interesting phenomenon is that the more decision variables, the more pronounced the advantages of the proposed VCAES.For example, in workflow instances with about 30 tasks, the VCAES obtains larger hypervolume values than algorithm LSMOF on 3 out of 5. Whereas, in workflow  To intuitively compare the convergence and diversity of the five multi-objective workflow scheduling algorithms, Fig. 4 illustrates the distribution of their output populations on workflow instances Inspiral, Montage, Sipht, and Cycles with large-scale tasks.
As illustrated in Fig. 4a, in the context of Inspiral with 100 tasks, the distribution of VCAES's output population is better than that of three baseline algorithms, i.e., MOELS, EMS-C, and WOF.More specifically, the output solutions of VCAES dominate that of the three algorithms.Although the diversity of solutions obtained by algorithm LSMOF is superior to that of the VCAES when the cost is less than 10, it is far inferior to algorithm VCAES during other range.To sum up, the proposed VCAES is superior to each baseline algorithm in terms of convergence and diversity.In the context of Montage with 100 tasks, the proposed VCAES has similar advantages over the comparison algorithms in terms of convergence and diversity, as shown in Fig. 4b.As can be observed from Fig. 4c, in the context of Sipht with 100 tasks, the solutions obtained by algorithm VCAES dominate most solutions obtained by comparison algorithms.This means that the VCAES outperforms all the four baseline algorithms in both convergence and diversity.In the context of Cycles with 657 tasks, the proposed VCAES has slight advantages over the comparison algorithms in terms of convergence and diversity, as shown in Fig. 4d.

Performance influence of different mechanisms
The VCAES mainly contains two new mechanisms to dynamically classify the decision variables and adaptively allocate evolution opportunities to each constructed group of decision variables.To measure the respective contributions of the two mechanisms to the overall performance, we construct three variants of the VCAES for comparison.The first variant, denoted as Variant 1, is constructed by replacing the decision variable classification mechanism with a random one and iterating each group of decision variables in a round-robin manner.The second variant, denoted as Variant 2, is constructed by replacing the decision variable classifica-tion mechanism with a random one and adaptively allocating evolution opportunities to each group of decision variables.The third variant, denoted as Variant 3, is constructed by removing the evolution opportunity allocation mechanism.
The main difference between the VCAES and its Variant 3 is that Variant 3 does not have the adaptive evolution opportunity allocation mechanism.Then, the improvement in the hypervolume value of the VCAES relative to Variant 3 can be attributed to the performance contribution of the adaptive evolution opportunity allocation mechanism.Similarly, the improvement in the hypervolume value of the VCAES relative to Variant 3 can be attributed to the performance contribution of the decision variable classification mechanism.The comparison results in Fig. 5 illustrate that the two proposed mechanisms contribute to the overall performance of the VCAES, with the decision variable classification mechanism contributing more.The performance improvement of the VCAES for Variant 1 can be attributed to the performance contribution of mixing the two proposed mechanisms.As shown in Fig. 5, in most workflow instances, the performance contribution of mixing the two mechanisms is better than that of any one.An exception is shown in Fig. 5b, where we can see that the overall performance contribution of mixing the two mechanisms is not as good as the contribution of the 123 decision variable classification mechanism.This is because a mechanism cannot be efficient in any scenario.Instead, it has certain advantages in some scenarios and inevitably has its disadvantages in other scenarios.

Conclusions
This paper focuses on two challenges in multi-objective cloud workflow scheduling: (1) large-scale decision variables; (2) and imbalance feature among variables regarding their contributions to objectives.To deal with these two challenges, this paper suggests a variable-contribution-based adaptive evolutionary cloud workflow scheduling approach that dynamically classifies the variables and adaptively allocates evolution opportunities to each constructed group of variables.Finally, in the context of real-world workflows and cloud platforms, this paper conducts comparison experiments to verify the effectiveness of the proposed adaptive mechanism in enhancing the population to approximate the Pareto-fronts of multi-objective cloud workflow scheduling problems.Cloud workflow scheduling is a representative grey-box problem, and it is interesting to mine the knowledge on the workflows and cloud resources to derive efficient scheduling algorithms.Another potential direction is to design a parallel evolutionary framework to shorten the time overhead of evolution optimization to support cloud workflow scheduling in real-time and uncertain situations.

Fig. 3
Fig. 3 DAG diagrams of workflows with about 30 tasks

3 )
instances with 100 and 1000 tasks, the algorithm proposed in this paper significantly outperforms all the four comparison algorithms.Similar to the baseline algorithms MOELS and EMS-C, the proposed VCAES also follows the framework of NSGA-II.Different from MOELS and EMS-C, the VCAES integrates a new adaptive strategy based on variable contribution.Although two large-scale multi-objective evolutionary algorithms based on problem transformation, i.e., WOF and LSMOF, exhibit competitive performance in solving continuous problems, their performance in multiobjective scheduling workflows is far inferior to the proposed VCAES.The above comparison results demonstrate that the proposed adaptive mechanism in this paper can effectively improve the performance of evolutionary algorithms in solving multi-objective cloud workflow scheduling problems with large-scale variables.

Fig. 4
Fig.4 Population distributions of the five algorithms on solving different workflow scheduling problems , it sorts the variables in a non-ascending

Table 1
Parameters for five types of cloud resource

Table 2
Comparison results for the five algorithms on 15 workflows in terms of the hypervolume metric

Table 3
Comparison results for the five algorithms on 15 workflows in terms of the hypervolume metric