Note
Go to the end to download the full example code.
Multi-objective Binh-Korn example with the mNBI algorithm#
In this example, the modified Normal Boundary Intersection (mNBI) algorithm is used
to solve the BinhKorn
problem [BK97]:
\[\begin{split}\begin{aligned}
\text{minimize the objective function } & f_1(x, y) = 4x^2 + 4y^2 \\
& f_2(x, y) = (x-5)^2 + (y-5)^2 \\
\text{with respect to the design variables }&x,\,y \\
\text{subject to the general constraints }
& g_1(x,y) = (x-5)^2 + y^2 \leq 25.0\\
& g_2(x, y) = (x-8)^2 + (y+3)^2 \geq 7.7\\
\text{subject to the bound constraints }
& 0 \leq x \leq 5.0\\
& 0 \leq y \leq 3.0
\end{aligned}\end{split}\]
from __future__ import annotations
from numpy import array
from gemseo import configure_logger
from gemseo import execute_algo
from gemseo import execute_post
from gemseo.algos.opt.mnbi.settings.mnbi_settings import MNBI_Settings
from gemseo.problems.multiobjective_optimization.binh_korn import BinhKorn
configure_logger()
<RootLogger root (INFO)>
Solve the Binh-Korn optimization problem#
The 50 sub-optimization problems of mNBI are solved with SLSQP, a gradient-based optimization algorithm from the NLOPT library, with a maximum of 200 iterations. The analytic gradients are provided.
problem = BinhKorn()
mnbi_settings = MNBI_Settings(
max_iter=10000,
sub_optim_max_iter=200,
n_sub_optim=50,
sub_optim_algo="NLOPT_SLSQP",
)
result = execute_algo(problem, settings_model=mnbi_settings)
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize compute_binhkorn(x, y) = (4*x**2+ 4*y**2, (x-5.)**2 + (y-5.)**2)
INFO - 11:02:22: with respect to x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 11:02:22: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1 | 5 | float |
INFO - 11:02:22: | y | 0 | 1 | 3 | float |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: Solving optimization problem with algorithm MNBI:
INFO - 11:02:22: Searching for the individual optimum of each objective
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize f_0
INFO - 11:02:22: with respect to x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 11:02:22: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1 | 5 | float |
INFO - 11:02:22: | y | 0 | 1 | 3 | float |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 6
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: 0.0
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: ineq1 = [0.]
INFO - 11:02:22: ineq2 = [-59.3]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0 | 5 | float |
INFO - 11:02:22: | y | 0 | 0 | 3 | float |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize f_1
INFO - 11:02:22: with respect to x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 11:02:22: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1 | 5 | float |
INFO - 11:02:22: | y | 0 | 1 | 3 | float |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 4
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: 4.0
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: ineq1 = [-22.]
INFO - 11:02:22: ineq2 = [-7.3]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.999999999999934 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [133.2244898 4.93877551]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1 | 5 | float |
INFO - 11:02:22: | y | 0 | 1 | 3 | float |
INFO - 11:02:22: | t | -inf | 0 | inf | float |
INFO - 11:02:22: +------+-------------+-------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.01998643812378128
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [8.48468495e-09 2.12135087e-09]
INFO - 11:02:22: wrapped_ineq1 = [-21.98060064]
INFO - 11:02:22: wrapped_ineq2 = [-8.15508884]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.860718419626012 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999996 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.01998643812378128 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [130.44897959 5.87755102]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.860718419626012 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999996 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.01998643812378128 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.03911696914200459
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.07045392e-06 2.67613453e-07]
INFO - 11:02:22: wrapped_ineq1 = [-21.92182929]
INFO - 11:02:22: wrapped_ineq2 = [-9.0557122]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.72040975072545 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.03911696914200459 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [127.67346939 6.81632653]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.72040975072545 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.03911696914200459 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 6
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.057372362603395656
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.11930942e-06 2.79827355e-07]
INFO - 11:02:22: wrapped_ineq1 = [-21.82280187]
INFO - 11:02:22: wrapped_ineq2 = [-10.00289263]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.57905091675667 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.05737236260339566 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [124.89795918 7.75510204]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.57905091675667 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.05737236260339566 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 6
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.074732674761726
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.17118362e-06 2.92795885e-07]
INFO - 11:02:22: wrapped_ineq1 = [-21.68260071]
INFO - 11:02:22: wrapped_ineq2 = [-10.99769138]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.436617985224222 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.074732674761726 | inf | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [122.12244898 8.69387755]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.436617985224222 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.074732674761726 | inf | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 5
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.09117718665696596
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.22631340e-06 3.06578345e-07]
INFO - 11:02:22: wrapped_ineq1 = [-21.50027273]
INFO - 11:02:22: wrapped_ineq2 = [-12.04121071]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.293086093376062 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.09117718665696596 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [119.34693878 9.63265306]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.293086093376062 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.09117718665696596 | inf | float |
INFO - 11:02:22: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.10668436745322656
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.28495811e-06 3.21239565e-07]
INFO - 11:02:22: wrapped_ineq1 = [-21.27482752]
INFO - 11:02:22: wrapped_ineq2 = [-13.13459605]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.148429404209126 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999998 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1066843674532266 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [116.57142857 10.57142857]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.148429404209126 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999998 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1066843674532266 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 6
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.12123182918922926
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.34740275e-06 3.36850684e-07]
INFO - 11:02:22: wrapped_ineq1 = [-21.00523523]
INFO - 11:02:22: wrapped_ineq2 = [-14.27903845]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.002621052169886 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1212318291892293 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [113.79591837 11.51020408]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 4.002621052169886 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1212318291892293 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.1347962782336722
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.41395998e-06 3.53490002e-07]
INFO - 11:02:22: wrapped_ineq1 = [-20.69042436]
INFO - 11:02:22: wrapped_ineq2 = [-15.47577713]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.855633084900768 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1347962782336722 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [111.02040816 12.44897959]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.855633084900768 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1347962782336722 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.1473534631473843
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.48497408e-06 3.71243520e-07]
INFO - 11:02:22: wrapped_ineq1 = [-20.32927934]
INFO - 11:02:22: wrapped_ineq2 = [-16.72610225]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.707436400674776 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1473534631473843 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [108.24489796 13.3877551 ]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.707436400674776 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1473534631473843 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.15887811862176143
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.56082414e-06 3.90206031e-07]
INFO - 11:02:22: wrapped_ineq1 = [-19.92063796]
INFO - 11:02:22: wrapped_ineq2 = [-18.03135795]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.558000681121583 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1588781186217614 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [105.46938776 14.32653061]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.558000681121583 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1588781186217614 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 51/10000 [00:00<00:24, 411.87 it/sec, obj=[86.63747507 6.07936207]]
INFO - 11:02:22: 1%| | 52/10000 [00:00<00:23, 415.00 it/sec, obj=[82.45608407 6.53467094]]
INFO - 11:02:22: 1%| | 53/10000 [00:00<00:23, 417.77 it/sec, obj=[82.4386183 6.53671139]]
INFO - 11:02:22: 1%| | 54/10000 [00:00<00:23, 420.53 it/sec, obj=[82.43861798 6.53671142]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.16934390512642403
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.64192927e-06 4.10482285e-07]
INFO - 11:02:22: wrapped_ineq1 = [-19.46328861]
INFO - 11:02:22: wrapped_ineq2 = [-19.39294547]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.407294318804734 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.169343905126424 | inf | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [102.69387755 15.26530612]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.407294318804734 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.169343905126424 | inf | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: 1%| | 55/10000 [00:00<00:24, 413.39 it/sec, obj=[82.43861798 6.53671142]]
INFO - 11:02:22: 1%| | 56/10000 [00:00<00:24, 413.55 it/sec, obj=[82.4386183 6.53671139]]
INFO - 11:02:22: 1%| | 57/10000 [00:00<00:23, 416.36 it/sec, obj=[78.4046267 7.04173918]]
INFO - 11:02:22: 1%| | 58/10000 [00:00<00:23, 418.77 it/sec, obj=[78.38750452 7.04403274]]
INFO - 11:02:22: 1%| | 59/10000 [00:00<00:23, 421.35 it/sec, obj=[78.38750419 7.04403278]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.1787233438577463
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.72875264e-06 4.32188164e-07]
INFO - 11:02:22: wrapped_ineq1 = [-18.95596726]
INFO - 11:02:22: wrapped_ineq2 = [-20.8123267]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.25528433915987 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1787233438577463 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [99.91836735 16.20408163]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.25528433915987 | 5 | float |
INFO - 11:02:22: | y | 0 | 3 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1787233438577463 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 60/10000 [00:00<00:24, 410.43 it/sec, obj=[78.38750419 7.04403278]]
INFO - 11:02:22: 1%| | 61/10000 [00:00<00:24, 413.28 it/sec, obj=[74.50478473 7.6000843 ]]
INFO - 11:02:22: 1%| | 62/10000 [00:00<00:23, 415.04 it/sec, obj=[74.48803564 7.60264575]]
INFO - 11:02:22: 1%| | 63/10000 [00:00<00:23, 417.49 it/sec, obj=[74.48803531 7.6026458 ]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 6
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.18698774653318131
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.82180674e-06 4.55451708e-07]
INFO - 11:02:22: wrapped_ineq1 = [-18.39735425]
INFO - 11:02:22: wrapped_ineq2 = [-22.29102785]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.101936316247944 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1869877465331813 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [97.14285714 17.14285714]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 3.101936316247944 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1869877465331813 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 64/10000 [00:00<00:24, 407.55 it/sec, obj=[74.48803531 7.6026458 ]]
INFO - 11:02:22: 1%| | 65/10000 [00:00<00:24, 410.11 it/sec, obj=[70.75658544 8.21055954]]
INFO - 11:02:22: 1%| | 66/10000 [00:00<00:24, 412.67 it/sec, obj=[70.74019853 8.21247029]]
INFO - 11:02:22: 1%| | 67/10000 [00:00<00:23, 415.34 it/sec, obj=[70.7397555 8.21235896]]
INFO - 11:02:22: 1%| | 68/10000 [00:00<00:23, 417.48 it/sec, obj=[70.73968192 8.21237138]]
INFO - 11:02:22: 1%| | 69/10000 [00:00<00:23, 419.25 it/sec, obj=[70.73968188 8.21237139]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 8
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.19414099622670417
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [2.67250769e-07 6.68126852e-08]
INFO - 11:02:22: wrapped_ineq1 = [-17.92018718]
INFO - 11:02:22: wrapped_ineq2 = [-23.5380477]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.973627520197292 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.973627389560434 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1941409962267042 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [94.36734694 18.08163265]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.973627520197292 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.973627389560434 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1941409962267042 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 70/10000 [00:00<00:24, 412.52 it/sec, obj=[70.73968188 8.21237139]]
INFO - 11:02:22: 1%| | 71/10000 [00:00<00:23, 414.17 it/sec, obj=[70.73968192 8.21237138]]
INFO - 11:02:22: 1%| | 72/10000 [00:00<00:23, 416.54 it/sec, obj=[67.10688455 8.85599137]]
INFO - 11:02:22: 1%| | 73/10000 [00:00<00:23, 418.92 it/sec, obj=[67.09022377 8.85523197]]
INFO - 11:02:22: 1%| | 74/10000 [00:00<00:23, 421.43 it/sec, obj=[67.08836417 8.8547516 ]]
INFO - 11:02:22: 1%| | 75/10000 [00:00<00:23, 423.40 it/sec, obj=[67.0880443 8.8548097]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.20058310774576443
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.59261866e-08 3.98157241e-09]
INFO - 11:02:22: wrapped_ineq1 = [-17.67673508]
INFO - 11:02:22: wrapped_ineq2 = [-24.24810451]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.895860068518814 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.895860068939334 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2005831077457644 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [91.59183673 19.02040816]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.895860068518814 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.895860068939334 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2005831077457644 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 76/10000 [00:00<00:23, 414.36 it/sec, obj=[67.0880443 8.8548097]]
INFO - 11:02:22: 1%| | 77/10000 [00:00<00:23, 415.99 it/sec, obj=[63.53594358 9.5256767 ]]
INFO - 11:02:22: 1%| | 78/10000 [00:00<00:23, 418.08 it/sec, obj=[63.51929686 9.52496401]]
INFO - 11:02:22: 1%| | 79/10000 [00:00<00:23, 420.29 it/sec, obj=[63.51742882 9.52448661]]
INFO - 11:02:22: 1%| | 80/10000 [00:00<00:23, 422.11 it/sec, obj=[63.51711681 9.52454701]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.20643176428316765
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.60213212e-08 4.00534717e-09]
INFO - 11:02:22: wrapped_ineq1 = [-17.4199899]
INFO - 11:02:22: wrapped_ineq2 = [-24.97359042]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.817736612798091 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.817736606382494 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2064317642831676 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [88.81632653 19.95918367]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.817736612798091 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.817736606382494 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2064317642831676 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 81/10000 [00:00<00:23, 414.05 it/sec, obj=[63.51711681 9.52454701]]
INFO - 11:02:22: 1%| | 82/10000 [00:00<00:23, 416.05 it/sec, obj=[60.04682873 10.22303763]]
INFO - 11:02:22: 1%| | 83/10000 [00:00<00:23, 418.14 it/sec, obj=[60.03019841 10.22237318]]
INFO - 11:02:22: 1%| | 84/10000 [00:00<00:23, 420.17 it/sec, obj=[60.02832235 10.22189905]]
INFO - 11:02:22: 1%| | 85/10000 [00:00<00:23, 421.97 it/sec, obj=[60.0280183 10.22196177]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.21167873709715218
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.60891247e-08 4.02227762e-09]
INFO - 11:02:22: wrapped_ineq1 = [-17.14976701]
INFO - 11:02:22: wrapped_ineq2 = [-25.71472011]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.739252146962356 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.739252133472566 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2116787370971522 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [86.04081633 20.89795918]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.739252146962356 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.739252133472566 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2116787370971522 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 86/10000 [00:00<00:23, 416.42 it/sec, obj=[60.0280183 10.22196177]]
INFO - 11:02:22: 1%| | 87/10000 [00:00<00:23, 417.80 it/sec, obj=[60.0280183 10.22196177]]
INFO - 11:02:22: 1%| | 88/10000 [00:00<00:23, 419.75 it/sec, obj=[56.64068487 10.94846138]]
INFO - 11:02:22: 1%| | 89/10000 [00:00<00:23, 421.76 it/sec, obj=[56.62407339 10.94784674]]
INFO - 11:02:22: 1%| | 90/10000 [00:00<00:23, 423.62 it/sec, obj=[56.62218975 10.9473762 ]]
INFO - 11:02:22: 1%| | 91/10000 [00:00<00:23, 425.18 it/sec, obj=[56.62189378 10.94744126]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.21631560710603964
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.61276361e-08 4.03189482e-09]
INFO - 11:02:22: wrapped_ineq1 = [-16.86587782]
INFO - 11:02:22: wrapped_ineq2 = [-26.47171246]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.660401619556324 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.660401598766817 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2163156071060396 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [83.26530612 21.83673469]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.660401619556324 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.660401598766817 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2163156071060396 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 92/10000 [00:00<00:23, 418.40 it/sec, obj=[56.62189377 10.94744126]]
INFO - 11:02:22: 1%| | 93/10000 [00:00<00:23, 420.30 it/sec, obj=[53.31868379 11.70234427]]
INFO - 11:02:22: 1%| | 94/10000 [00:00<00:23, 422.09 it/sec, obj=[53.30209371 11.70178107]]
INFO - 11:02:22: 1%| | 95/10000 [00:00<00:23, 423.96 it/sec, obj=[53.30020299 11.70131446]]
INFO - 11:02:22: 1%| | 96/10000 [00:00<00:23, 425.35 it/sec, obj=[53.29991518 11.70138187]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22033375703666092
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.61356368e-08 4.03389500e-09]
INFO - 11:02:22: wrapped_ineq1 = [-16.5681293]
INFO - 11:02:22: wrapped_ineq2 = [-27.24479154]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.581179860208634 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.581179831908805 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2203337570366609 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [80.48979592 22.7755102 ]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.581179860208634 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.581179831908805 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2203337570366609 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 97/10000 [00:00<00:23, 418.82 it/sec, obj=[53.29991518 11.70138187]]
INFO - 11:02:22: 1%| | 98/10000 [00:00<00:23, 420.52 it/sec, obj=[50.08202514 12.48509203]]
INFO - 11:02:22: 1%| | 99/10000 [00:00<00:23, 422.19 it/sec, obj=[50.06545915 12.48458197]]
INFO - 11:02:22: 1%| | 100/10000 [00:00<00:23, 423.93 it/sec, obj=[50.06356185 12.48411964]]
INFO - 11:02:22: 1%| | 101/10000 [00:00<00:23, 425.47 it/sec, obj=[50.06328232 12.48418943]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22372436483038086
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.61119189e-08 4.02796374e-09]
INFO - 11:02:22: wrapped_ineq1 = [-16.25632384]
INFO - 11:02:22: wrapped_ineq2 = [-28.03418671]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.501581575675593 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.501581539671863 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2237243648303809 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [77.71428571 23.71428571]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.501581575675593 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.501581539671863 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2237243648303809 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 102/10000 [00:00<00:23, 418.67 it/sec, obj=[50.06328232 12.48418943]]
INFO - 11:02:22: 1%| | 103/10000 [00:00<00:23, 420.34 it/sec, obj=[46.93193733 13.29712012]]
INFO - 11:02:22: 1%| | 104/10000 [00:00<00:23, 421.95 it/sec, obj=[46.91539824 13.29666496]]
INFO - 11:02:22: 1%| | 105/10000 [00:00<00:23, 423.52 it/sec, obj=[46.91349492 13.2962073 ]]
INFO - 11:02:22: 1%| | 106/10000 [00:00<00:23, 424.99 it/sec, obj=[46.91322377 13.29627947]]
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22647839676434958
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.60554023e-08 4.01386124e-09]
INFO - 11:02:22: wrapped_ineq1 = [-15.93025908]
INFO - 11:02:22: wrapped_ineq2 = [-28.84013285]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.421601345713955 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.421601301831813 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2264783967643496 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [74.93877551 24.65306122]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.421601345713955 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.421601301831813 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2264783967643496 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 107/10000 [00:00<00:23, 418.72 it/sec, obj=[46.91322377 13.29627947]]
INFO - 11:02:22: 1%| | 108/10000 [00:00<00:23, 420.25 it/sec, obj=[43.86967855 14.13885408]]
INFO - 11:02:22: 1%| | 109/10000 [00:00<00:23, 421.78 it/sec, obj=[43.8531693 14.13845566]]
INFO - 11:02:22: 1%| | 110/10000 [00:00<00:23, 423.40 it/sec, obj=[43.85126056 14.13800305]]
INFO - 11:02:22: 1%| | 111/10000 [00:00<00:23, 424.77 it/sec, obj=[43.85099789 14.13807762]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22858660027279698
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.59653339e-08 3.99136368e-09]
INFO - 11:02:22: wrapped_ineq1 = [-15.58972776]
INFO - 11:02:22: wrapped_ineq2 = [-29.66287052]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.341233618773758 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.341233566859697 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.228586600272797 | inf | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [72.16326531 25.59183673]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.341233618773758 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.341233566859697 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.228586600272797 | inf | float |
INFO - 11:02:22: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:22: 1%| | 112/10000 [00:00<00:23, 419.29 it/sec, obj=[43.85099789 14.13807762]]
INFO - 11:02:22: 1%| | 113/10000 [00:00<00:23, 420.79 it/sec, obj=[40.89653773 15.01072985]]
INFO - 11:02:22: 1%| | 114/10000 [00:00<00:23, 422.26 it/sec, obj=[40.88006141 15.01039006]]
INFO - 11:02:22: 1%| | 115/10000 [00:00<00:23, 423.85 it/sec, obj=[40.87814788 15.00994293]]
INFO - 11:02:22: 1%| | 116/10000 [00:00<00:23, 425.16 it/sec, obj=[40.8778938 15.0100199]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.2300394964522294
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.58414366e-08 3.96036626e-09]
INFO - 11:02:22: wrapped_ineq1 = [-15.23451757]
INFO - 11:02:22: wrapped_ineq2 = [-30.50264619]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.260472707501503 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.260472647424826 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2300394964522294 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [69.3877551 26.53061224]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.260472707501503 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.260472647424826 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2300394964522294 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 117/10000 [00:00<00:23, 419.33 it/sec, obj=[40.8778938 15.0100199]]
INFO - 11:02:22: 1%| | 118/10000 [00:00<00:23, 420.73 it/sec, obj=[38.01383568 15.91319413]]
INFO - 11:02:22: 1%| | 119/10000 [00:00<00:23, 422.16 it/sec, obj=[37.99739551 15.91291495]]
INFO - 11:02:22: 1%| | 120/10000 [00:00<00:23, 423.66 it/sec, obj=[37.99547788 15.91247373]]
INFO - 11:02:22: 1%| | 121/10000 [00:00<00:23, 424.90 it/sec, obj=[37.99523249 15.91255313]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.23082737223334543
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.56837494e-08 3.92091692e-09]
INFO - 11:02:22: wrapped_ineq1 = [-14.86441091]
INFO - 11:02:22: wrapped_ineq2 = [-31.35971238]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.179312784043387 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.179312715697836 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2308273722333454 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [66.6122449 27.46938776]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.179312784043387 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.179312715697836 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2308273722333454 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%| | 122/10000 [00:00<00:23, 419.64 it/sec, obj=[37.99523249 15.91255313]]
INFO - 11:02:22: 1%| | 123/10000 [00:00<00:23, 420.99 it/sec, obj=[35.22292611 16.84670478]]
INFO - 11:02:22: 1%| | 124/10000 [00:00<00:23, 422.40 it/sec, obj=[35.20652549 16.84648825]]
INFO - 11:02:22: 1%|▏ | 125/10000 [00:00<00:23, 423.82 it/sec, obj=[35.20460449 16.84605342]]
INFO - 11:02:22: 1%|▏ | 126/10000 [00:00<00:23, 424.98 it/sec, obj=[35.20436789 16.84613524]]
INFO - 11:02:22: 1%|▏ | 127/10000 [00:00<00:23, 426.11 it/sec, obj=[35.20436789 16.84613524]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 8
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.2309402722013747
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.54922830e-08 3.87308674e-09]
INFO - 11:02:22: wrapped_ineq1 = [-14.47918481]
INFO - 11:02:22: wrapped_ineq2 = [-32.23432794]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.097747875137614 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.097747798442793 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2309402722013747 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [63.83673469 28.40816327]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.097747875137614 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.097747798442793 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2309402722013747 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%|▏ | 128/10000 [00:00<00:23, 421.12 it/sec, obj=[35.20436789 16.84613524]]
INFO - 11:02:22: 1%|▏ | 129/10000 [00:00<00:23, 422.43 it/sec, obj=[32.52519687 17.81173117]]
INFO - 11:02:22: 1%|▏ | 130/10000 [00:00<00:23, 423.73 it/sec, obj=[32.50883935 17.81157943]]
INFO - 11:02:22: 1%|▏ | 131/10000 [00:00<00:23, 425.04 it/sec, obj=[32.50691576 17.81115147]]
INFO - 11:02:22: 1%|▏ | 132/10000 [00:00<00:23, 426.13 it/sec, obj=[32.50668806 17.81123573]]
INFO - 11:02:22: 1%|▏ | 133/10000 [00:00<00:23, 427.28 it/sec, obj=[32.50668806 17.81123573]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.23036799004535574
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.52679220e-08 3.81699117e-09]
INFO - 11:02:22: wrapped_ineq1 = [-14.07861062]
INFO - 11:02:22: wrapped_ineq2 = [-33.12675824]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.015771856984091 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.015771771886646 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2303679900453557 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [61.06122449 29.34693878]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 2.015771856984091 | 5 | float |
INFO - 11:02:22: | y | 0 | 2.015771771886646 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2303679900453557 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%|▏ | 134/10000 [00:00<00:23, 422.00 it/sec, obj=[32.50668806 17.81123573]]
INFO - 11:02:22: 1%|▏ | 135/10000 [00:00<00:23, 423.24 it/sec, obj=[29.92207113 18.80875463]]
INFO - 11:02:22: 1%|▏ | 136/10000 [00:00<00:23, 424.43 it/sec, obj=[29.90576043 18.8086699 ]]
INFO - 11:02:22: 1%|▏ | 137/10000 [00:00<00:23, 425.63 it/sec, obj=[29.90383509 18.80824933]]
INFO - 11:02:22: 1%|▏ | 138/10000 [00:00<00:23, 426.62 it/sec, obj=[29.9036164 18.80833604]]
INFO - 11:02:22: 1%|▏ | 139/10000 [00:00<00:23, 427.61 it/sec, obj=[29.9036164 18.80833604]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 8
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22910005961555657
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.50115760e-08 3.75289311e-09]
INFO - 11:02:22: wrapped_ineq1 = [-13.66245391]
INFO - 11:02:22: wrapped_ineq2 = [-34.03727539]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.933378449878997 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.933378356353534 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2291000596155566 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [58.28571429 30.28571429]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.933378449878997 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.933378356353534 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2291000596155566 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%|▏ | 140/10000 [00:00<00:23, 424.33 it/sec, obj=[29.9036164 18.80833604]]
INFO - 11:02:22: 1%|▏ | 141/10000 [00:00<00:23, 425.14 it/sec, obj=[29.9036164 18.80833604]]
INFO - 11:02:22: 1%|▏ | 142/10000 [00:00<00:23, 426.31 it/sec, obj=[27.41500866 19.83826888]]
INFO - 11:02:22: 1%|▏ | 143/10000 [00:00<00:23, 427.48 it/sec, obj=[27.39874869 19.83825346]]
INFO - 11:02:22: 1%|▏ | 144/10000 [00:00<00:22, 428.72 it/sec, obj=[27.39682249 19.83784083]]
INFO - 11:02:22: 1%|▏ | 145/10000 [00:00<00:22, 429.76 it/sec, obj=[27.3966129 19.83792999]]
INFO - 11:02:22: 1%|▏ | 146/10000 [00:00<00:22, 430.79 it/sec, obj=[27.3966129 19.83792999]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 8
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22712574556687226
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.47244741e-08 3.68110342e-09]
INFO - 11:02:22: wrapped_ineq1 = [-13.23047421]
INFO - 11:02:22: wrapped_ineq2 = [-34.96615851]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.850561212600975 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.850561110650689 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2271257455668723 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [55.51020408 31.2244898 ]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.850561212600975 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.850561110650689 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2271257455668723 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 1%|▏ | 147/10000 [00:00<00:23, 426.32 it/sec, obj=[27.3966129 19.83792999]]
INFO - 11:02:22: 1%|▏ | 148/10000 [00:00<00:23, 427.39 it/sec, obj=[25.00550716 20.90078046]]
INFO - 11:02:22: 1%|▏ | 149/10000 [00:00<00:22, 428.53 it/sec, obj=[24.98930201 20.90083674]]
INFO - 11:02:22: 2%|▏ | 150/10000 [00:00<00:22, 429.66 it/sec, obj=[24.98737592 20.90043263]]
INFO - 11:02:22: 2%|▏ | 151/10000 [00:00<00:22, 429.07 it/sec, obj=[24.98717553 20.90052426]]
INFO - 11:02:22: 2%|▏ | 152/10000 [00:00<00:22, 429.92 it/sec, obj=[24.98717553 20.90052426]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.22443403356452646
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.44083359e-08 3.60209285e-09]
INFO - 11:02:22: wrapped_ineq1 = [-12.7824248]
INFO - 11:02:22: wrapped_ineq2 = [-35.91369398]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.767313536534702 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.767313426191696 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2244340335645265 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [52.73469388 32.16326531]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.767313536534702 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.767313426191696 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2244340335645265 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 153/10000 [00:00<00:23, 426.98 it/sec, obj=[24.98717553 20.90052426]]
INFO - 11:02:22: 2%|▏ | 154/10000 [00:00<00:23, 427.73 it/sec, obj=[24.98717553 20.90052426]]
INFO - 11:02:22: 2%|▏ | 155/10000 [00:00<00:22, 428.77 it/sec, obj=[22.69510365 21.99680919]]
INFO - 11:02:22: 2%|▏ | 156/10000 [00:00<00:22, 429.85 it/sec, obj=[22.67895762 21.99693969]]
INFO - 11:02:22: 2%|▏ | 157/10000 [00:00<00:22, 430.95 it/sec, obj=[22.67703266 21.9965447 ]]
INFO - 11:02:22: 2%|▏ | 158/10000 [00:00<00:22, 431.78 it/sec, obj=[22.67684157 21.99663879]]
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.2210136200267818
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.40650158e-08 3.51628948e-09]
INFO - 11:02:22: wrapped_ineq1 = [-12.31805248]
INFO - 11:02:22: wrapped_ineq2 = [-36.88017568]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.683628639516674 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.683628520841981 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2210136200267818 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [49.95918367 33.10204082]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.683628639516674 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.683628520841981 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2210136200267818 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 159/10000 [00:00<00:23, 427.42 it/sec, obj=[22.67684157 21.99663879]]
INFO - 11:02:22: 2%|▏ | 160/10000 [00:00<00:22, 428.38 it/sec, obj=[20.48537592 23.12688872]]
INFO - 11:02:22: 2%|▏ | 161/10000 [00:00<00:22, 429.42 it/sec, obj=[20.46929353 23.12709601]]
INFO - 11:02:22: 2%|▏ | 162/10000 [00:00<00:22, 430.43 it/sec, obj=[20.46737079 23.1267108 ]]
INFO - 11:02:22: 2%|▏ | 163/10000 [00:00<00:22, 431.32 it/sec, obj=[20.4671891 23.12680736]]
INFO - 11:02:22: 2%|▏ | 164/10000 [00:00<00:22, 432.27 it/sec, obj=[20.4671891 23.12680736]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.21685290137763716
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.36969263e-08 3.42422446e-09]
INFO - 11:02:22: wrapped_ineq1 = [-11.83709732]
INFO - 11:02:22: wrapped_ineq2 = [-37.86590532]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.59949955938698 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.599499432470258 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2168529013776372 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [47.18367347 34.04081633]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.59949955938698 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.599499432470258 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2168529013776372 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 165/10000 [00:00<00:22, 429.31 it/sec, obj=[20.4671891 23.12680736]]
INFO - 11:02:22: 2%|▏ | 166/10000 [00:00<00:22, 429.92 it/sec, obj=[20.4671891 23.12680736]]
INFO - 11:02:22: 2%|▏ | 167/10000 [00:00<00:22, 430.94 it/sec, obj=[18.37794411 24.29156696]]
INFO - 11:02:22: 2%|▏ | 168/10000 [00:00<00:22, 431.91 it/sec, obj=[18.36193008 24.29185378]]
INFO - 11:02:22: 2%|▏ | 169/10000 [00:00<00:22, 432.92 it/sec, obj=[18.36001075 24.29147902]]
INFO - 11:02:22: 2%|▏ | 170/10000 [00:00<00:22, 433.77 it/sec, obj=[18.35983854 24.29157804]]
INFO - 11:02:22: 2%|▏ | 171/10000 [00:00<00:22, 434.70 it/sec, obj=[18.35983854 24.29157804]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 8
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.21193996278057933
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.33062201e-08 3.32653372e-09]
INFO - 11:02:22: wrapped_ineq1 = [-11.33929244]
INFO - 11:02:22: wrapped_ineq2 = [-38.87119268]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.514919147229693 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.514919012188628 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2119399627805793 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [44.40816327 34.97959184]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.514919147229693 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.514919012188628 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2119399627805793 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 172/10000 [00:00<00:22, 431.96 it/sec, obj=[18.35983854 24.29157804]]
INFO - 11:02:22: 2%|▏ | 173/10000 [00:00<00:22, 432.59 it/sec, obj=[18.35983854 24.29157804]]
INFO - 11:02:22: 2%|▏ | 174/10000 [00:00<00:22, 433.53 it/sec, obj=[16.37447224 25.49140673]]
INFO - 11:02:22: 2%|▏ | 175/10000 [00:00<00:22, 434.44 it/sec, obj=[16.35853156 25.49177588]]
INFO - 11:02:22: 2%|▏ | 176/10000 [00:00<00:22, 435.43 it/sec, obj=[16.35661687 25.4914123 ]]
INFO - 11:02:22: 2%|▏ | 177/10000 [00:00<00:22, 436.17 it/sec, obj=[16.35645426 25.49151379]]
INFO - 11:02:22: 2%|▏ | 178/10000 [00:00<00:22, 437.06 it/sec, obj=[16.35645426 25.49151379]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.2062625663224441
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.28956579e-08 3.22393312e-09]
INFO - 11:02:22: wrapped_ineq1 = [-10.8243637]
INFO - 11:02:22: wrapped_ineq2 = [-39.89635594]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.42988006028335 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.429879917262757 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2062625663224441 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [41.63265306 35.91836735]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.42988006028335 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.429879917262757 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.2062625663224441 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 179/10000 [00:00<00:22, 434.50 it/sec, obj=[16.35645426 25.49151379]]
INFO - 11:02:22: 2%|▏ | 180/10000 [00:00<00:22, 435.03 it/sec, obj=[16.35645426 25.49151379]]
INFO - 11:02:22: 2%|▏ | 181/10000 [00:00<00:22, 435.99 it/sec, obj=[14.47666997 26.72698623]]
INFO - 11:02:22: 2%|▏ | 182/10000 [00:00<00:22, 436.88 it/sec, obj=[14.46080786 26.72744066]]
INFO - 11:02:22: 2%|▏ | 183/10000 [00:00<00:22, 437.77 it/sec, obj=[14.45889916 26.72708902]]
INFO - 11:02:22: 2%|▏ | 184/10000 [00:00<00:22, 438.59 it/sec, obj=[14.45874622 26.72719297]]
INFO - 11:02:22: 2%|▏ | 185/10000 [00:00<00:22, 439.39 it/sec, obj=[14.45874622 26.72719297]]
INFO - 11:02:22: 2%|▏ | 186/10000 [00:00<00:22, 440.16 it/sec, obj=[14.45874622 26.72719297]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.19980813861420837
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.24679964e-08 3.11698045e-09]
INFO - 11:02:22: wrapped_ineq1 = [-10.29202946]
INFO - 11:02:22: wrapped_ineq2 = [-40.94172201]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.344374754501563 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.344374603672217 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1998081386142084 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [38.85714286 36.85714286]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.344374754501563 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.344374603672217 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1998081386142084 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 187/10000 [00:00<00:22, 436.15 it/sec, obj=[14.45874622 26.72719297]]
INFO - 11:02:22: 2%|▏ | 188/10000 [00:00<00:22, 437.01 it/sec, obj=[12.68629432 27.99889971]]
INFO - 11:02:22: 2%|▏ | 189/10000 [00:00<00:22, 437.87 it/sec, obj=[12.67051626 27.99944249]]
INFO - 11:02:22: 2%|▏ | 190/10000 [00:00<00:22, 438.76 it/sec, obj=[12.66861499 27.99910359]]
INFO - 11:02:22: 2%|▏ | 191/10000 [00:00<00:22, 439.49 it/sec, obj=[12.66847181 27.99921 ]]
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.19256375777317306
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.20261276e-08 3.00653902e-09]
INFO - 11:02:22: wrapped_ineq1 = [-9.74200027]
INFO - 11:02:22: wrapped_ineq2 = [-42.00762687]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.258395476742429 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.258395318299658 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1925637577731731 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [36.08163265 37.79591837]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.258395476742429 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.258395318299658 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1925637577731731 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 192/10000 [00:00<00:22, 435.76 it/sec, obj=[12.66847181 27.99921 ]]
INFO - 11:02:22: 2%|▏ | 193/10000 [00:00<00:22, 436.55 it/sec, obj=[11.00515155 29.30775805]]
INFO - 11:02:22: 2%|▏ | 194/10000 [00:00<00:22, 437.46 it/sec, obj=[10.98946331 29.30839237]]
INFO - 11:02:22: 2%|▏ | 195/10000 [00:00<00:22, 438.30 it/sec, obj=[10.98757099 29.30806707]]
INFO - 11:02:22: 2%|▏ | 196/10000 [00:00<00:22, 439.00 it/sec, obj=[10.98743766 29.30817594]]
INFO - 11:02:22: 2%|▏ | 197/10000 [00:00<00:22, 439.78 it/sec, obj=[10.98743766 29.30817594]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.18451613974837508
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.15728334e-08 2.89322699e-09]
INFO - 11:02:22: wrapped_ineq1 = [-9.17397857]
INFO - 11:02:22: wrapped_ineq2 = [-43.09441589]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.171934256563914 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.171934090725966 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1845161397483751 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [33.30612245 38.73469388]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.171934256563914 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.171934090725966 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1845161397483751 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 198/10000 [00:00<00:22, 436.39 it/sec, obj=[10.98743766 29.30817594]]
INFO - 11:02:22: 2%|▏ | 199/10000 [00:00<00:22, 437.25 it/sec, obj=[ 9.43509909 30.65418945]]
INFO - 11:02:22: 2%|▏ | 200/10000 [00:00<00:22, 438.07 it/sec, obj=[ 9.41950673 30.65491867]]
INFO - 11:02:22: 2%|▏ | 201/10000 [00:00<00:22, 438.89 it/sec, obj=[ 9.41762501 30.65460787]]
INFO - 11:02:22: 2%|▏ | 202/10000 [00:00<00:22, 439.36 it/sec, obj=[ 9.4175016 30.65471918]]
INFO - 11:02:22: 2%|▏ | 203/10000 [00:00<00:22, 440.03 it/sec, obj=[ 9.4175016 30.65471918]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.17565162394828981
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.11113643e-08 2.77783130e-09]
INFO - 11:02:22: wrapped_ineq1 = [-8.58765836]
INFO - 11:02:22: wrapped_ineq2 = [-44.20244425]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.084982897600546 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.084982724606786 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1756516239482898 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [30.53061224 39.67346939]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 1.084982897600546 | 5 | float |
INFO - 11:02:22: | y | 0 | 1.084982724606786 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1756516239482898 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 204/10000 [00:00<00:22, 436.56 it/sec, obj=[ 9.4175016 30.65471918]]
INFO - 11:02:22: 2%|▏ | 205/10000 [00:00<00:22, 437.33 it/sec, obj=[ 7.97804763 32.03884013]]
INFO - 11:02:22: 2%|▏ | 206/10000 [00:00<00:22, 438.05 it/sec, obj=[ 7.96255753 32.03966771]]
INFO - 11:02:22: 2%|▏ | 207/10000 [00:00<00:22, 438.82 it/sec, obj=[ 7.96068815 32.03937238]]
INFO - 11:02:22: 2%|▏ | 208/10000 [00:00<00:22, 439.50 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 11:02:22: 2%|▏ | 209/10000 [00:00<00:22, 440.17 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.16595615812678902
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.06442215e-08 2.66106248e-09]
INFO - 11:02:22: wrapped_ineq1 = [-7.98272487]
INFO - 11:02:22: wrapped_ineq2 = [-45.33207732]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.9975329684950597 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.9975327886040215 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.165956158126789 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [27.75510204 40.6122449 ]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.9975329684950597 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.9975327886040215 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.165956158126789 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 210/10000 [00:00<00:22, 436.92 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 11:02:22: 2%|▏ | 211/10000 [00:00<00:22, 437.64 it/sec, obj=[ 6.63596326 33.46237501]]
INFO - 11:02:22: 2%|▏ | 212/10000 [00:00<00:22, 438.44 it/sec, obj=[ 6.62058212 33.46330461]]
INFO - 11:02:22: 2%|▏ | 213/10000 [00:00<00:22, 439.21 it/sec, obj=[ 6.61872696 33.46302576]]
INFO - 11:02:22: 2%|▏ | 214/10000 [00:00<00:22, 439.82 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 11:02:22: 2%|▏ | 215/10000 [00:00<00:22, 440.13 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 11:02:22: 2%|▏ | 216/10000 [00:00<00:22, 440.77 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 10
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.15541528248005268
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [1.01742508e-08 2.54357690e-09]
INFO - 11:02:22: wrapped_ineq1 = [-7.3588542]
INFO - 11:02:22: wrapped_ineq2 = [-46.48369104]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.9095757933566083 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.9095756068439402 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1554152824800527 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [24.97959184 41.55102041]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.9095757933566083 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.9095756068439402 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1554152824800527 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 217/10000 [00:00<00:22, 438.36 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 11:02:22: 2%|▏ | 218/10000 [00:00<00:22, 438.78 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 11:02:22: 2%|▏ | 219/10000 [00:00<00:22, 439.53 it/sec, obj=[ 5.41086973 34.92547856]]
INFO - 11:02:22: 2%|▏ | 220/10000 [00:00<00:22, 440.29 it/sec, obj=[ 5.39560462 34.92651398]]
INFO - 11:02:22: 2%|▏ | 221/10000 [00:00<00:22, 441.02 it/sec, obj=[ 5.39376567 34.92625268]]
INFO - 11:02:22: 2%|▏ | 222/10000 [00:00<00:22, 441.65 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 11:02:22: 2%|▏ | 223/10000 [00:00<00:22, 442.28 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.14401411290347046
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [9.70410241e-09 2.42604159e-09]
INFO - 11:02:22: wrapped_ineq1 = [-6.71571295]
INFO - 11:02:22: wrapped_ineq2 = [-47.6576724]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.8211024417149247 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.8211022488712774 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1440141129034705 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [22.20408163 42.48979592]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.8211024417149247 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.8211022488712774 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1440141129034705 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 224/10000 [00:00<00:22, 440.18 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 11:02:22: 2%|▏ | 225/10000 [00:00<00:22, 440.64 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 11:02:22: 2%|▏ | 226/10000 [00:00<00:22, 441.36 it/sec, obj=[ 4.30485089 36.42885552]]
INFO - 11:02:22: 2%|▏ | 227/10000 [00:00<00:22, 442.14 it/sec, obj=[ 4.28970922 36.43000074]]
INFO - 11:02:22: 2%|▏ | 228/10000 [00:00<00:22, 442.84 it/sec, obj=[ 4.28788864 36.42975814]]
INFO - 11:02:22: 2%|▏ | 229/10000 [00:00<00:22, 443.46 it/sec, obj=[ 4.28780567 36.42987905]]
INFO - 11:02:22: 2%|▏ | 230/10000 [00:00<00:21, 444.09 it/sec, obj=[ 4.28780567 36.42987905]]
ERROR - 11:02:22: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 8
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.13173732335367283
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [9.23635923e-09 2.30908626e-09]
INFO - 11:02:22: wrapped_ineq1 = [-6.05295781]
INFO - 11:02:22: wrapped_ineq2 = [-48.85441989]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.7321037179375566 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.7321035190664417 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1317373233536728 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [19.42857143 43.42857143]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.7321037179375566 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.7321035190664417 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1317373233536728 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 231/10000 [00:00<00:22, 440.88 it/sec, obj=[ 4.28780567 36.42987905]]
INFO - 11:02:22: 2%|▏ | 232/10000 [00:00<00:22, 441.61 it/sec, obj=[ 3.32005317 37.97323185]]
INFO - 11:02:22: 2%|▏ | 233/10000 [00:00<00:22, 442.31 it/sec, obj=[ 3.30504278 37.97449102]]
INFO - 11:02:22: 2%|▏ | 234/10000 [00:00<00:22, 443.04 it/sec, obj=[ 3.30324284 37.97426832]]
INFO - 11:02:22: 2%|▏ | 235/10000 [00:00<00:22, 443.64 it/sec, obj=[ 3.30317013 37.97439158]]
INFO - 11:02:22: 2%|▏ | 236/10000 [00:00<00:21, 444.28 it/sec, obj=[ 3.30317013 37.97439158]]
INFO - 11:02:22: Optimization result:
INFO - 11:02:22: Optimizer info:
INFO - 11:02:22: Status: None
INFO - 11:02:22: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:22: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:22: Solution:
INFO - 11:02:22: The solution is feasible.
INFO - 11:02:22: Objective: -0.1185691272565207
INFO - 11:02:22: Standardized constraints:
INFO - 11:02:22: beta_sub_optim_constraint = [8.77297524e-09 2.19327045e-09]
INFO - 11:02:22: wrapped_ineq1 = [-5.37023516]
INFO - 11:02:22: wrapped_ineq2 = [-50.07434394]
INFO - 11:02:22: Design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.6425701500746512 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.6425699454903127 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1185691272565207 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: Solving mNBI sub-problem for phi_beta = [16.65306122 44.36734694]
INFO - 11:02:22: Optimization problem:
INFO - 11:02:22: minimize -t_extraction
INFO - 11:02:22: with respect to t, x, y
INFO - 11:02:22: subject to constraints:
INFO - 11:02:22: wrapped_ineq1 <= 0.0
INFO - 11:02:22: wrapped_ineq2 <= 0.0
INFO - 11:02:22: beta_sub_optim_constraint <= 0.0
INFO - 11:02:22: over the design space:
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: | x | 0 | 0.6425701500746512 | 5 | float |
INFO - 11:02:22: | y | 0 | 0.6425699454903127 | 3 | float |
INFO - 11:02:22: | t | -inf | 0.1185691272565207 | inf | float |
INFO - 11:02:22: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:22: 2%|▏ | 237/10000 [00:00<00:22, 441.11 it/sec, obj=[ 3.30317013 37.97439158]]
INFO - 11:02:22: 2%|▏ | 238/10000 [00:00<00:22, 441.77 it/sec, obj=[ 2.45868826 39.55935553]]
INFO - 11:02:22: 2%|▏ | 239/10000 [00:00<00:22, 442.47 it/sec, obj=[ 2.44381739 39.56073302]]
INFO - 11:02:22: 2%|▏ | 240/10000 [00:00<00:22, 443.16 it/sec, obj=[ 2.44204056 39.56053151]]
INFO - 11:02:22: 2%|▏ | 241/10000 [00:00<00:21, 443.75 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 11:02:23: 2%|▏ | 242/10000 [00:00<00:21, 444.38 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.10449325789723277
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [8.31663449e-09 2.07913864e-09]
INFO - 11:02:23: wrapped_ineq1 = [-4.66718063]
INFO - 11:02:23: wrapped_ineq2 = [-51.3178675]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.552491978093023 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.552491768118351 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.1044932578972328 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [13.87755102 45.30612245]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.552491978093023 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.552491768118351 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.1044932578972328 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: 2%|▏ | 243/10000 [00:00<00:22, 441.48 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 11:02:23: 2%|▏ | 244/10000 [00:00<00:22, 442.16 it/sec, obj=[ 1.72303596 41.18799761]]
INFO - 11:02:23: 2%|▏ | 245/10000 [00:00<00:22, 442.81 it/sec, obj=[ 1.70831329 41.18949797]]
INFO - 11:02:23: 2%|▏ | 246/10000 [00:00<00:21, 443.49 it/sec, obj=[ 1.70656219 41.18931901]]
INFO - 11:02:23: 2%|▏ | 247/10000 [00:00<00:21, 444.07 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 11:02:23: 2%|▏ | 248/10000 [00:00<00:21, 444.64 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 11:02:23: 2%|▏ | 249/10000 [00:00<00:21, 445.23 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 10
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.08949294772366469
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [7.86864796e-09 1.96713668e-09]
INFO - 11:02:23: wrapped_ineq1 = [-3.94341862]
INFO - 11:02:23: wrapped_ineq2 = [-52.58542653]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.4618591414580615 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.4618589264225835 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.08949294772366469 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [11.10204082 46.24489796]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.4618591414580615 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.4618589264225835 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.08949294772366469 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: 2%|▎ | 250/10000 [00:00<00:22, 442.33 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 11:02:23: 3%|▎ | 251/10000 [00:00<00:22, 442.93 it/sec, obj=[ 1.11544708 42.85995314]]
INFO - 11:02:23: 3%|▎ | 252/10000 [00:00<00:21, 443.58 it/sec, obj=[ 1.10088177 42.86158114]]
INFO - 11:02:23: 3%|▎ | 253/10000 [00:00<00:21, 444.19 it/sec, obj=[ 1.09915921 42.86142618]]
INFO - 11:02:23: 3%|▎ | 254/10000 [00:00<00:21, 444.73 it/sec, obj=[ 1.09911754 42.86155626]]
INFO - 11:02:23: 3%|▎ | 255/10000 [00:00<00:21, 445.33 it/sec, obj=[ 1.09911754 42.86155626]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.07355090648824657
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [7.43066231e-09 1.85767357e-09]
INFO - 11:02:23: wrapped_ineq1 = [-3.19856184]
INFO - 11:02:23: wrapped_ineq2 = [-53.87747056]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.3706612660188607 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.3706610462568307 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.07355090648824657 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [ 8.32653061 47.18367347]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.3706612660188607 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.3706610462568307 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.07355090648824657 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: 3%|▎ | 256/10000 [00:00<00:22, 442.30 it/sec, obj=[ 1.09911754 42.86155626]]
INFO - 11:02:23: 3%|▎ | 257/10000 [00:00<00:21, 442.89 it/sec, obj=[ 0.63834664 44.57604226]]
INFO - 11:02:23: 3%|▎ | 258/10000 [00:00<00:21, 443.54 it/sec, obj=[ 0.62394835 44.57780289]]
INFO - 11:02:23: 3%|▎ | 259/10000 [00:00<00:21, 444.17 it/sec, obj=[ 0.62225734 44.57767347]]
INFO - 11:02:23: 3%|▎ | 260/10000 [00:00<00:21, 444.71 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 11:02:23: 3%|▎ | 261/10000 [00:00<00:21, 445.24 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.05664929814788619
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [7.00442992e-09 1.75109216e-09]
INFO - 11:02:23: wrapped_ineq1 = [-2.43221075]
INFO - 11:02:23: wrapped_ineq2 = [-55.19446335]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.2788876501480995 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.2788874259967131 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.05664929814788619 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [ 5.55102041 48.12244898]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.2788876501480995 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.2788874259967131 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.05664929814788619 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: 3%|▎ | 262/10000 [00:00<00:21, 443.35 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 11:02:23: 3%|▎ | 263/10000 [00:00<00:21, 443.45 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 11:02:23: 3%|▎ | 264/10000 [00:00<00:21, 444.04 it/sec, obj=[ 0.29423714 46.33711132]]
INFO - 11:02:23: 3%|▎ | 265/10000 [00:00<00:21, 444.66 it/sec, obj=[ 0.28001607 46.33900984]]
INFO - 11:02:23: 3%|▎ | 266/10000 [00:00<00:21, 445.29 it/sec, obj=[ 0.27835982 46.3389076 ]]
INFO - 11:02:23: 3%|▎ | 267/10000 [00:00<00:21, 445.84 it/sec, obj=[ 0.27833898 46.33904203]]
INFO - 11:02:23: 3%|▎ | 268/10000 [00:00<00:21, 446.39 it/sec, obj=[ 0.27833898 46.33904203]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 9
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.038769716434479316
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [6.59059385e-09 1.64767378e-09]
INFO - 11:02:23: wrapped_ineq1 = [-1.64395306]
INFO - 11:02:23: wrapped_ineq2 = [-56.53688344]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.1865272500842602 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.1865270218820104 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.03876971643447932 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [ 2.7755102 49.06122449]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.1865272500842602 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.1865270218820104 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.03876971643447932 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: 3%|▎ | 269/10000 [00:00<00:21, 443.71 it/sec, obj=[ 0.27833898 46.33904203]]
INFO - 11:02:23: 3%|▎ | 270/10000 [00:00<00:21, 444.31 it/sec, obj=[ 0.08570212 48.14403408]]
INFO - 11:02:23: 3%|▎ | 271/10000 [00:00<00:21, 444.93 it/sec, obj=[ 0.07166903 48.14607599]]
INFO - 11:02:23: 3%|▎ | 272/10000 [00:00<00:21, 445.53 it/sec, obj=[ 0.07005101 48.14600266]]
INFO - 11:02:23: 3%|▎ | 273/10000 [00:00<00:21, 446.06 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 11:02:23: 3%|▎ | 274/10000 [00:00<00:21, 446.62 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 11:02:23: 3%|▎ | 275/10000 [00:00<00:21, 447.16 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 11:02:23: 3%|▎ | 276/10000 [00:00<00:21, 447.49 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 11
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.01989315900133865
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [6.19085672e-09 1.54771107e-09]
INFO - 11:02:23: wrapped_ineq1 = [-0.83336312]
INFO - 11:02:23: wrapped_ineq2 = [-57.9052249]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0.09356866441943447 | 5 | float |
INFO - 11:02:23: | y | 0 | 0.0935684325046217 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.01989315900133865 | inf | float |
INFO - 11:02:23: +------+-------------+---------------------+-------------+-------+
INFO - 11:02:23: 3%|▎ | 277/10000 [00:00<00:22, 425.24 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: None
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 410
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: 30.28780069872395
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: ineq1 = [-10.82436312]
INFO - 11:02:23: ineq2 = [-39.89635695]
INFO - 11:02:23: Pareto efficient solutions:
INFO - 11:02:23: Pareto optimal points : 203 / 277
INFO - 11:02:23: Utopia point : [0. 4.]
INFO - 11:02:23: Compromise solution (closest to utopia) : [[14.45889916 26.72708902]]
INFO - 11:02:23: Distance from utopia : 26.936598526003834
INFO - 11:02:23: Objective values:
INFO - 11:02:23: +------+----+-----+---------+
INFO - 11:02:23: | name | 1 | 2 | 3 |
INFO - 11:02:23: +------+----+-----+---------+
INFO - 11:02:23: | 1 | 0 | 136 | 14.4589 |
INFO - 11:02:23: | 2 | 50 | 4 | 26.7271 |
INFO - 11:02:23: +------+----+-----+---------+
INFO - 11:02:23: Constraint values:
INFO - 11:02:23: +------+-------+------+----------+
INFO - 11:02:23: | name | 1 | 2 | 3 |
INFO - 11:02:23: +------+-------+------+----------+
INFO - 11:02:23: | 1 | 0 | -22 | -10.2918 |
INFO - 11:02:23: | 2 | -59.3 | -7.3 | -40.9422 |
INFO - 11:02:23: +------+-------+------+----------+
INFO - 11:02:23: Design space:
INFO - 11:02:23: +-------+-------------+---+---+---------+-------------+-------+
INFO - 11:02:23: | name | lower_bound | 1 | 2 | 3 | upper_bound | type |
INFO - 11:02:23: +-------+-------------+---+---+---------+-------------+-------+
INFO - 11:02:23: | x (1) | 0 | 0 | 5 | 1.34434 | 5 | float |
INFO - 11:02:23: | y (1) | 0 | 0 | 3 | 1.34442 | 3 | float |
INFO - 11:02:23: +-------+-------------+---+---+---------+-------------+-------+
Display the Pareto front#
GEMSEO detects the Pareto optimal points and the dominated ones.
execute_post(problem, post_name="ParetoFront", save=False, show=True)
<gemseo.post.pareto_front.ParetoFront object at 0x7f63527f8c50>
Refine the Pareto front in the user specified area#
The Pareto front is then refined with 5 new sub-optimizations.
The custom_anchor_points
argument corresponds to the bounds of
both objectives in order to define the refinement area.
mnbi_settings = MNBI_Settings(
max_iter=10000,
sub_optim_max_iter=200,
n_sub_optim=5,
sub_optim_algo="NLOPT_SLSQP",
custom_anchor_points=[array([44.5, 14]), array([29.4, 19])],
)
execute_algo(problem, settings_model=mnbi_settings)
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize compute_binhkorn(x, y) = (4*x**2+ 4*y**2, (x-5.)**2 + (y-5.)**2)
INFO - 11:02:23: with respect to x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 11:02:23: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 1.429879988721573 | 5 | float |
INFO - 11:02:23: | y | 0 | 1.429879988721184 | 3 | float |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: Solving optimization problem with algorithm MNBI:
WARNING - 11:02:23: Option `custom_anchor_points` was set. The resulting Pareto front might be incomplete.
INFO - 11:02:23: Searching for the individual optimum of each objective
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize f_0
INFO - 11:02:23: with respect to x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 11:02:23: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 1.429879988721573 | 5 | float |
INFO - 11:02:23: | y | 0 | 1.429879988721184 | 3 | float |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
ERROR - 11:02:23: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 6
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: 0.0
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: ineq1 = [0.]
INFO - 11:02:23: ineq2 = [-59.3]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+-------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+-------+-------------+-------+
INFO - 11:02:23: | x | 0 | 0 | 5 | float |
INFO - 11:02:23: | y | 0 | 0 | 3 | float |
INFO - 11:02:23: +------+-------------+-------+-------------+-------+
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize f_1
INFO - 11:02:23: with respect to x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 11:02:23: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 1.429879988721573 | 5 | float |
INFO - 11:02:23: | y | 0 | 1.429879988721184 | 3 | float |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
ERROR - 11:02:23: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.12/site-packages/nlopt/nlopt.py", line 332, in optimize
return _nlopt.opt_optimize(self, *args)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 5
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: 4.0
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: ineq1 = [-22.]
INFO - 11:02:23: ineq2 = [-7.3]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+-------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+-------+-------------+-------+
INFO - 11:02:23: | x | 0 | 5 | 5 | float |
INFO - 11:02:23: | y | 0 | 3 | 3 | float |
INFO - 11:02:23: +------+-------------+-------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [63.01630435 28.68566176]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 1.429879988721573 | 5 | float |
INFO - 11:02:23: | y | 0 | 1.429879988721184 | 3 | float |
INFO - 11:02:23: | t | -inf | 0 | inf | float |
INFO - 11:02:23: +------+-------------+-------------------+-------------+-------+
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 10
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.23006624766887837
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [2.35958401e-05 5.89896004e-06]
INFO - 11:02:23: wrapped_ineq1 = [-13.95717111]
INFO - 11:02:23: wrapped_ineq2 = [-33.39411594]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 1.991452157953277 | 5 | float |
INFO - 11:02:23: | y | 0 | 1.991468771474438 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2300662476688784 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [65.50652174 27.84338235]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 1.991452157953277 | 5 | float |
INFO - 11:02:23: | y | 0 | 1.991468771474438 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2300662476688784 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.23079497225296655
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [3.84877907e-09 9.62177893e-10]
INFO - 11:02:23: wrapped_ineq1 = [-14.32145361]
INFO - 11:02:23: wrapped_ineq2 = [-32.58771023]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.06513935996432 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.065139414923385 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2307949722529666 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [67.99673913 27.00110294]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.06513935996432 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.065139414923385 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2307949722529666 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.23096895870405437
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [3.93921695e-09 9.84808679e-10]
INFO - 11:02:23: wrapped_ineq1 = [-14.67324895]
INFO - 11:02:23: wrapped_ineq2 = [-31.7958384]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.138485442599787 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.1384854919883 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2309689587040544 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [70.48695652 26.15882353]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.138485442599787 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.1384854919883 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2309689587040544 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.2305955098019026
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [4.02403089e-09 1.00601660e-09]
INFO - 11:02:23: wrapped_ineq1 = [-15.01278173]
INFO - 11:02:23: wrapped_ineq2 = [-31.01819989]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.211503064031537 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.211503107833047 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2305955098019026 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Solving mNBI sub-problem for phi_beta = [72.97717391 25.31654412]
INFO - 11:02:23: Optimization problem:
INFO - 11:02:23: minimize -t_extraction
INFO - 11:02:23: with respect to t, x, y
INFO - 11:02:23: subject to constraints:
INFO - 11:02:23: wrapped_ineq1 <= 0.0
INFO - 11:02:23: wrapped_ineq2 <= 0.0
INFO - 11:02:23: beta_sub_optim_constraint <= 0.0
INFO - 11:02:23: over the design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.211503064031537 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.211503107833047 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2305955098019026 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 7
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: -0.22968191496902265
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: beta_sub_optim_constraint = [4.09995593e-09 1.02497921e-09]
INFO - 11:02:23: wrapped_ineq1 = [-15.34021524]
INFO - 11:02:23: wrapped_ineq2 = [-30.25460517]
INFO - 11:02:23: Design space:
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | Name | Lower bound | Value | Upper bound | Type |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: | x | 0 | 2.284196597883743 | 5 | float |
INFO - 11:02:23: | y | 0 | 2.284196636140467 | 3 | float |
INFO - 11:02:23: | t | -inf | 0.2296819149690227 | inf | float |
INFO - 11:02:23: +------+-------------+--------------------+-------------+-------+
INFO - 11:02:23: Optimization result:
INFO - 11:02:23: Optimizer info:
INFO - 11:02:23: Status: None
INFO - 11:02:23: Message: None
INFO - 11:02:23: Number of calls to the objective function by the optimizer: 464
INFO - 11:02:23: Solution:
INFO - 11:02:23: The solution is feasible.
INFO - 11:02:23: Objective: 30.28780069872395
INFO - 11:02:23: Standardized constraints:
INFO - 11:02:23: ineq1 = [-10.82436312]
INFO - 11:02:23: ineq2 = [-39.89635695]
INFO - 11:02:23: Pareto efficient solutions:
INFO - 11:02:23: Pareto optimal points : 226 / 312
INFO - 11:02:23: Utopia point : [0. 4.]
INFO - 11:02:23: Compromise solution (closest to utopia) : [[14.45889916 26.72708902]]
INFO - 11:02:23: Distance from utopia : 26.936598526003834
INFO - 11:02:23: Objective values:
INFO - 11:02:23: +------+----+-----+---------+
INFO - 11:02:23: | name | 1 | 2 | 3 |
INFO - 11:02:23: +------+----+-----+---------+
INFO - 11:02:23: | 1 | 0 | 136 | 14.4589 |
INFO - 11:02:23: | 2 | 50 | 4 | 26.7271 |
INFO - 11:02:23: +------+----+-----+---------+
INFO - 11:02:23: Constraint values:
INFO - 11:02:23: +------+-------+------+----------+
INFO - 11:02:23: | name | 1 | 2 | 3 |
INFO - 11:02:23: +------+-------+------+----------+
INFO - 11:02:23: | 1 | 0 | -22 | -10.2918 |
INFO - 11:02:23: | 2 | -59.3 | -7.3 | -40.9422 |
INFO - 11:02:23: +------+-------+------+----------+
INFO - 11:02:23: Design space:
INFO - 11:02:23: +-------+-------------+---+---+---------+-------------+-------+
INFO - 11:02:23: | name | lower_bound | 1 | 2 | 3 | upper_bound | type |
INFO - 11:02:23: +-------+-------------+---+---+---------+-------------+-------+
INFO - 11:02:23: | x (1) | 0 | 0 | 5 | 1.34434 | 5 | float |
INFO - 11:02:23: | y (1) | 0 | 0 | 3 | 1.34442 | 3 | float |
INFO - 11:02:23: +-------+-------------+---+---+---------+-------------+-------+
Display the Pareto front#
We can clearly see the effect of the local refinement.
execute_post(problem, post_name="ParetoFront", save=False, show=True)
<gemseo.post.pareto_front.ParetoFront object at 0x7f6361b854f0>
Total running time of the script: (0 minutes 1.221 seconds)