For each genetic run in Algorithm 1, P0 is initialized from the search space S and the fitness value f(xi,t), ∀xi,t ∈ P0 is calculated. At each time step t, the population diversity of Pt is measured by calculating the standard deviation ζt as follows