ORIGINAL RESEARCH article

Front. Appl. Math. Stat., 17 October 2025

Sec. Optimization

Volume 11 - 2025 | https://doi.org/10.3389/fams.2025.1640044

EESB-FDO: enhancing the fitness-dependent optimizer through a modified boundary handling mechanism

  • 1. Department of Computer Science, College of Science, Charmo University, Sulaimani, Chamchamal, KR, Iraq

  • 2. Department of IT, Technical College of Informatics, Sulaimani Polytechnic University, Sulaimani, KR, Iraq

  • 3. Department of Software Engineering, College of Engineering and Computational Science, Charmo University, Chamchamal, Iraq

Abstract

The fitness-dependent optimizer (FDO) has recently gained attention as an effective metaheuristic for solving different optimization problems. However, it faces limitations in exploitation and convergence speed. To overcome these challenges, this study introduces two enhanced variants: enhancing exploitation through stochastic boundary for FDO (EESB-FDO) and enhancing exploitation through boundary carving for FDO (EEBC-FDO). In addition, the ELFS strategy is proposed to constrain Levy flight steps, ensuring more stable exploration. Experimental results show that these modifications significantly improve the performance of FDO compared to the original version. To evaluate the performance of the EESB-FDO and EEBC-FDO, three primary categories of benchmark test functions were utilized: classical, CEC 2019, and CEC 2022. The assessment was further supported by the application of statistical analysis methods to ensure a comprehensive and rigorous performance evaluation. The performance of the proposed EESB-FDO and EEBC-FDO algorithms was evaluated through comparative analysis with several existing FDO modifications, as well as with other well-established metaheuristic algorithms, including the Arithmetic Optimization Algorithm (AOA), the Learner Performance-Based Behavior Algorithm (LPB), the Whale Optimization Algorithm (WOA), and the Fox-inspired Optimization Algorithm (FOX). The statistical analysis indicated that both EESB-FDO and EEBC-FDO exhibit better performance compared to the aforementioned algorithms. Furthermore, a final evaluation involved applying EESB-FDO and EEBC-FDO to four real-world optimization problems: the gear train design problem, the three-bar truss problem, the pathological igg fraction in the nervous system, and the integrated cyber-physical attack on a manufacturing system. The results demonstrate that both proposed variants significantly outperform both the FDO and the modified fitness-dependent optimizer (MFDO) in solving these complex problems.

1 Introduction

In real-world scenarios, optimization plays a vital role in solving complex problems across diverse fields such as engineering design, healthcare, transportation, energy systems, and machine learning. Many of these problems are difficult to solve using traditional mathematical approaches, particularly when the search space is large, non-linear, or contains multiple local optima. Metaheuristic algorithms have therefore emerged as powerful tools, providing efficient and flexible methods to approximate optimal solutions in practical applications []. Over the past two decades, metaheuristic algorithms have gained prominence due to their ability to handle complex, high-dimensional problems. They generally operate by applying strategic rules to candidate solutions, starting from randomly generated points and iteratively improving them through evaluation and comparison []. One of the earliest milestones in this field was John Holland's genetic algorithm (GA), developed in the 1960s and later published in the 1970s and 1980s []. This was followed by simulated annealing (SA) [], inspired by the annealing process in metallurgy, and particle swarm optimization (PSO) [, ], inspired by the collective behavior of birds. These pioneering algorithms demonstrated strong potential for addressing real-world optimization challenges by finding high-quality solutions, reducing computational time, and overcoming local optima.

Since then, numerous algorithms have been introduced, such as the pathfinder algorithm (PFA) [], the mayfly algorithm (MA) [], the bear smell search algorithm (BSSA) [], the gradient-based optimizer (GBO) [], the group teaching optimization algorithm (GTOA) [], the heap-based optimizer (HBO) [], the Henry gas solubility optimization (HGSO) [], the marine predators algorithm (MPA) [], the political optimizer (PO) [] the atom search optimization (ASO) [], the Aquila optimizer (AO) [], the herd optimization algorithm (HOA) [], the Ebola optimization search algorithm (EOSA) [] and the single candidate optimizer (SCO) []. The large number of algorithms can be attributed to the fact that no specific algorithm can be applied to all optimization problems, as pointed out in the no free lunch (NFL) theorem []. This theorem has logically demonstrated that no single algorithm is universally optimal for solving all types of optimization problems. The outcomes of these algorithms show that some achieve better results than others for certain specific problems [].

A recurring theme in the design of these algorithms is the balance between exploration and exploitation []. Achieving a balance between these concepts is critical for improving algorithmic performance, as it enhances the ability of algorithms to converge efficiently while avoiding local optima. One of the primary limitations of the fitness-dependent optimizer (FDO) algorithm, as identified in previous studies, is its underdeveloped exploration capability, which has prompted significant research efforts aimed at addressing this issue. This paper emphasizes improving exploitation while maintaining a strong focus on exploration. To confine the algorithm within a specified search space defined by lower and upper bounds after updating the bees' positions, the FDO algorithm utilizes a specialized function called the “getBoundary” function [, ]. This function effectively transforms outlier values of bee positions into feasible values within the defined boundaries. By introducing modifications to the equations within the “getBoundary” function, this study aims to enhance the FDO algorithm's performance, contributing to more efficient optimization outcomes.

The main contributions of this study can be summarized as follows:

  • Proposed two novel boundary handling strategies for the Fitness fitness-dependent optimizer (FDO):

  • EESB-FDO: Enhancing exploitation through stochastic repositioning, which introduces random values within bounds when scout bees exceed search space limits.

  • EEBC-FDO: Enhancing exploitation through boundary carving, which redirects bees toward feasible regions using boundary carving equations.

  • 2. Introduced the ELFS strategy, which modifies the Levy flight mechanism to restrict step sizes within a bounded range, preventing instability due to excessive jumps.

  • 3. Conducted extensive benchmarking using classical, CEC 2019, and CEC 2022 test functions, demonstrating significant improvements in exploitation capability and convergence behavior.

  • 4. Applied the proposed variants to four real-world engineering problems, proving their practical its effectiveness in solving complex, constrained optimization scenarios.

  • 5. Performed comparative statistical analysis with state-of-the-art algorithms, confirming the competitiveness and reliability of EESB-FDO and EEBC-FDO under various optimization challenges.

The remaining sections of this paper are organized as follows: Section 2 presents a review of related work on the FDO. Section 3 provides a detailed description of the FDO algorithm. In Section 4, the proposed methodology is explained, followed by the introduction of the FDO modifications in Section 5. Section 6 presents the experimental results and comparisons with other algorithms using a set of benchmark test functions. Finally, Section 7 summarizes the main findings and offers directions for future research.

2 Related work

Advancing the FDO recently and effectiveness in solving real-world problems demonstrate that it is one of the most powerful metaheuristic algorithms in recent years []. Additionally, it uses random walk techniques, such as Levy flight [], which have been employed in prior studies to improve the performance of the optimization of other algorithms. Selecting an appropriate equation or method is crucial for efficiently reaching the global optimum while avoiding local optima. Following the development of any new algorithm, the modification and enhancement of existing algorithms using various methods and techniques are critical steps commonly undertaken by researchers to improve performance and adaptability. For example, the GA [] has been extended to the adaptive genetic algorithm (AGA) [], and PSO [] has evolved into the comprehensive learning PSO (CLPSO) []. Similarly, the ant colony optimization (ACO) [] method was improved to the max–min ant system (MMAS) [], while the FA [] was adapted into the adaptive firefly algorithm (AFA) [], and AOA [] was refined into the enhanced AOA (EAOA) [].

The enhancement and generation of algorithms are increasing rapidly; this rapid expansion has necessitated the need to manage and categorize them. There are various ways to classify these algorithms, and one common approach is based on their source of inspiration, which can be either natural or non-natural. Nature-inspired metaheuristic algorithms can be broadly categorized into four main classes: physics-based algorithms, evolutionary algorithms, human-based algorithms, and swarm-based algorithms []. The FDO is a metaheuristic algorithm that belongs to the swarm-based algorithms group, inspired by the bee swarming reproductive process. More specifically, it can be seen as an extension of PSO []. Various methods and techniques have been employed to enhance and improve the performance of the FDO algorithm, and these improvements must be acknowledged. Adding two parameters (alignment and cohesion) to update the position of search agents in addition to the existing pace factor []. Ten chaotic maps have been used for the initialization population, which is another improvement implemented in FDO []. After that, updating pace (velocity) by utilizing the sine–cosine scheme, modified pace-updating equations in the search phase, random weight factor, global fitness weight strategy, conversion parameter strategy, and the best solution-updating strategy[]. Then, modifying the fitness-dependent optimizer by changing waited factor range from [0,1] to [0,0.2] is another mechanism [].

With these developments and improvements, the FDO algorithm has been effectively applied to solve and address a variety of problems, demonstrating significant success in finding optimal solutions. For instance, the adaptive fitness-dependent optimizer has been applied to the one-dimensional bin packing problem, an NP-hard problem, producing more promising results compared to other approaches []. Subsequently, the automatic generation control (AGC) of a multi-source interconnected power system (IPS) was more effectively managed using FDO, reducing frequency overshoot/undershoot and settling time []. Then Improved-FDO was applied to control the same problem and provide better results than FDO []. Building on the advantages of FDO in controlling Automatic Generation Control (AGC), a new controller called the “modified PID controller” is introduced, which achieves more effective control than before, again utilizing FDO []. The performance of this algorithm was further leveraged to enhance another area, specifically in multilayer perceptron (MLP) neural networks, resulting in an increase in the average accuracy of the model []. FDO has demonstrated superior efficiency in solving non-linear optimal control problems (NOCPs) compared to the GA algorithm, effectively minimizing absolute error and providing better solutions [].

A new technique called generalized regression neural network combined with fitness-dependent optimization (GRNNFDO) has been used to control the behavior of the thermoelectric generator (TEG) system, and this technique produces more power compared to other algorithms []. The economic load dispatch (ELD) problem is another area where the improved FDO has been used to effectively reduce emission allocation and fuel costs []. Another effective method for classifying COVID-19 cases as positive or negative, where FDO combined with neural networks was tested on three different datasets, enhanced accuracy rates and minimized error rates []. A hybrid approach combining genetic algorithm and FDO (GA-MFDO) has been used to solve the NP-hard problem of workflow scheduling, achieving a significant p-value compared to other methods []. Similarly, the edge server allocation problem, an NP-hard challenge, was effectively addressed using the effective fitness-dependent optimizer (EFDO). When tested on Shanghai Telecom's dataset, which includes factors such as access delay, energy consumption, and workload balance, EFDO outperformed competing algorithms []. Finally, FDO integrated with neural networks has also been employed for detecting steel plate faults, achieving higher classification accuracy than other methods [57].

Levy flight is a type of random walk where step lengths follow a Levy distribution, characterized by its heavy-tailed probability distribution. It is one of the most effective approaches for generating a diverse and well-distributed random initialization set, thereby enhancing the exploration capabilities of the algorithm when metaheuristic algorithms begin with a randomly initialized set of values, structured according to the algorithm's design and procedural steps []. Levy flight is a fundamental component of the cuckoo search (CS) algorithm, which mimics the foraging behavior of certain animals and insects []. In this work, Levy flight is utilized to simulate insect movement, significantly boosting the global search capability of the flower pollination algorithm (FPA) []. For algorithms like the dragonfly algorithm (DA), Levy flight guides stochastic behavior, allowing for better traversal of the search space []. Additionally, the combination of differential evolution (DE) and Levy flight (LF) techniques has been shown to improve the performance of the harmony search algorithm (HSA) []. The slap swarm algorithm (SSA) also benefits from Levy flight mutation, which enhances its randomness and exploration capabilities during the search process []. Similarly, Levy flight-based mutation has been introduced to improve the exploration ability and overall efficiency of the slime mould algorithm (SMA) []. Furthermore, the use of Levy flight in the bat algorithm (BA) has been shown to improve its ability to escape local optima and achieve faster convergence [] Similarly, the incorporation of Levy flight in the gray wolf optimizer (GWO) enhances its global search ability, making it more effective for solving complex optimization problems [].

3 FDO algorithm

The FDO is inspired by the swarming behavior of bees during colony reproduction. When bees seek to establish a new hive, they follow a systematic set of rules to identify a suitable location. These decision-making processes within the bee community can be categorized into two fundamental components: searching for an appropriate site and moving toward it. Further details regarding FDO are provided through its parameter settings, as summarized in Table 1. These parameters significantly influence the algorithm's performance, and deviations beyond their predefined values may adversely affect the quality of the results.

Table 1

ParametersParameter value
Search agent30
Weight factorEither 0 or 1, FDO use “0” in all testing
Fitness weightRange between (0,1)
Levy flightRange between (1,−1)

FDO parameter settings.

In optimization algorithms, the search process is typically conducted within a predefined search space to identify optimal solutions. To ensure that the search remains within these boundaries, some algorithms implement a boundary handling mechanism. This mechanism involves iterating through potential solutions and adjusting any values that exceed the predefined limits by reassigning them to a random boundary value. This section provides a detailed exploration of these three key aspects: searching, movement, and boundary handling. Showed Pseudo code of FDO in Figure 1.

Figure 1

].

3.1 Bee searching

Scout bees initiate the search for an optimal hive by exploring various potential locations. This process represents a fundamental stage in the algorithm, which begins with the random generation of the scout bee population within a predefined search space, denoted as Xi(where i = 1, 2, 3, …, n). Each scout bee's position corresponds to a newly discovered solution within the search space. To identify a more suitable hive, scout bees perform a stochastic search for alternative locations. If a newly identified position offers a more favorable solution, the previous locations are disregarded. Conversely, if the new position does not yield an improvement, the scout bees continue moving in their prior direction while disregarding the less optimal locations.

3.2 Bee movement

In the natural world, scout bees search for a hive randomly to identify a more suitable location for reproduction. Similarly, artificial scout bees adopt this principle by moving and updating their current position to a potentially improved one through the addition of a pace value. If the newly discovered solution proves to be superior to the previous one, the scouts proceed in the current direction; otherwise, they continue following their prior trajectory. The movement of artificial scout bees is determined using Equation 1.

where i represents the current search agent, t represents the current iteration, x represents an artificial scout bee (search agent), and pace is the movement rate and direction of the artificial scout bee. The fitness weight (fw) as expressed in Equation 2 is used to calculate and manage the pace when the direction of pace is completely dependent on a random mechanism.

where denotes the best global solution's fitness function value that has been revealed so far; denotes the current solution's value of the fitness function; and wf denotes a weight factor, randomly set between 0 and 1, which is used for controlling the fw.

Later, the algorithm considers some settings for fw, for instance, if fw = 1 or 0, and = 0, the algorithm sets the pace randomly according to Equation 3. On the other hand, if fw > 0 and fw < 1, then the algorithm generates a random number in the (−1, 1) range to make the scout search in every direction; when r < 0, pace is calculated according to Equation 4, and when r >= 1, pace is calculated according to Equation 5.

where r denotes a random number in the range of [−1, 1], Xi, t denotes the current solution, and denotes the global best solution achieved thus far. Among various applications for random numbers, the FDO selects Levy flight because it considers further stable movement via its fair distribution curve [].

3.3 Bee boundary

In optimization algorithms, the search space refers to the set of all possible solutions within which the algorithm seeks to identify the optimal solution. The algorithm explores this space to determine the most suitable outcome. Each algorithm operates within a defined search space, constrained by specific boundaries and limitations to ensure both efficiency and feasibility. The performance of an algorithm is significantly influenced by the structure of the search space, and maintaining the algorithm's operations within the prescribed boundaries is crucial for its success. For instance, the FDO utilizes mechanisms, such as Equations 6 and 7 to adjust values that exceed the defined range, ensuring that all solutions remain within the limited search space. This approach enhances the algorithm's ability to converge toward optimal solutions effectively.

where Xi, t+1 represent the new bee value after the update in each iteration, the upper bound represents the upper boundary value of search space, the lower bound represent the lower boundary value of search space, and Levy shows Levy flight to generate a random value between (−1,1).

4 Methodology

This section of the research outlines the study aimed at enhancing the FDO algorithm. It encompasses addressing exploitation limitations, improving FDO through adaptive bee update strategies, experimental setup, performance evaluation, and effectiveness assessment. The primary objective of this study is to improve the algorithm's exploitation capability, which has been identified as one of its main weaknesses, illustrated in Figure 2.

Figure 2

4.1 Addressing exploitation limitations

Optimization problems are prevalent across various specialized domains, where the primary objective is to identify the optimal solution among numerous possible alternatives. In the field of metaheuristic algorithms, enhancing both exploration and exploitation is essential for improving performance. Exploitation focuses on refining solutions within promising regions. The FDO has been observed to exhibit weaknesses in exploitation; theoretically, exploitation can be improved by concentrating the search around high-quality solutions. In the context of swarm-based algorithms, this involves reinforcing local search behavior while maintaining feasibility. The proposed modifications aim to increase solution density within promising regions, thereby enhancing the probability of convergence toward global optima. According to convergence theory in metaheuristics, increasing selective pressure around fitter solutions helps refine accuracy in late-stage optimization, limiting the algorithm's ability to converge efficiently toward optimal solutions. Therefore, this research aims to address and enhance the exploitation capability of the FDO algorithm to improve its overall effectiveness in solving complex optimization problems.

4.2 Improving FDO through adaptive bee update strategies

As discussed in the preceding sections, various methods and solutions have been employed to address this issue, contributing to a relative improvement in the exploitation capability of the FDO. This enhancement of the algorithm adds new progression to previous improvements by introducing some novel equations to update bees: that exceed the boundary of the search space. In this research, two new modifications are proposed for updating these exceeding bees, the first one, EESB-FDO, is designed to update bees, and the second, EEBC-FDO, also aims to update bees, both contributing to improving the exploitation capability of the FDO algorithm. From a theoretical standpoint, boundary handling mechanisms play a critical role in shaping the topology of the solution space. The proposed modifications replace rigid boundary clipping with stochastic repositioning or controlled correction strategies. These approaches improve the algorithm's capacity to recover from constraint violations while maintaining stable and coherent search dynamics. Such strategies are consistent with the principles of constrained optimization and contribute to preserving the continuity and efficiency of the search trajectory.

4.3 Experimental setup and performance evaluation

The algorithm was implemented in the MATLAB 2022 environment on a laptop equipped with a 10th-generation Intel Core i5 processor and 8 GB of RAM. After incorporating modifications, the algorithm was executed more than 30 times, each for 500 iterations. The performance was evaluated using a set of benchmark functions and statistical methods. The results demonstrate the improved efficiency of the algorithm with the proposed updates compared to its previous version.

4.4 Effectiveness assessment

The testing stage is crucial for evaluating the effectiveness and capability of the algorithm. This process consists of several key steps, including benchmark test functions, comparative analysis, statistical test methods, and real-world problem testing. Benchmark testing functions, such as classical, CEC2019, and CEC2022, are used to assess the performance and capability of the new improvements. Comparative analysis measures the enhanced algorithm's results against the original version and other existing algorithms. Additionally, statistical methods, including the Wilcoxon rank-sum test, standard deviation, and mean, are used to verify the significance, performance, consistency, and reliability of the algorithm. Real-world testing is conducted by applying the algorithm to practical problems. The testing phase ultimately confirms that the modifications have led to meaningful enhancements.

5 FDO modifications

In this section, this study focuses on three modifications: EESB-FDO, EEBC-FDO, and ELFS, which are proposed as extensions of the FDO. These modifications and FDO are substantive to parameter setting and change results by changing one of them []. The fundamental concept of FDO is inspired by the reproductive process and collective decision-making behavior exhibited by bees. Specifically, it mimics the way bees search for multiple available hives and collaboratively select the most optimal ones.

5.1 Enhancing FDO performance through scout bee boundary repositioning

The artificial scout bee operates on the same principle, searching for the optimal solution among multiple available options within the defined boundary search space. Maintaining the search within a restricted area is a crucial factor influencing the algorithm's performance and effectiveness. To enforce this condition, the FDO employs Equations 6 and 7. When the position of an artificial scout bee exceeds the predefined boundaries, the FDO updates its value according to these equations, ensuring that the search remains within the permitted limits. In this enhancement, two new modifications are proposed to update scout bees that exceed the upper or lower boundary, ensuring they remain within the predefined search space as described below.

5.1.1 Enhancing exploitation through stochastic boundary for FDO

The first modification employs Equation 8 to regulate values that surpass either the upper or lower boundary, meaning using just one Equation 8 instead of Equations 6, 7. By using this equation, EESB-FDO updates the bee value that exceeds the boundaries (upper and lower); the proposed approach utilizes a randomized value within a predefined fixed range. This involves generating random values within the defined search space limits, thereby introducing greater diversity and adaptability into the search process, as shown in Figure 3A. This approach leverages uniform random sampling within the valid bounds, which ensures all feasible regions are equally likely to be explored upon violation. This strategy preserves population diversity and prevents premature convergence caused by agents clustering near hard boundary limits. It also avoids stagnation by introducing controlled randomness, a key requirement for ergodicity in metaheuristic convergence.

Here, update Xi, t+1 represents the new bee value after update in each iteration, random (boundary) mean generate any random value in the predefined search space between the upper boundary and lower boundary, Levy flight to generate a random value between (−1,1).

Figure 3

5.1.2 Enhancing exploitation through boundary carving for FDO

The second modification utilizes Equation 9 to adjust values that exceed the upper boundary and Equation 10 to regulate values that fall below the lower boundary. By formulating two equations for updating scout bees, notable progress has been achieved in enhancing the algorithm's performance. These equations contribute to more effective handling of boundary violations and improve the overall exploitation capability of the algorithm, as explained in Figure 3B.

where updateXi, t+1 represents the new bee value after the update in each iteration, Xi, t+1 represents the old bee value that exceed boundaries, the upper bound represents the upper boundary value of the search space, the lower bound represents the lower boundary value of the search space, and Levy shows Levy flight to generate random values between (−1,1).

5.2 Eliminating levy flight shortcomings

Levy flight is a type of random walk where step lengths follow a Levy distribution, known for its heavy-tailed probability distribution, as shown in the pseudo code in Figure 4. This property makes Levy flight highly effective for exploring large and unknown search spaces [].

Figure 4

Levy flight represents a type of random walk and plays a crucial role in metaheuristic algorithms by generating random values within the range of (−1,1). However, these values should not exceed the defined limits. To assess the behavior of Levy flight, it was executed 30 times, each with 10,000 iterations. The results, as presented in Figure 5, indicate that in each execution, some randomly generated values exceeded the specified range of (−1,1). To address this issue, Equation 11 was incorporated into the Levy flight mechanism. This equation ensures that any generated value exceeding the upper limit of 1 is adjusted to 1, while any value falling below the lower limit of −1 is set to −1. In Figure 4, explaining Levy flight before enhancement, it was shown that O was represented by (0.01 × step) and then changed to Equation 11 after the update.

where O represents the output of the Levy flight function, Equation 11 ensures that the generated values remain within the defined range of (−1,1) without exceeding these limits.

Figure 5

6 Result and discussion

The key point of algorithm evaluation is the comparison with other established optimization methods. To ensure a comprehensive assessment, various evaluation metrics are employed, including benchmark functions, statistical methods, and real-world problem scenarios. In this study, the proposed EESB-FDO and EEBC-FDO were compared with seven metaheuristic algorithms: FDO [], arithmetic optimization algorithm (AOA) [], gravitational search algorithm (GSA) [66], GWO [], learner performance-based behavior algorithm (LPB) [67], whale optimization algorithm (WOA) [68], and fox-inspired optimization algorithm (FOX) [69] using classical benchmark functions, CEC-2019, and CEC-2022. The performance evaluation is conducted based on statistical measures such as mean, standard deviation, and the Wilcoxon rank-sum test to ensure the performance and capability. To ensure fair and unbiased comparisons, all algorithms were evaluated under identical experimental conditions. Each algorithm was executed independently for 30 runs to account for the stochastic nature of metaheuristic methods. In every run, a population of 30 search agents was employed, and the optimization process was conducted for a maximum of 500 iterations. This consistent configuration across all algorithms ensures that observed performance differences are attributable to algorithmic effectiveness rather than variations in experimental settings. Following the execution, the average and standard deviation were calculated, and the results are presented in this section.

6.1 Classical benchmark test functions

The performance of EESB-FDO and EEBC-FDO was evaluated using three categories of test function groups, as outlined in Zervoudakis and Tsafarakis []. These test function groups possess distinct characteristics that assess different aspects of the algorithms' performance. Unimodal test functions (f1–f7) are employed to evaluate the algorithms' exploitation capability and convergence efficiency, as they contain only a single global minimum (or maximum) without any local optima, as shown in Table 2. In contrast, multimodal test functions (f8–f13) are utilized to assess the algorithms' exploration ability, ensuring they can escape local optima and effectively search for the global solution in complex landscapes with multiple optimal points, as presented in Table 3. Additionally, composite benchmark functions (f14–f23) are represent advanced test functions designed to simulate real-world optimization challenges, providing a more rigorous evaluation of the algorithms' robustness and adaptability in highly complex and irregular search spaces, as presented in Table 4.

Table 2

FunctionRangefmin
.[−100,100]0
.[−10,10]0
.[−100,100]0
TF4(x) = max {|x|, 1 ≤ in} [−100,100]0
.[−30,30]0
.[−100,100]0
.[−1.28,1.28]0

Unimodal benchmark functions [80].

Table 3

FunctionRangeShift positionfmin
.[−500, 500][−300, …-300]−418.9829
.[−5.12, 5.12][−2, −2, …-2]0
[−32, 32]0
.[−600, 600][−400, …−400]0
. u(xi, a, k, m)[−50, 50][−30, 30, …30]0
.[−50, 50][−100, …−100]0

Multi-modal benchmark functions [80].

Table 4

FunctionDimRangefmin
.2[−65,65]1
.4[−5,5]0.00030
.2[−5,5]−1.0316
.2[−5,5]0.398
.2[−2,2]3
.3[, ]−3.86
.6[0,1]−3.32
.4[0,10]−10.1532
.4[0,10]−10.4028
.4[0,10]−10.5363

Fixed-dimension multimodal benchmark functions [80].

All algorithms presented in Tables 510 were evaluated over 30 independent runs, each consisting of 500 iterations. The optimization process was initiated with 30 randomly generated search agents in 10, 20, and 30-dimensional search spaces for the first two test function groups (unimodal and multimodal functions). For the composite test function group, the dimensionality was set according to the structural requirements of each function. The performance of the algorithms was assessed by computing statistical measures, including the mean and standard deviation, to ensure a comprehensive evaluation of their optimization capabilities.

Table 5

TFMetricsEESB-FDOFDOLPBWOAGSAGWO
F1Mean1.53E-221.74E-186.27E-029.53E-761.83E-021.34E-27
SD7.14E-229.39E-185.06E-025.11E-759.86E-021.71E-27
F2Mean2.68E+002.34E+003.58E-021.02E-525.97E-029.16E-17
SD7.49E-016.92E-011.44E-025.15E-521.40E-016.35E-17
F3Mean2.23E-087.27E-082.06E+021.38E+029.95E+021.25E-05
SD4.13E-082.61E-071.43E+022.92E+024.74E+022.83E-05
F4Mean6.43E-086.90E-087.17E-014.05E+006.96E+001.01E-06
SD1.96E-075.88E-072.37E-016.31E+002.12E+001.09E-06
F5Mean1.46E+012.62E+015.43E+016.95E+009.31E+012.70E+01
SD2.96E+013.88E+013.78E+016.51E-011.22E+028.67E-01
F6Mean9.49E-191.30E-174.83E-029.19E-038.10E+007.64E-01
SD5.19E-187.13E-174.44E-024.13E-021.17E+013.50E-01
F7Mean6.29E-016.33E-011.57E-022.69E-038.01E-021.86E-03
SD3.12E-013.36E-019.18E-033.24E-034.31E-029.14E-04
F8Mean−4.23E+03−8.30E+03−3.71E+03−3.17E+03−2.55E+03−5.80E+03
SD3.96E+031.14E+042.43E+021.28E+033.79E+029.72E+02
F9Mean7.34E+007.05E+003.51E-022.69E+002.80E+012.32E+00
SD4.00E+003.11E+002.73E-028.84E+008.81E+003.70E+00
F10Mean7.07E-156.95E-151.19E-019.35E-151.29E-081.05E-13
SD1.21E-152.61E-156.49E-021.92E-143.83E-091.82E-14
F11Mean1.31E-011.56E-011.96E-016.21E-022.80E+012.84E- 03
SD9.74E-029.93E-021.75E-011.05E-016.13E+007.18E- 03
F12Mean8.66E+008.02E+003.01E-034.64E-031.62E+004.44E-02
SD4.29E+005.60E+006.13E-037.75E-031.01E+002.29E-02
F13Mean3.37E+003.70E+007.09E-031.27E-018.03E+006.46E-01
SD5.45E+004.10E+006.12E-031.82E-014.63E+002.01E-01
F14Mean2.52E+002.72E+009.98E-014.75E+005.09E+004.65E+00
SD1.32E+001.46E+004.44E-164.07E+004.24E+004.37E+00
F15Mean8.24E-041.35E-038.63E-031.04E-034.49E-033.80E-03
SD5.86E-042.06E-039.32E-031.63E-032.24E-037.41E-03
F16Mean−1.03E+00−1.03E+00−1.03E+00−9.97E-01−1.03E+00−1.03E+00
SD6.66E-166.66E-166.66E-161.85E-016.66E-166.66E-16
F17Mean3.98E-013.98E-013.98E-013.98E-013.98E-013.98E-01
SD1.67E-161.67E-163.16E-061.91E-051.67E-162.49E-06
F18Mean3.00E+003.00E+003.00E+003.00E+003.00E+003.00E+00
SD0.00E+000.00E+001.23E-045.12E-050.00E+004.00E-05
F19Mean−3.86E+00−3.86E+00−3.20E+00−3.85E+00−3.86E+00−3.86E+00
SD3.11E-153.11E-151.33E-151.89E-023.11E-152.32E-03
F20Mean−3.25E+00−2.86E+00−3.25E+00−3.22E+00−3.32E+00−3.27E+00
SD8.41E-021.14E+005.89E-028.60E-021.78E-156.69E-02
F21Mean−5.05E+00−5.05E+00−5.39E+00−8.94E+00−5.94E+00−8.55E+00
SD5.33E-043.61E-043.25E+002.23E+003.69E+002.48E+00
F22Mean−5.24E+00−4.70E+00−5.56E+00−8.86E+00−1.04E+01−1.02E+01
SD3.09E+003.38E+003.23E+002.78E+000.00E+009.54E-01
F23Mean−4.66E+00−4.32E+00−6.56E+00−6.75E+00−9.93E+00−1.02E+01
SD3.30E+002.91E+003.76E+004.60E+002.49E+001.35E+00

Results of comparing EESB-FDO with FDO and other chosen algorithms using classical benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

6.1.1 Classical benchmark test for enhancing exploitation through stochastic boundary for FDO

The results of EESB-FDO, FDO, LPB, WOA, GSA, and GWO, as shown in Table 5 for test functions TF1 to TF7, indicate that EESB-FDO generally outperforms FDO, with the exception of TF2. Specifically, EESB-FDO demonstrates superior performance over FDO, LPB, and GSA for TF1. Additionally, for TF3, TF4, and TF6, EESB-FDO outperforms all other algorithms. However, in TF2 and TF5, EESB-FDO shows comparatively poor performance, with other algorithms outperforming it. These results, by enhancing F1 to F7 static values, mean minimizing values and prove that the new modification in EESB-FDO significantly improves the algorithm's exploitation capability.

On the other hand, in TF8, TF11, TF14, and TF15, the results indicate that EESB-FDO outperformed FDO. However, when compared to other algorithms, the performance of these alternatives was slightly superior. Furthermore, in TF16, TF17, TF18, and TF19, both EESB-FDO and FDO achieved the optimal value and demonstrated comparable performance. In general, their results were more favorable than those of other algorithms.

In addition to the previous comparisons, Table 6 presents the results of EESB-FDO in comparison with M-IFDO, IFDO, SC-FDO, and CFDO. The findings indicate that EESB-FDO performed better than M-IFDO and IFDO in TF2 and TF4. In TF5, it outperformed all algorithms except SC-FDO, while in TF6, it demonstrated superior performance compared to the others. However, for TF1, TF3, and TF7, its results were the weakest among the algorithms. Furthermore, regarding multimodal functions, EESB-FDO achieved better performance than M-IFDO and IFDO in TF9, TF12, and TF13. However, in TF8, TF11, TF14, and TF15, its results fluctuated, showing slight weaknesses. Additionally, from TF16 to TF19, EESB-FDO successfully attained the optimal values, demonstrating its effectiveness in these test cases.

Table 6

TFMetricsEESB-FDOM_IFDOIFDOSC-FDOCFDO
F1Mean1.53E-223.44E-245.38E-240.00E+006.44E-51
SD7.14E-221.12E-232.74E-230.00E+003.46E-50
F2Mean5.95E-025.17E-017.34E-010.00E+004.02E-34
SD1.22E-012.71E-011.62E+000.00E+002.02E-33
F3Mean2.23E-081.06E-132.88E-070.00E+001.81E-06
SD4.20E-083.93E-136.90E-070.00E+005.60-06
F4Mean6.43E-085.00E-053.60E-040.00E+007.95E-06
SD1.96E-074.27E-038.11E-040.00E+004.27E-05
F5Mean1.46E+013.10E+011.84E+022.69E-011.68E+01
SD2.96E+014.00E+013.45E+022.69E-013.14E+01
F6Mean9.49E-194.15E+064.22E+062.02E-022.01E-01
SD5.19E-181.10E+038.15E-095.53E+025.02E-01
F7Mean6.29E-017.20E-015.68E-014.49-025.75E-01
SD3.12E-013.17E-013.14E-015.53E+023.54E-01
F8Mean−4.23E+03−3.00E+06−2.92E+06−1.02E+04−1.04E+04
SD3.96E+031.48E+052.34E+052.90E+041.29E+04
F9Mean7.34E+008.98E+001.35E+010.00E+005.64E-01
SD4.00E+009.85E+006.66E+000.00E+003.04E+00
F10Mean7.07E-153.89E-155.18E-153.26E-151.48E-15
SD1.21E-155.88E-161.67E-153.66E-151.32E-15
F11Mean1.31E-017.35E-025.26E-010.00E+009.26E-02
SD9.74E-023.91E-028.90E-020.00E+001.25E-01
F12Mean8.66E+001.75E+011.81E+014.54E-021.19E+00
SD4.29E+001.86E+012.57E+015.61E-027.97E-01
F13Mean3.37E+004.18E+094.10E+091.73E-016.79E-01
SD5.45E+003.13E+071.50E-051.90E-022.87E-01
F14Mean2.52E+008.50E-072.68E-076.36E+015.50E+00
SD1.32E+003.45E-054.68E-079.01E+013.92E+00
F15Mean8.24E-042.00E-194.03E-162.28E+023.10E-03
SD5.86E-046.73E-049.25E-162.37E+021.12E-02
F16Mean−1.03E+001.94E-169.14E-163.55E+02−1.03E+00
SD6.66E-163.34E-023.61E-163.61E+022.80E-03
F17Mean3.98E-012.20E+012.38E+015.17E+02−1.03E+00
SD1.67E-163.21E-011.24E-015.25E+023.01E+01
F18Mean3.00E+002.23E+ 022.24E+021.56E+021.93E+01
SD0.00E+001.34E-022.68E-051.68E+023.51E-03
F19Mean−3.86E+003.15E+013.15E+017.09E+02−2.64E+00
SD3.11E-157.90E-021.32E-037.31E+025.17E-01

Results of comparing EESB-FDO with modifications Of FDO using classical benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

The final evaluation of EESB-FDO, in comparison with FDO, was conducted using both 20- and 30-dimensional settings across unimodal and multimodal benchmark test functions, as presented in Table 7. In the 20-dimensional case, EESB-FDO outperformed FDO on several functions, specifically F1, F2, F3, F5, F6, F10, and F11, indicating a clear advantage in terms of both convergence speed and exploitation capability within this dimensional context. In the 30-dimensional case, EESB-FDO demonstrated superior performance compared to FDO on several benchmark functions, namely F1, F2, F5, F6, F9, F12, and F13. These results highlight the algorithm's enhanced convergence speed and exploitation capability, underscoring its effectiveness in handling complex optimization tasks within this dimensional setting. The results clearly indicate that several benchmark functions showed significant improvement, approaching near-optimal solutions, following dimension modification and the discovery of new optimal points.

Table 7

TFMetrics20-Dimensional30-Dimensional
EESB-FDOFDOEESB-FDOFDO
F1Mean5.51E-211.22E-197.13E-058.70E-05
SD1.5E-204.16E-191.08E-041.43E-04
F2Mean5.24E+005.71E+001.45E+011.65E+01
SD2.91E+003.00E+003.43E+003.49E+00
F3Mean7.40E-011.01E+004.40E+024.30E+02
SD4.75E-011.02E+002.93E+022.54E+02
F4Mean1.86E+001.45E+001.40E+011.47E+01
SD1.08E+007.90E-013.84E+003.08E+00
F5Mean3.43E+013.66E+017.33E+017.94E+01
SD2.89E+012.94E+014.66E+016.51E+01
F6Mean1.71E-202.01E-206.31E-047.95E-04
SD5.56E-208.10E-205.44E-046.60E-04
F7Mean6.07E-017.68E-018.97E-018.47E-01
SD2.92E-012.82E-013.32E-012.49E-01
F8Mean−5.19E+03−8.87E+03−8.20E+03−9.50E+03
SD4.46E+032.66E+049.83E+031.11E+04
F9Mean2.15E+012.27E+013.45E+013.70E+01
SD1.14E+018.05E+001.16E+011.18E+01
F10Mean1.78E+001.84E+002.18E+001.98E+00
SD9.21E-011.21E+001.18E+009.69E-01
F11Mean2.20E-023.65E-021.76E-021.53E-02
SD2.02E-022.80E-021.74E-021.73E-02
F12Mean6.55E+005.12E+003.13E+004.79E+00
SD4.60E+003.87E+002.79E+002.81E+00
F13Mean1.89E+002.53E+001.61E+003.98E+00
SD4.61E+005.68E+004.94E+006.29E+00

Results of comparing EESB-FDO with FDO for (20 and 30) dimensions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

6.1.2 Classical benchmark test for enhancing exploitation through boundary carving for FDO

For the second modification, EEBC-FDO, the results presented in Table 8 clearly demonstrate its superiority over the FDO in TF1 to TF7, with the exception of TF3, where it outperformed other algorithms but performed worse than FDO. These outcomes indicate an improvement in the algorithm's exploration capability. Additionally, in TF1 and TF2, EEBC-FDO performed better than LPB and GSA but was outperformed by WOA and GWO. In TF4, TF5, and TF6, EEBC-FDO outperformed all other algorithms, demonstrating its competitive advantage. However, from TF9 to TF15, its performance showed slight weaknesses, except for TF10, where it performed favorably compared to other algorithms. Nonetheless, in TF16 to TF19, both FDO and EEBC-FDO achieved equal performance, obtaining the optimal value. On the other hand, Table 9 presents a comparison of EEBC-FDO with previous modifications of FDO, including M-IFDO, IFDO, SC-FDO, and CFDO. The results indicate that the performance of EEBC-FDO is closely aligned with that of EESB-FDO when compared to these modifications.

Table 8

TFMetricsEEBC-FDOFDOLPBWOAGSAGWO
F1Mean2.78E-211.74E-186.27E-029.53E-761.83E-021.34E-27
SD1.50E-209.39E-185.06E-025.11E-759.86E-021.71E-27
F2Mean1.87E-032.34E+003.58E-021.02E-525.97E-029.16E-17
SD5.98E-036.92E-011.44E-025.15E-521.40E-016.35E-17
F3Mean5.67E-077.27E-082.06E+021.38E+029.95E+021.25E-05
SD2.83E-062.61E-071.43E+022.92E+024.74E+022.83E-05
F4Mean3.12E-085.09E-087.17E-014.05E+006.96E+001.01E-06
SD4.04E-085.88E-082.37E-016.31E+002.12E+001.09E-06
F5Mean1.96E+012.62E+015.43E+016.95E+019.31E+012.70E+01
SD3.15E+013.88E+013.78E+016.51E+ 011.22E+028.67E-01
F6Mean2.61E-201.30E-174.83E-029.19E-038.10E+007.64E-01
SD1.43E-197.13E-174.44E-024.13E-021.17E+013.50E-01
F7Mean5.53E-016.33E-011.57E-022.69E-038.01E-021.86E-03
SD3.35E-013.36E-019.18E-033.24E-034.31E-029.14E-04
F8Mean−4.20E+03−8.30E+03−3.71E+03−3.17E+03−2.55E+03−5.80E+03
SD3.86E+031.14E+042.43E+021.28E+033.79E+029.72E+02
F9Mean7.50E+007.05E+003.51E-022.69E+002.80E+012.32E+00
SD3.77E+003.11E+002.73E-028.84E+008.81E+003.70E+00
F10Mean6.73E-156.95E-151.19E-019.35E-151.29E-081.05E-13
SD3.03E-152.61E-156.49E-021.92E-143.83E-091.82E-14
F11Mean1.92E-011.56E-011.96E-016.21E-022.80E+012.84E- 03
SD1.42E-019.93E-021.75E-011.05E-016.13E+007.18E- 03
F12Mean7.72E+008.02E+003.01E-034.64E-031.62E+004.44E-02
SD5.53E+005.60E+006.13E-037.75E-031.01E+002.29E-02
F13Mean5.14E+003.70E+007.09E-031.27E-018.03E+006.46E-01
SD4.95E+004.10E+006.12E-031.82E-014.63E+002.01E-01
F14Mean2.47E+002.72E+009.98E-014.75E+005.09E+004.65E+00
SD2.90E+001.46E+004.44E-164.07E+004.24E+004.37E+00
F15Mean1.72E-031.35E-038.63E-031.04E-034.49E-033.80E-03
SD3.91E-032.06E-039.32E-031.63E-032.24E-037.41E-03
F16Mean−1.03E+00−1.03E+00−1.03E+00−9.97E-01−1.03E+00−1.03E+00
SD6.78E-166.66E-166.66E-161.85E-016.66E-166.66E-16
F17Mean3.98E-013.98E-013.98E-013.98E-013.98E-013.98E-01
SD1.69E-161.67E-163.16E-061.91E-051.67E-162.49E-06
F18Mean3.00E+003.00E+003.00E+003.00E+003.00E+003.00E+00
SD0.00E+000.00E+001.23E-045.12E-050.00E+004.00E-05
F19Mean−3.86E+00−3.86E+00−3.20E+00−3.85E+00−3.86E+00−3.86E+00
SD3.16E-153.11E-151.33E-151.89E-023.11E-152.32E-03
F20Mean−2.89E+00−2.86E+00−3.25E+00−3.22E+00−3.32E+00−3.27E+00
SD1.16E+001.14E+005.89E-028.60E-021.78E-156.69E-02
F21Mean−5.92E+00−5.05E+00−5.39E+00−8.94E+00−5.94E+00−8.55E+00
SD2.50E+003.61E-043.25E+002.23E+003.69E+002.48E+00
F22Mean−5.39E+00−4.70E+00−5.56E+00−8.86E+00−1.04E+01−1.02E+01
SD2.73E+003.38E+003.23E+002.78E+000.00E+009.54E-01
F23Mean−5.69E+00−4.32E+00−6.56E+00−6.75E+00−9.93E+00−1.02E+01
SD3.12E+002.91E+003.76E+004.60E+002.49E+001.35E+00

Results of comparing EEBC-FDO with FDO and other chosen algorithms Using classical benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

Table 9

TFMetricsEEBC-FDOM_IFDOIFDOSC-FDOCFDO
F1Mean2.78E-213.44E-245.38E-240.00E+006.44E-51
SD1.50E-201.12E-232.74E-230.00E+003.46E-50
F2Mean1.87E-035.17E-017.34E-010.00E+004.02E+34
SD5.98E-032.71E-011.62E+000.00E+002.02E-33
F3Mean5.67E-071.06E-132.88E-070.00E+001.81E-06
SD2.83E-063.93E-136.90E-070.00E+005.60-06
F4Mean3.12E-085.00E-053.60E-040.00E+007.95E-06
SD4.04E-084.27E-038.11E-040.00E+004.27E-05
F5Mean1.96E+013.10E+011.84E+022.69E-011.68E+01
SD3.15E+014.00E+013.45E+022.69E-013.14E+01
F6Mean2.61E-204.15E+064.22E+062.02E-022.01E-01
SD1.43E-191.10E+038.15E-095.53E+025.02E-01
F7Mean5.53E-017.20E-015.68E-014.49-025.75E-01
SD3.35E-013.17E-013.14E-015.53E+023.54E-01
F8Mean−4.20E+03−3.00E+06−2.92E+06−1.02E+04−1.04E+04
SD3.86E+031.48E+052.34E+052.90E+041.29E+04
F9Mean7.50E+008.98E+001.35E+010.00E+005.64E-01
SD3.77E+009.85E+006.66E+000.00E+003.04E+00
F10Mean6.73E-153.89E-155.18E-153.26E-151.48E-15
SD3.03E-155.88E-161.67E-153.66E-151.32E-15
F11Mean1.92E-017.35E-025.26E-010.00E+009.26E-02
SD1.42E-013.91E-028.90E-020.00E+001.25E-01
F12Mean7.72E+001.75E+011.81E+014.54E-021.19E+00
SD5.53E+001.86E+012.57E+015.61E-027.97E-01
F13Mean5.14E+004.18E+094.10E+091.73E-016.79E-01
SD4.95E+003.13E+071.50E-051.90E-022.87E-01
F14Mean2.47E+008.50E-072.68E-076.36E+015.50E+00
SD2.90E+003.45E-054.68E-079.01E+013.92E+00
F15Mean1.72E-032.00E-194.03E-162.28E+023.10E-03
SD3.91E-036.73E-049.25E-162.37E+021.12E-02
F16Mean−1.03E+001.94E-169.14E-163.55E+02−1.03E+00
SD6.78E-163.34E-023.61E-163.61E+022.80E-03
F17Mean3.98E-012.20E+012.38E+015.17E+02−1.03E+00
SD1.69E-163.21E-011.24E-015.25E+023.01E+01
F18Mean3.00E+002.23E+ 022.24E+021.56E+021.93E+01
SD0.00E+001.34E-022.68E-051.68E+023.51E-03
F19Mean−3.86E+003.15E+013.15E+017.09E+02−2.64E+00
SD3.16E-157.90E-021.32E-037.31E+025.17E-01

Comparing EEBC-FDO with modifications Of FDO using classical benchmark functions.

The use of higher-dimensional evaluations is a critical aspect of validating the effectiveness of EEBC-FDO. To this end, both 20- and 30-dimensional benchmark functions were employed to compare the performance of EEBC-FDO against FDO. In the 20-dimensional case, as shown in Table 10, EEBC-FDO outperformed FDO on functions F1, F4, F6, and F11, while also demonstrating superior average performance on functions F2, F3, F5, F7, F9, F10, and F13. These results confirm the robustness and improved optimization capabilities of EEBC-FDO in more complex, high-dimensional problem spaces. In the second case, using the 30-dimensional benchmark functions, EEBC-FDO exhibited some limitations. While it outperformed FDO on only three functions—F5, F6, and F13—it still achieved better average performance on functions F4, F7, F9, and F10. These results suggest that although EEBC-FDO maintains competitive performance in certain high-dimensional scenarios, its effectiveness may diminish as problem dimensionality increases further.

Table 10

TFMetrics20-Dimensional30-Dimensional
EEBC-FDOFDOEEBC-FDOFDO
F1Mean3.33E-211.22E-198.71E-058.70E-05
SD1.15E-204.16E-191.05E-041.43E-04
F2Mean5.08E+005.71E+001.41E+011.65E+01
SD3.25E+003.00E+005.06E+003.49E+00
F3Mean9.65E-011.01E+004.94E+024.30E+02
SD1.15E+001.02E+003.37E+022.54E+02
F4Mean9.92E-011.45E+001.28E+011.47E+01
SD5.19E-017.90E-014.04E+003.08E+00
F5Mean3.32E+013.66E+017.29E+017.94E+01
SD3.83E+012.94E+016.16E+016.51E+01
F6Mean1.60E-202.01E-206.30E-047.95E-04
SD2.55E-208.10E-205.81E-046.60E-04
F7Mean7.39E-017.68E-017.75E-018.47E-01
SD2.87E-012.82E-013.12E-012.49E-01
F8Mean−1.37E+07−8.87E+03−1.98E+09−9.50E+03
SD6.19E+072.66E+045.17E+091.11E+04
F9Mean2.25E+012.27E+013.56E+013.70E+01
SD9.44E+008.05E+001.42E+011.18E+01
F10Mean1.55E+001.84E+001.92E+001.98E+00
SD1.26E+001.21E+001.21E+009.69E-01
F11Mean2.25E-023.65E-022.09E-021.53E-02
SD2.30E-022.80E-022.05E-021.73E-02
F12Mean4.13E+005.12E+004.88E+004.79E+00
SD5.09E+003.87E+003.23E+002.81E+00
F13Mean2.04E+002.53E+001.53E+003.98E+00
SD6.86E+005.68E+005.14E+006.29E+00

Results of comparing EEBC-FDO with FDO for 20- and 30-dimensional algorithms using classical benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

6.2 CEC2019 benchmark test functions

Testing EESB-FDO and EEBC-FDO on a single benchmark test function is not a sufficient measure to ensure their effectiveness. Therefore, to thoroughly evaluate their performance and capability, the algorithm was tested on the CEC2019 benchmark functions presented in Table 11, which are commonly used for assessing optimization algorithms. This benchmark set includes ten multimodal test functions, which serve as enhanced evaluation functions for optimization purposes [70].

Table 11

No.FunctionsDimRangefmin
CEC01Storn's Chebyshev polynomial fitting problem9[−8,192, 8,192]1
CEC02Inverse Hilbert matrix problem16[−16,384, 16,384]1
CEC03Lennard–Jones minimum energy cluster18[−4, 4]1
CEC04Rastrigin's function10[−100, 100]1
CEC05Griewangk's function10[−100, 100]1
CEC06Weierstrass function10[−100, 100]1
CEC07Modified Schwefel's function10[−100,100]1
CEC08expanded Schaffer's F6 function10[−100,100]1
CEC09Happy cat function10[−100, 100]1
CEC10Ackley function10[−100, 100]1

CEC-C06 2019 benchmark functions [80].

The results of the proposed EESB-FDO, presented in Table 12, demonstrate its success in outperforming other algorithms, achieving the highest ranking in CEC02, CEC04, CEC05, CEC06, and CEC09. Additionally, in CEC03, CEC08, and CEC10, EESB-FDO obtained the best performance based on the average calculation. Furthermore, the LPB algorithm secured second place and showed competitive performance, particularly in CEC03 and CEC07. Meanwhile, WOA and FOX exhibited similar performance, ranking third overall. Finally, AOA demonstrated the weakest performance among the compared algorithms. On the other hand, Table 13 presents the results of EEBC-FDO, clearly indicating that FDO and EEBC-FDO did not achieve superior performance in the competition, as their results were nearly identical. However, overall, the outcomes demonstrated better performance compared to the other algorithms discussed in this section.

Table 12

CECMetricsEESB-FDOFDOAOALPBWOAFOX
CEC01Mean1.91E+091.74E+097.67E+097.03E+103.81E+101.83E+04
SD2.49E+091.77E+092.68E+107.40E+105.07E+102.25E+04
CEC02Mean1.73E+011.73E+011.93E+013.76E+011.74E+011.83E+01
SD1.08E-141.08E-143.76E-012.00E+019.95E-033.95E-04
CEC03Mean1.27E+011.27E+011.27E+011.27E+011.27E+011.37E+01
SD5.42E-155.42E-159.08E-045.33E-155.33E-157.11E-15
CEC04Mean3.73E+014.21E+011.09E+047.14E+013.62E+021.70E+03
SD1.26E+012.15E+013.35E+033.00E+011.28E+029.84E+02
CEC05Mean1.17E+001.20E+003.84E+001.27E+001.78E+006.48E+00
SD1.07E-011.35E-019.36E-011.29E-013.64E-011.32E+00
CEC06Mean1.08E+001.09E+019.02E+006.34E+009.43E+004.45E+00
SD8.50E-019.21E-011.22E+007.50E-011.20E+001.04E+00
CEC07Mean3.30E+023.47E+022.52E+022.68E+026.15E+023.86E+02
SD1.93E+022.12E+021.63E+021.66E+022.62E+022.48E+02
CEC08Mean5.40E+005.23E+005.57E+005.57E+005.91E+005.63E+00
SD6.67E-016.41E-015.31E-015.60E-016.13E-015.88E-01
CEC09Mean2.58E+002.62E+009.36E+023.09E+004.82E+003.77E+00
SD1.65E-012.98E-014.66E+022.97E-019.07E-013.55E-01
CEC10Mean1.93E+012.00E+012.01E+012.01E+012.01E+012.10E+01
SD3.59E+001.21E-026.09E-022.79E-021.10E+005.06E-03

Results of comparing EESB-FDO with FDO and other chosen algorithms Using CEC 2019 benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

Table 13

CECMetricsEEBC-FDOFDOAOALPBWOAFOX
CEC01Mean1.80E+091.74E+097.67E+097.03E+103.81E+101.83E+04
SD2.21E+091.77E+092.68E+107.40E+105.07E+102.25E+04
CEC02Mean1.73E+011.73E+011.93E+013.76E+011.74E+011.83E+01
SD1.08E-141.08E-143.76E-012.00E+019.95E-033.95E-04
CEC03Mean1.27E+011.27E+011.27E+011.27E+011.27E+011.37E+01
SD5.42E-155.42E-159.08E-045.33E-155.33E-157.11E-15
CEC04Mean4.51E+014.21E+011.09E+047.14E+013.62E+021.70E+03
SD2.76E+012.15E+013.35E+033.00E+011.28E+029.84E+02
CEC05Mean1.21E+001.20E+003.84E+001.27E+001.78E+006.48E+00
SD1.26E-011.35E-019.36E-011.29E-013.64E-011.32E+00
CEC06Mean1.10E+001.09E+019.02E+006.34E+009.43E+004.45E+00
SD1.10E+009.21E-011.22E+007.50E-011.20E+001.04E+00
CEC07Mean3.84E+023.47E+022.52E+022.68E+026.15E+023.86E+02
SD1.80E+022.12E+021.63E+021.66E+022.62E+022.48E+02
CEC08Mean5.14E+005.23E+005.57E+005.57E+005.91E+005.63E+00
SD6.20E-016.41E-015.31E-015.60E-016.13E-015.88E-01
CEC09Mean2.58E+002.62E+009.36E+023.09E+004.82E+003.77E+00
SD1.30E-012.98E-014.66E+022.97E-019.07E-013.55E-01
CEC10Mean2.00E+012.00E+012.01E+012.01E+012.01E+012.10E+01
SD3.15E-021.21E-026.09E-022.79E-021.10E+005.06E-03

Results of comparing EEBC-FDO with FDO and other chosen algorithms using CEC 2019 benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

6.3 CEC2022 benchmark test functions

To further evaluate the feasibility of the proposed algorithm, the CEC2022 benchmark test functions were utilized, as presented in Table 14. This benchmark set comprises 12 functions, labeled CEC01 to CEC12, which are organized into four groups: unimodal, multimodal, hybrid, and composite functions. The search space for these benchmark problems is defined within specific ranges of decision variables, constrained within the interval [−100,100] for all test functions. Therefore, d = 10 and 20 dimensions were used to assess the proposed modification [71].

Table 14

TypeNo.FunctionsDimRangefmin
Unimodal function1Shifted and full rotated Zakharov function10[−100,100]300
Basic functions2Shifted and full rotated Rosenbrock's function10[−100,100]400
3Shifted and full rotated expanded Schaffer's F6 function10[−100,100]600
4Shifted and full rotated non-continuous Rastrigin's function10[−100,100]800
5Shifted and full rotated Levy function10[−100,100]900
Hybrid functions6Hybrid function 1 (N = 3)10[−100,100]1,800
7Hybrid function 2 (N = 6)10[−100,100]2,000
8Hybrid function 3 (N = 5)10[−100,100]2,200
Composition functions9Composition function 1 (N = 5)10[−100,100]2,300
10Composition function 2 (N = 4)10[−100,100]2,400
11Composition function 3 (N = 5)10[−100,100]2,600
12Composition function 4 (N = 6)10[−100,100]2,700

CEC-C06 2022 benchmark functions [87].

The results of EESB-FDO, FDO, AOA, LPB, WOA, and FOX were evaluated by d = 10 dimension and are presented in Table 15. The outcomes of the proposed EESB-FDO modification indicate that it achieved the highest rank, outperforming other algorithms in CEC01, CEC02, CEC07, CEC09, and CEC12. Additionally, it obtained the best average performance in CEC05, CEC08, and CEC10. Meanwhile, the LPB algorithm secured second place, demonstrating superior performance in CEC03, CEC04, and CEC06. The FDO ranked third, while the remaining algorithms exhibited the lowest performance. To further evaluate the proposed modifications, Table 16 presents comparative results for both 10- and 20-dimensional problem instances. The outcomes demonstrate that EESB-FDO continues to perform successfully in the 20-dimensional setting, particularly on benchmark functions CEC01, CEC02, CEC05, CEC07, CEC08, CEC11, and CEC12, thereby confirming the robustness and scalability of the algorithm across higher-dimensional search spaces.

Table 15

CECMetricsEESB-FDOFDOAOALPBWOAFOX
CEC01Mean3.52E+033.93E+031.45E+048.77E+032.75E+043.99E+03
SD3.31E+033.39E+035.77E+035.59E+031.30E+043.56E+03
CEC02Mean4.07E+024.09E+021.33E+034.21E+024.57E+024.43E+02
SD1.77E+012.13E+016.85E+023.05E+011.04E+023.79E+01
CEC03Mean6.13E+026.14E+026.42E+026.00E+026.23E+026.55E+02
SD7.37E+008.94E+009.24E+002.25E-011.07E+029.58E+00
CEC04Mean8.25E+028.26E+028.37E+028.22E+028.44E+028.36E+02
SD7.73E+007.99E+007.36E+007.35E+001.70E+011.16E+01
CEC05Mean9.05E+029.57E+021.39E+031.14E+031.60E+031.48E+03
SD1.72E+021.06E+021.66E+022.63E+025.00E+025.49E+01
CEC06Mean4.41E+035.08E+031.57E+073.91E+036.84E+034.36E+03
SD3.16E+032.75E+034.11E+071.96E+034.91E+032.12E+03
CEC07Mean2.04E+032.04E+032.11E+032.04E+032.08E+032.15E+03
SD2.06E+012.32E+014.05E+014.60E+013.52E+018.26E+01
CEC08Mean2.23E+032.23E+032.28E+032.25E+032.23E+032.40E+03
SD2.26E+013.15E+018.27E+015.03E+015.94E+001.51E+02
CEC09Mean2.53E+032.53E+032.76E+032.56E+032.60E+032.63E+03
SD1.04E+012.68E+016.70E+013.43E+015.91E+014.91E+01
CEC10Mean2.55E+032.57E+032.78E+032.57E+032.71E+033.10E+03
SD6.35E+016.20E+012.82E+026.37E+013.73E+026.87E+02
CEC11Mean2.71E+032.68E+033.49E+032.80E+032.73E+033.10E+03
SD1.26E+021.30E+023.66E+021.44E+023.90E+026.87E+02
CEC12Mean2.88E+032.88E+033.03E+032.88E+032.90E+033.02E+03
SD1.06E+011.65E+018.58E+011.61E+014.39E+011.15E+02

Results of comparing EESB-FDO with FDO and Other chosen algorithms using CEC 2022 benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

Table 16

CECMetrics10-Dimensional20-Dimensional
EESB-FDOFDOEESB-FDOFDO
CEC01Mean3.52E+033.93E+034.88E+045.11E+04
SD3.31E+033.39E+031.55E+041.57E+04
CEC02Mean4.07E+024.09E+024.63E+024.67E+02
SD1.77E+012.13E+012.12E+012.33E+01
CEC03Mean6.13E+026.14E+026.47E+026.40E+02
SD7.37E+008.94E+008.71E+009.67E+00
CEC04Mean8.25E+028.26E+028.70E+028.85E+02
SD7.73E+007.99E+001.28E+012.20E+01
CEC05Mean9.05E+029.57E+022.16E+032.29E+03
SD1.72E+021.06E+026.04E+026.20E+02
CEC06Mean4.41E+035.08E+036.35E+035.45E+03
SD3.16E+032.75E+035.88E+034.62E+03
CEC07Mean2.04E+032.04E+032.14E+032.20E+03
SD2.06E+012.32E+016.33E+011.10E+02
CEC08Mean2.23E+032.23E+032.42E+032.48E+03
SD2.26E+013.15E+013.24E+024.01E+02
CEC09Mean2.53E+032.53E+032.49E+032.42E+03
SD1.04E+012.68E+016.20E+003.24E+02
CEC10Mean2.55E+032.57E+033.57E+033.37E+03
SD6.35E+016.20E+016.96E+029.67E+02
CEC11Mean2.71E+032.68E+032.94E+032.97E+03
SD1.26E+021.30E+024.78E+017.40E+02
CEC12Mean2.88E+032.88E+032.99E+033.09E+03
SD1.06E+011.65E+015.57E+021.09E+02

Results of comparing EESB-FDO with FDO for 10 and 20 dimensions using CEC 2022 benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

On the other hand, as illustrated in Table 17, the results indicate that the EEBC-FDO algorithm outperforms the other algorithms on the CEC03, CEC05, and CEC09 benchmark functions. Conversely, the remaining functions show superior performance by the other algorithms compared to EEBC-FDO. Subsequently, EEBC-FDO was evaluated using 20-dimensional benchmark functions, with the results presented in Table 18. The outcomes indicate that EEBC-FDO achieved superior performance, outperforming all competing algorithms on CEC01, CEC02, CEC06, CEC07, CEC08, and CEC11, thereby demonstrating its effectiveness and competitiveness in high-dimensional optimization tasks.

Table 17

CECMetricsEEBC-FDOFDOAOALPBWOAFOX
CEC01Mean3.84E+033.93E+031.45E+048.77E+032.75E+043.99E+03
SD4.17E+033.39E+035.77E+035.59E+031.30E+043.56E+03
CEC02Mean4.16E+024.09E+021.33E+034.21E+024.57E+024.43E+02
SD2.79E+012.13E+016.85E+023.05E+011.04E+023.79E+01
CEC03Mean6.12E+026.14E+026.42E+026.00E+026.23E+026.55E+02
SD7.94E+008.94E+009.24E+002.25E-011.07E+029.58E+00
CEC04Mean8.28E+028.26E+028.37E+028.22E+028.44E+028.36E+02
SD1.03E+017.99E+007.36E+007.35E+001.70E+011.16E+01
CEC05Mean9.47E+029.57E+021.39E+031.14E+031.60E+031.48E+03
SD5.13E+011.06E+021.66E+022.63E+025.00E+025.49E+01
CEC06Mean5.56E+035.08E+031.57E+073.91E+036.84E+034.36E+03
SD3.96E+032.75E+034.11E+071.96E+034.91E+032.12E+03
CEC07Mean2.04E+032.04E+032.11E+032.04E+032.08E+032.15E+03
SD2.68E+012.32E+014.05E+014.60E+013.52E+018.26E+01
CEC08Mean2.24E+032.23E+032.28E+032.25E+032.23E+032.40E+03
SD5.17E+013.15E+018.27E+015.03E+015.94E+001.51E+02
CEC09Mean2.53E+032.53E+032.76E+032.56E+032.60E+032.63E+03
SD9.25E-131.04E+016.70E+013.43E+015.91E+014.91E+01
CEC10Mean2.57E+032.55E+032.78E+032.57E+032.71E+033.10E+03
SD6.15E+016.20E+012.82E+026.37E+013.73E+026.87E+02
CEC11Mean2.70E+032.68E+033.49E+032.80E+032.73E+033.10E+03
SD1.23E+021.30E+023.66E+021.44E+023.90E+026.87E+02
CEC12Mean2.88E+032.88E+033.03E+032.88E+032.90E+033.02E+03
SD2.72E+011.65E+018.58E+011.61E+014.39E+011.15E+02

Results of comparing EEBC-FDO with FDO and Other chosen algorithms using CEC 2022 benchmark function.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

Table 18

CECMetrics10 Dim20 Dim
EEBC-FDOFDOEEBC-FDOFDO
CEC01Mean3.84E+033.93E+034.94E+045.11E+04
SD4.17E+033.39E+031.56E+041.57E+04
CEC02Mean4.16E+024.09E+024.64E+024.67E+02
SD2.79E+012.13E+011.72E+012.33E+01
CEC03Mean6.12E+026.14E+026.40E+026.40E+02
SD7.94E+008.94E+001.13E+019.67E+00
CEC04Mean8.28E+028.26E+028.86E+028.85E+02
SD1.03E+017.99E+001.98E+012.20E+01
CEC05Mean9.47E+029.57E+022.51E+032.29E+03
SD5.13E+011.06E+027.02E+026.20E+02
CEC06Mean5.56E+035.08E+034.57E+035.45E+03
SD3.96E+032.75E+033.65E+034.62E+03
CEC07Mean2.04E+032.04E+032.17E+032.20E+03
SD2.68E+012.32E+017.81E+011.10E+02
CEC08Mean2.24E+032.23E+032.39E+032.48E+03
SD5.17E+013.15E+011.36E+024.01E+02
CEC09Mean2.53E+032.53E+032.48E+032.42E+03
SD9.25E-131.04E+012.81E+003.24E+02
CEC10Mean2.57E+032.55E+033.46E+033.37E+03
SD6.15E+016.20E+019.32E+029.67E+02
CEC11Mean2.70E+032.68E+032.92E+032.97E+03
SD1.23E+021.30E+024.05E+017.40E+02
CEC12Mean2.88E+032.88E+033.16E+033.09E+03
SD2.72E+011.65E+011.93E+021.09E+02

Results of comparing EEBC-FDO with FDO for (10 and 20) dimensions using CEC 2022 benchmark function.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

6.4 Wilcoxon rank-sum test

One of the primary statistical methods used for hypothesis testing in this study is the Wilcoxon rank-sum test, which evaluates whether the performance differences between two algorithms are statistically significant. This non-parametric test is well-suited for optimization problems, as it does not assume a normal distribution of results and instead relies on rank comparisons. The test produces a p-value (probability value) that quantifies the strength of evidence against the null hypothesis. A smaller p-value (typically p < 0.05) indicates strong evidence to reject the null hypothesis, suggesting that one algorithm significantly outperforms the other, whereas larger p-values imply that the performance difference is not statistically significant [7274].

Table 19 presents the Wilcoxon test results comparing FDO, EESB-FDO, and EEBC-FDO against MFDO on the CEC2022 benchmark suite. The findings show that, for many functions (e.g., CEC01, CEC03, CEC04, CEC05), the extremely small p-values confirm that the proposed variants achieve statistically significant improvements. In contrast, higher p-values in functions such as CEC02 and CEC09 indicate no significant difference, suggesting that MFDO retains competitiveness in these landscapes, where exploration dominates over exploitation. Notably, EESB-FDO and EEBC-FDO reduce the number of non-significant cases compared to the original FDO, confirming the effectiveness of the proposed modifications in enhancing exploitation and convergence speed. These variations highlight that algorithm performance is problem-dependent: functions requiring stronger exploitation benefit from EESB-FDO and EEBC-FDO, while functions emphasizing broad exploration (e.g., CEC02, CEC09) show less pronounced differences. Overall, the Wilcoxon rank-sum test not only validates the robustness of the proposed variants but also provides insights into the problem characteristics where they achieve the most significant advantages.

Table 19

CECFDO_vs_MFDOEESB-FDO_vs_MFDOEEBC-FDO_vs_MFDO
CEC010.0000020.0000020.000002
CEC020.9426040.7343250.110878
CEC030.0000020.0000020.000002
CEC040.0000110.0000060.000003
CEC050.0000030.0000370.000002
CEC060.0000020.0000310.000004
CEC070.0000630.0000280.000022
CEC080.0000120.0000310.000002
CEC090.6262210.1940590.367185
CEC100.2989440.0284860.068714
CEC110.6620040.0078160.007734
CEC120.1779070.0407020.393334

WILCOXON rank-sum test for CEC2022 benchmark functions.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

6.5 Quantitative measurement metrics

Quantitative analysis is a crucial metric that provides detailed insights and in-depth observations of new algorithms and their modifications. Figures 68 illustrate this analysis, where the convergence of iterations represents the measurement of an agent's global best performance. Among the unimodal test functions (TF1–TF7), TF1 was selected as the representative function. For multimodal test functions (TF8–TF13), TF11 was chosen, while for composite test functions (TF14–TF19), TF15 was nominated. To further investigate performance, a two-dimensional search space was explored over 150 iterations using 10 search agents. As the number of iterations increased, (the global best agent) demonstrated greater precision. Additionally, when the scout bee emphasized exploitation and local search, the proposed modifications achieved superior results within fewer iterations, thereby demonstrating the algorithm's adaptability and effectiveness.

Figure 6

Figure 7

Figure 8

6.6 Real-world applications of EESB-FDO and EEBC-FDO

6.6.1 Gear train design problem

The primary objective of the gear train design (GTD) problem is to minimize the cost associated with the transmission ratio of the gear train configuration illustrated in Figure 9 [75, 76]. This optimization problem involves four discrete decision variables, Ta, Tb,Td, and Tf, each corresponding to the number of teeth on a specific gear within the train. The problem is formulated mathematically as follows:

Figure 9

6.6.2 Three-bar truss problem

As illustrated in Figure 10, the three-bar truss structure is considered for structural optimization. The objective is to minimize the total volume of the truss while ensuring that the design adheres to specified stress constraints [77, 78]. The corresponding mathematical formulation is provided below.

where .

Figure 10

6.6.3 The pathological IgG fraction in the nervous system

This real-world problem, described in Ismaeel Ghareb et al. [79], aims to determine the optimal solution for pathological conditions in humans or animals. It is known as the dimmer solution and is used to identify the IgG fraction in the nervous system. The mathematical formulation is presented below.

where Xi represents the input albumin quotient values, and Y(Xi) denotes the corresponding optimized IgG values. These equations, supported by Figure 11, demonstrate a real-world application of statistical modeling and optimization in medical diagnostics, offering an improved framework for evaluating immune dysfunctions in the nervous system.

Figure 11

6.6.4 Integrated cyber-physical attack for manufacturing system

This real-world application, supported by Figure 12, is designed to identify potentially harmful nodes within network systems, with a particular focus on software and hardware engineering design phases. The primary objective is to optimize the detection and regulation of these nodes within a normalized range of [0, 1], as their presence can introduce critical vulnerabilities, including system bugs and risk-prone components [80].

Figure 12

The fitness function is defined as

where X represents the node state values, and the coefficients A, B, and C are functions of system parameters, calculated as

in these equations, d denotes the total number of nodes (ranging from 15 to 36), while k1∈[0, 1] and k2∈[0.1, 0.5] are system-specific constants reflecting infection rate and risk thresholds. The aim of this fitness function is to identify configurations in which F(I(p, t)) → 0, indicating an optimal system state with minimal harmful node influence.

6.6.5 Real-world application comparison and discussion

The final step in evaluating any newly proposed algorithm or enhancement involves its application to real-world problems. To assess performance, standard evaluation metrics such as the mean, standard deviation, and p-value are employed for comparative analysis. In this section, four real-world applications were utilized to fulfill this objective. The corresponding results are presented in Tables 2022.

Table 20

Real-world problemsMetricsFDOEESB-FDOMFDO
Pathological IgG fraction in the nervous systemMean4.91E+004.90E+004.91E+00
SD1.98E-021.50E-021.47E-02
Integrated cyber-physical attack for manufacturing systemMean3.46E-04−3.72E+054.52E-04
SD1.44E-041.91E+061.86E-04
Three-bar trussMean1.62E+026.03E+01−2.39E+02
SD1.97E+021.48E+029.23E+02
Gear train designMean1.08E-074.31E-08−2.88E+04
SD3.02E-072.32E-073.27E+04

The results of comparison FDO with EESB-FDO and MFDO.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

Table 21

Real-world problemsMetricsFDOEEBC-FDOMFDO
Pathological IgG fraction in the nervous systemMean4.91E+004.79E+004.91E+00
SD1.98E-021.93E-011.47E-02
Integrated cyber-physical attack for manufacturing systemMean3.46E-043.30E-044.52E-04
SD1.44E-042.96E-041.86E-04
Three-bar trussMean1.62E+021.61E+02−2.39E+02
SD1.97E+021.15E+029.23E+02
Gear train designMean1.08E-074.14E-11−2.88E+04
SD3.02E-076.57E-113.27E+04

The results of comparison FDO with EEBC-FDO and MFDO.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

Table 22

Real-world problemFDO_vs_ EESB-FDOFDO_vs_ EEBC-FDO
Pathological IgG fraction in the nervous system0.0131940.000174
Integrated cyber-physical attack for manufacturing system0.0000020.149443
Three-bar truss0.0009790.005386
Gear train design0.2058880.049836

WILCOXON rank-sum test for four real-world problems.

The bold values in each row indicate the best performance among all compared algorithms for that test function (TF). Specifically, the lowest mean value represents the best solution quality*, while for the standard deviation (SD), the lowest value reflects the highest stability/consistency of results. *The lowest mean value represents the best solution quality.

The results of EESB-FDO and EEBC-FDO, presented in Tables 20, 21, demonstrate the performance of these two modified versions in comparison with the FDO and MFDO algorithms. In three real-world applications, both EESB-FDO and EEBC-FDO consistently achieved the best rankings, with lower mean and standard deviation values, indicating improved stability and reliability, outperforming the other algorithms. However, one of the modified versions exhibited a slight weakness in a specific case. On the other hand, the p-value, which quantifies the strength of evidence against the null hypothesis, indicates that a smaller p-value, typically p < 0.05, indicates strong evidence to reject the null hypothesis, suggesting that the observed results are unlikely to have occurred by chance. Table 22 shows the results of the comparison between FDO and EESB-FDO, which has just one case that fails to reject the null hypothesis for the gear train design, and FDO with EEBC-FDO, which shows the same result but with a different real-world problem's integrated cyber-physical attack for manufacturing system.

7 Conclusion

This study proposed two enhanced variants of the fitness-dependent optimizer (FDO) designed to address its limitations in exploitation and convergence speed. The first variant, EESB-FDO, introduces a stochastic boundary repositioning mechanism that reassigns scout bees that exceed the search space to random positions within the valid range. The second variant, EEBC-FDO, employs a boundary carving strategy that redirects out-of-bound solutions toward feasible regions using customized correction equations. In addition, a third enhancement, the ELFS strategy, constrains the Levy flight step range to (−1, 1), ensuring more stable and bounded exploration behavior. To validate the effectiveness of these modifications, extensive experiments were conducted using classical benchmark functions and CEC 2019 and CEC 2022 test suites. The results demonstrated that both EESB-FDO and EEBC-FDO significantly improve exploitation capability and convergence performance compared to the FDO and other state-of-the-art algorithms.

In the future, various enhancements and extensions of the current limited work can be developed. An important aspect to improve is the application of multi-objective vector optimization, utilizing functions specifically related to multi-objective optimization techniques [81]. A major limitation of this study lies in the inability to test certain benchmark functions in higher dimensions, which limits the generalizability of the findings. Future work could explore dimensionality reduction techniques or alternative benchmark suites that support high-dimensional testing. While classical benchmark functions can typically be tested across various dimensions, including 50, many CEC2019 functions are restricted and do not support high-dimensional testing. This constraint limits the generalizability and robustness analysis of the proposed algorithm. However, to partially address this issue, the CEC2022 benchmark set, which allows testing with 20 dimensions, can be used as an alternative for future evaluations. Despite these limitations, the enhancements should be applied to classification tasks and parameter tuning to improve the new model. This would also provide an opportunity to test the model on datasets related to education, healthcare, and other fields [79, 82]. Furthermore, the modification of this algorithm has not yet been validated on classification tasks involving real-world healthcare datasets, such as those related to sleep disorders or chronic disease diagnosis [83, 84]. As a result, its adaptability and robustness in managing noisy, imbalanced, and high-dimensional medical data remain uncertain due to the absence of dedicated benchmarking in medical classification scenarios. Although improvements have been demonstrated across several real-world applications, further development is needed by verifying dataset authenticity and incorporating security-focused algorithms to better guide the exploitation and exploration strategies in searching for new solutions [85, 86].

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

AF: Formal analysis, Funding acquisition, Methodology, Writing – original draft, Writing – review & editing, Conceptualization. AMA: Conceptualization, Methodology, Writing – review & editing, Resources. AAA: Conceptualization, Data curation, Writing – review & editing, Project administration.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Appendix

This appendix provides a detailed of the benchmark test functions employed to assess the performance of the proposed algorithms. The evaluation includes a comprehensive set of functions categorized into three main groups: classical benchmark functions (unimodal, multimodal, and composite types), as well as the CEC 2019 and CEC 2022 benchmark suites.

Summary

Keywords

stochastic boundary, optimization, fitness dependent optimizer, metaheuristic algorithm, Levy flight

Citation

Faraj AK, Aladdin AM and Ameen AA (2025) EESB-FDO: enhancing the fitness-dependent optimizer through a modified boundary handling mechanism. Front. Appl. Math. Stat. 11:1640044. doi: 10.3389/fams.2025.1640044

Received

03 June 2025

Accepted

17 September 2025

Published

17 October 2025

Volume

11 - 2025

Edited by

Betül Yildiz, Bursa Uludag Universitesi, Türkiye

Reviewed by

Zakariya Yahya Algamal, University of Mosul, Iraq

Lucas Wekesa, Strathmore University Institute of Mathematical Sciences, Kenya

Updates

Copyright

*Correspondence: Aram Kamal Faraj Aso M. Aladdin

†Present address: Aso M. Aladdin Department of Computer Engineering, Tishk International University, Erbil, KR, Iraq

†ORCID: Aso M. Aladdin orcid.org/0000-0002-8734-0811

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics