Many real-world optimization problems involve noisy objective evaluations and probabilistic constraints, particularly in the form of joint chance constraints, which are computationally expensive to evaluate. In this work, we propose CR-EA-C, a confidence-driven evolutionary algorithm for solving noisy black-box optimization problems under joint chance constraints. CR-EA-C introduces three key components: (1) analytical feasibility estimation for joint chance constraints, (2) a pairwise statistical ranking mechanism for robust comparison under noise, and (3) a modified infeasibility-driven survival strategy to accelerate convergence. These components enable statistically reliable decision-making while improving the efficiency of function evaluations. The proposed method is evaluated against four recent metaheuristic algorithms under various uncertainty distributions. Furthermore, its practical effectiveness is also assessed on two additional real-world optimization problems and compared with conventional static sampling methods. Experimental results show that CR-EA-C consistently satisfies the prescribed joint chance constraints while achieving competitive objective values overall. This demonstrates that CR-EA-C is an effective general-purpose approach for noisy optimization.
Figures & tables
Figure 1: Flowchart of CR-EA-C. The green boxes correspond to standard EA components, whereas the other elements correspond to the novel mechanisms introduced in this work.
Figure 2: Affect of sampling budget adaptation. n is sample size and p is probability of feasibility.
Parameter
Symbol
Value
Population size
N
10D
SBX crossover index
ηc
30
Probability of crossover
Pc
1.0
Polynomial mutation index
ηm
20
Probability of mutation
Pm
0.1
Minimum sample size
s0
3
Table 1: Parameter settings for CR-EA-C
Problem
Algorithm
Best
Worst
Mean
STDEV
Lower CI
Upper CI
AV
Example 1 Open Storage Network (minimize)
HPSO
64.085
109.347
73.883
10.476
71.431
76.335
0.0470
SSGA-I
61.592
73.436
64.693
2.734
64.053
65.333
0.0530
SSGA-II
61.898
67.721
64.420
1.659
64.032
64.808
0.0530
IOM
129.973
180.823
147.552
12.268
144.680
150.423
0.0000
ASIA
108.948
123.329
118.511
3.307
117.737
119.286
0.0000
CR-EA-C
105.692
131.212
116.7941
6.472
114.673
118.916
0.0004
Table 2: Performance summary for all algorithms on three problems [ 29 ] .
Figure 3: CR-EA-C median run’s surviving population of Example 2.
Figure 4: CR-EA-C median run’s sampling visualization of Example 2.
Figure 5: True rank of CR-EA-C’s first-rank solution in each generation for Examples 1–3.
Figure 6: CR-EA-C’s median run’s true rank of the first-rank solution in each generation for Examples 1–3.
Figure 7: CR-EA-C’s median run’s first-rank solution’s noiseless value in each generation for Examples 1–3.
Figure 8: CR-EA-C’s first-rank solution’s noiseless value in each generation of all runs for Examples 1–3.
Figure 9: The initial and last generation comparison with Examples 2.
Problem
ηc
Best
Worst
Mean
STDEV
Lower CI
Upper CI
AV
Example 1 Open Storage Network (minimize)
15
107.6604
131.8255
118.2605
6.8941
115.9546
120.5669
0.0008
20
110.1702
130.4887
119.3880
5.4513
117.6675
121.1079
0.0001
25
108.0109
129.6900
118.8864
7.1020
116.5107
121.2619
0.0012
30
105.6923
131.2129
116.7941
6.4725
114.6733
118.9167
0.0004
Example 2 Oil Production Planning (minimize)
15
136.4214
142.7890
138.3309
1.3243
137.9130
138.7490
0.0001
20
136.9151
149.9942
139.5526
2.8642
138.6643
140.4389
0.0000
Table 3: Parametric study results.
Figure 10: Median run’s first-rank solution’s p(x) value in each generation for Transportation problem (computed using large sample size of 300,000 evaluations). Dashed red line is at p(x)=0.95
Real-world optimization problems often involve black-box functions and uncertainties in their evaluation, widely referred to as noisy optimization problems (NOPs). Evolutionary algorithms (EA), including Evolutionary Strategies (ES) and genetic algorithms (GA) have been commonly adopted to solve these problems in the contemporary literature. An ongoing challenge is the computational expense involved, given the number of evaluations required for good fitness estimation and ranking. Two fundamental methods commonly used for fitness estimation for NOPs are implicit averaging and explicit averaging. Explicit averaging uses resampling of solutions to improve the estimates, while implicit averaging typically uses a large population size with low resampling. Implicit averaging has been shown to have theoretical advantages for certain cases, which has motivated some recent approaches to use them. However, a recent study demonstrated that its performance is highly dependent on certain assumptions about the function, such as steepness and constant noise level, which may not apply for majority of the real world problems. Moreover, most existing algorithms have only considered homoscedastic noise, where the amplitude of variation is uniform across the entire search space, as opposed to more generic case of heteroscedastic noise. To address these issues, we introduce a set of heteroscedastic test problems and propose a novel confidence ranking method that employs a computationally efficient explicit averaging strategy with sampling budget adaptation. It is implemented within the Covariance Matrix Adaptation ES (CMA-ES) and GA frameworks to demonstrate its effectiveness and versatility. The resulting algorithm is evaluated on a range of problems with both homoscedastic and heteroscedastic noise, and it demonstrates superior performance compared to state-of-the-art approaches.
Enrico Halim, Hemant Kumar Singh, Tapabrata Ray
The University of New South Wales, Canberra, Australia
Chance-constrained stochastic programs (CCSPs) require solutions that satisfy uncertain constraints with high probability. Standard approaches either embed scenarios directly into optimization by considering sample average approximation (SAA)-based reformulations, growing prohibitively expensive and often failing to find feasible solutions; or employ conservative reformulations that guarantee feasibility but sacrifice objective quality by design. We propose ACCESS (Adaptive Chance-Constrained Exploration via Structured Sampling), a structured sampling framework that decouples optimization from feasibility evaluation. First, high-quality candidate solutions are generated by solving many inexpensive single-scenario subproblems, avoiding sampling directly in the constrained high-dimensional solution space. Second, feasibility is assessed separately via large-scale batched simulation over a massive number of scenarios, which is naturally amenable to GPU parallelism. That reduces the search from a high-dimensional solution space to a two-dimensional parameter space controlling conservatism and scenario diversity; while concentrating on candidates near the feasibility boundary, where high-quality feasible solutions lie. We also introduce a variant, ACCESS(δ), which provides certified feasibility guarantees with probability 1−δ. Experiments across three benchmark classes and five scenario distributions show that ACCESS matches or exceeds the best objective found by any method in nearly every setting, while ACCESS(δ) is empirically feasible on every instance. The conservative reformulations remain several percent worse in objective, and SAA is frequently infeasible while taking an order of magnitude longer on the harder benchmarks.
Justin Dumouchelle, Shunyu Yao, Beste Basciftci +1
This paper presents novel method for distribution-free robust trajectory optimization and control of discrete-time, nonlinear, and non-Gaussian stochastic systems, with closed-loop guarantees on chance constraint satisfaction. Our framework employs conformal inference to generate coverage-based confidence sets for the closed-loop dynamics around arbitrary reference trajectories, by constructing a joint nonconformity score to quantify both the validity of contraction (i.e., incremental stability) conditions and the impact of external stochastic disturbance on the closed-loop dynamics, without any distributional assumptions. Via appropriate constraint tightening, chance constraints can be reformulated into tractable, statistically valid deterministic constraints on the reference trajectories. This enables a formal pathway to leverage and validate learning-based motion planners and controllers, such as those with neural contraction metrics, in safety-critical real-world applications. Notably, our statistical guarantees are non-diverging and can be computed with finite samples of the underlying uncertainty, without overly conservative structural priors. We demonstrate our approach in motion planning problems for designing safe, dynamically feasible trajectories in both numerical simulation and hardware experiments.
Rihan Aaron D'Silva, Hiroyasu Tsukamoto
Department of Aerospace Engineering, The Grainger College of Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois 61801