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 gemseo import configure_logger
from gemseo import execute_algo
from gemseo import execute_post
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()
result = execute_algo(
problem,
"MNBI",
max_iter=10000,
sub_optim_max_iter=200,
n_sub_optim=50,
sub_optim_algo="NLOPT_SLSQP",
)
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize compute_binhkorn(x, y) = (4*x**2+ 4*y**2, (x-5.)**2 + (y-5.)**2)
INFO - 10:40:03: with respect to x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 10:40:03: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | x | 0 | 1 | 5 | float |
INFO - 10:40:03: | y | 0 | 1 | 3 | float |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: Solving optimization problem with algorithm MNBI:
INFO - 10:40:03: Searching for the individual optimum of each objective
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize f_0
INFO - 10:40:03: with respect to x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 10:40:03: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | x | 0 | 1 | 5 | float |
INFO - 10:40:03: | y | 0 | 1 | 3 | float |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 4
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: 2.6137650460903986e-22
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: ineq1 = [3.61310981e-12]
INFO - 10:40:03: ineq2 = [-59.3]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+-----------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-----------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.459210778110538e-13 | 5 | float |
INFO - 10:40:03: | y | 0 | 8.071262636330439e-12 | 3 | float |
INFO - 10:40:03: +------+-------------+-----------------------+-------------+-------+
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize f_1
INFO - 10:40:03: with respect to x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 10:40:03: ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | x | 0 | 1 | 5 | float |
INFO - 10:40:03: | y | 0 | 1 | 3 | float |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
ERROR - 10:40:03: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 4
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: 4.0
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: ineq1 = [-22.]
INFO - 10:40:03: ineq2 = [-7.3]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.999999999999934 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [133.2244898 4.93877551]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: | x | 0 | 1 | 5 | float |
INFO - 10:40:03: | y | 0 | 1 | 3 | float |
INFO - 10:40:03: | t | -inf | 0 | inf | float |
INFO - 10:40:03: +------+-------------+-------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.019986438123775534
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [8.48413073e-09 2.12111273e-09]
INFO - 10:40:03: wrapped_ineq1 = [-21.98060064]
INFO - 10:40:03: wrapped_ineq2 = [-8.15508884]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.860718419626017 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999998 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.01998643812377553 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [130.44897959 5.87755102]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.860718419626017 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999998 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.01998643812377553 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.03911696914200083
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.07045384e-06 2.67613419e-07]
INFO - 10:40:03: wrapped_ineq1 = [-21.92182929]
INFO - 10:40:03: wrapped_ineq2 = [-9.0557122]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.720409750725461 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.03911696914200083 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [127.67346939 6.81632653]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.720409750725461 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.03911696914200083 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.05737236260338867
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.11930945e-06 2.79827344e-07]
INFO - 10:40:03: wrapped_ineq1 = [-21.82280187]
INFO - 10:40:03: wrapped_ineq2 = [-10.00289263]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.579050916756698 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.05737236260338867 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [124.89795918 7.75510204]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.579050916756698 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.05737236260338867 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.07473267476171336
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.1711835e-06 2.9279584e-07]
INFO - 10:40:03: wrapped_ineq1 = [-21.68260071]
INFO - 10:40:03: wrapped_ineq2 = [-10.99769138]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.436617985224268 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.07473267476171336 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [122.12244898 8.69387755]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.436617985224268 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.07473267476171336 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.09117718665694834
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.22631330e-06 3.06578367e-07]
INFO - 10:40:03: wrapped_ineq1 = [-21.50027273]
INFO - 10:40:03: wrapped_ineq2 = [-12.04121071]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.29308609337613 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.09117718665694834 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [119.34693878 9.63265306]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.29308609337613 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.09117718665694834 | inf | float |
INFO - 10:40:03: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.10668436745320037
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.28495816e-06 3.21239511e-07]
INFO - 10:40:03: wrapped_ineq1 = [-21.27482752]
INFO - 10:40:03: wrapped_ineq2 = [-13.13459605]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.148429404209233 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1066843674532004 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [116.57142857 10.57142857]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.148429404209233 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1066843674532004 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.12123182918919541
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.34740281e-06 3.36850676e-07]
INFO - 10:40:03: wrapped_ineq1 = [-21.00523523]
INFO - 10:40:03: wrapped_ineq2 = [-14.27903845]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.002621052170031 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1212318291891954 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [113.79591837 11.51020408]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 4.002621052170031 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1212318291891954 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.1347962782336287
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.41395994e-06 3.53489989e-07]
INFO - 10:40:03: wrapped_ineq1 = [-20.69042436]
INFO - 10:40:03: wrapped_ineq2 = [-15.47577713]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.85563308490096 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1347962782336287 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [111.02040816 12.44897959]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.85563308490096 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1347962782336287 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.14735346314733033
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.48497396e-06 3.71243503e-07]
INFO - 10:40:03: wrapped_ineq1 = [-20.32927934]
INFO - 10:40:03: wrapped_ineq2 = [-16.72610225]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.707436400675021 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1473534631473303 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [108.24489796 13.3877551 ]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.707436400675021 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999999 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1473534631473303 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.15887811862169762
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.56082422e-06 3.90206059e-07]
INFO - 10:40:03: wrapped_ineq1 = [-19.92063796]
INFO - 10:40:03: wrapped_ineq2 = [-18.03135795]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.558000681121891 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1588781186216976 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [105.46938776 14.32653061]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.558000681121891 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1588781186216976 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 51/10000 [00:00<01:06, 148.83 it/sec, obj=[82.4386183 6.53671139]]
INFO - 10:40:03: 1%| | 52/10000 [00:00<01:06, 150.68 it/sec, obj=[82.43861798 6.53671142]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 6
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.1693439051263492
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.64192933e-06 4.10482344e-07]
INFO - 10:40:03: wrapped_ineq1 = [-19.46328861]
INFO - 10:40:03: wrapped_ineq2 = [-19.39294547]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.407294318805109 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1693439051263492 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [102.69387755 15.26530612]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.407294318805109 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1693439051263492 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 53/10000 [00:00<01:08, 145.21 it/sec, obj=[82.43861798 6.53671142]]
INFO - 10:40:03: 1%| | 54/10000 [00:00<01:07, 146.79 it/sec, obj=[82.4386183 6.53671139]]
INFO - 10:40:03: 1%| | 55/10000 [00:00<01:06, 148.57 it/sec, obj=[78.4046267 7.04173918]]
INFO - 10:40:03: 1%| | 56/10000 [00:00<01:06, 150.28 it/sec, obj=[78.38750452 7.04403274]]
ERROR - 10:40:03: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 5
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.17872334385765837
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.72875256e-06 4.32188164e-07]
INFO - 10:40:03: wrapped_ineq1 = [-18.95596726]
INFO - 10:40:03: wrapped_ineq2 = [-20.8123267]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.255284339160327 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999998 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1787233438576584 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [99.91836735 16.20408163]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.255284339160327 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.999999999999998 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.1787233438576584 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 57/10000 [00:00<01:08, 145.07 it/sec, obj=[78.38750419 7.04403278]]
INFO - 10:40:03: 1%| | 58/10000 [00:00<01:07, 146.72 it/sec, obj=[74.50478473 7.6000843 ]]
INFO - 10:40:03: 1%| | 59/10000 [00:00<01:07, 148.32 it/sec, obj=[74.48803564 7.60264575]]
INFO - 10:40:03: 1%| | 60/10000 [00:00<01:06, 149.94 it/sec, obj=[74.48803531 7.6026458 ]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.186987746533081
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.82180677e-06 4.55451675e-07]
INFO - 10:40:03: wrapped_ineq1 = [-18.39735425]
INFO - 10:40:03: wrapped_ineq2 = [-22.29102785]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.101936316248493 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.186987746533081 | inf | float |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [97.14285714 17.14285714]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 3.101936316248493 | 5 | float |
INFO - 10:40:03: | y | 0 | 3 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.186987746533081 | inf | float |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: 1%| | 61/10000 [00:00<01:08, 145.12 it/sec, obj=[74.48803531 7.6026458 ]]
INFO - 10:40:03: 1%| | 62/10000 [00:00<01:07, 146.68 it/sec, obj=[70.75658544 8.21055954]]
INFO - 10:40:03: 1%| | 63/10000 [00:00<01:07, 148.27 it/sec, obj=[70.74019853 8.21247029]]
INFO - 10:40:03: 1%| | 64/10000 [00:00<01:06, 149.81 it/sec, obj=[70.7397555 8.21235896]]
INFO - 10:40:03: 1%| | 65/10000 [00:00<01:05, 151.29 it/sec, obj=[70.73968192 8.21237138]]
INFO - 10:40:03: 1%| | 66/10000 [00:00<01:05, 152.78 it/sec, obj=[70.73968188 8.21237139]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.19414099622659098
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [2.67250751e-07 6.68126585e-08]
INFO - 10:40:03: wrapped_ineq1 = [-17.92018718]
INFO - 10:40:03: wrapped_ineq2 = [-23.5380477]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.973627520197615 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.973627389560757 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.194140996226591 | inf | float |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [94.36734694 18.08163265]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.973627520197615 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.973627389560757 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.194140996226591 | inf | float |
INFO - 10:40:03: +------+-------------+-------------------+-------------+-------+
INFO - 10:40:03: 1%| | 67/10000 [00:00<01:07, 148.14 it/sec, obj=[70.73968188 8.21237139]]
INFO - 10:40:03: 1%| | 68/10000 [00:00<01:06, 149.57 it/sec, obj=[67.10688455 8.85599137]]
INFO - 10:40:03: 1%| | 69/10000 [00:00<01:05, 151.02 it/sec, obj=[67.09022377 8.85523197]]
INFO - 10:40:03: 1%| | 70/10000 [00:00<01:05, 152.44 it/sec, obj=[67.08836417 8.8547516 ]]
INFO - 10:40:03: 1%| | 71/10000 [00:00<01:04, 153.81 it/sec, obj=[67.0880443 8.8548097]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.20058310774563637
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.59261901e-08 3.98151023e-09]
INFO - 10:40:03: wrapped_ineq1 = [-17.67673508]
INFO - 10:40:03: wrapped_ineq2 = [-24.24810451]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.89586006851919 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.89586006893971 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2005831077456364 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [91.59183673 19.02040816]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.89586006851919 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.89586006893971 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2005831077456364 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 72/10000 [00:00<01:06, 149.22 it/sec, obj=[67.0880443 8.8548097]]
INFO - 10:40:03: 1%| | 73/10000 [00:00<01:05, 150.53 it/sec, obj=[63.53594358 9.5256767 ]]
INFO - 10:40:03: 1%| | 74/10000 [00:00<01:05, 151.84 it/sec, obj=[63.51929686 9.52496401]]
INFO - 10:40:03: 1%| | 75/10000 [00:00<01:04, 153.18 it/sec, obj=[63.51742882 9.52448661]]
INFO - 10:40:03: 1%| | 76/10000 [00:00<01:04, 154.45 it/sec, obj=[63.51711681 9.52454701]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.20643176428302623
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.60213460e-08 4.00534361e-09]
INFO - 10:40:03: wrapped_ineq1 = [-17.4199899]
INFO - 10:40:03: wrapped_ineq2 = [-24.97359042]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.817736612798519 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.817736606382921 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2064317642830262 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [88.81632653 19.95918367]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.817736612798519 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.817736606382921 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2064317642830262 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 77/10000 [00:00<01:05, 150.71 it/sec, obj=[63.51711681 9.52454701]]
WARNING - 10:40:03: Optimization found no feasible point; the least infeasible point is selected.
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 1
INFO - 10:40:03: Solution:
WARNING - 10:40:03: The solution is not feasible.
INFO - 10:40:03: Objective: -0.20643176428302623
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = None
INFO - 10:40:03: wrapped_ineq1 = None
INFO - 10:40:03: wrapped_ineq2 = None
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.817736612798519 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.817736606382921 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2064317642830262 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
WARNING - 10:40:03: No feasible optimum has been found for phi_beta = [88.81632653 19.95918367]
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [86.04081633 20.89795918]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.817736612798519 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.817736606382921 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2064317642830262 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 78/10000 [00:00<01:07, 146.76 it/sec, obj=[63.51711681 9.52454701]]
INFO - 10:40:03: 1%| | 79/10000 [00:00<01:07, 147.96 it/sec, obj=[56.69636778 10.95148504]]
INFO - 10:40:03: 1%| | 80/10000 [00:00<01:06, 149.19 it/sec, obj=[56.63052837 10.94904598]]
INFO - 10:40:03: 1%| | 81/10000 [00:00<01:05, 150.42 it/sec, obj=[56.62304869 10.94718774]]
INFO - 10:40:03: 1%| | 82/10000 [00:00<01:05, 151.58 it/sec, obj=[56.6218939 10.94744123]]
INFO - 10:40:03: 1%| | 83/10000 [00:00<01:04, 152.74 it/sec, obj=[56.62189377 10.94744126]]
INFO - 10:40:03: 1%| | 84/10000 [00:00<01:04, 153.91 it/sec, obj=[56.62189377 10.94744126]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.21631561210408246
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [8.23538862e-07 2.05884694e-07]
INFO - 10:40:03: wrapped_ineq1 = [-16.86587948]
INFO - 10:40:03: wrapped_ineq2 = [-26.47170905]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.660401912334083 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.660401311987996 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2163156121040825 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [83.26530612 21.83673469]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.660401912334083 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.660401311987996 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2163156121040825 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 85/10000 [00:00<01:06, 149.88 it/sec, obj=[56.62189377 10.94744126]]
INFO - 10:40:03: 1%| | 86/10000 [00:00<01:05, 151.01 it/sec, obj=[53.31868376 11.70234425]]
INFO - 10:40:03: 1%| | 87/10000 [00:00<01:05, 152.15 it/sec, obj=[53.30209371 11.70178107]]
INFO - 10:40:03: 1%| | 88/10000 [00:00<01:04, 153.28 it/sec, obj=[53.30020299 11.70131446]]
INFO - 10:40:03: 1%| | 89/10000 [00:00<01:04, 154.37 it/sec, obj=[53.29991518 11.70138187]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.2203337570364798
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.61359033e-08 4.03399092e-09]
INFO - 10:40:03: wrapped_ineq1 = [-16.5681293]
INFO - 10:40:03: wrapped_ineq2 = [-27.24479154]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.58117986020856 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.581179831910084 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2203337570364798 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [80.48979592 22.7755102 ]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.58117986020856 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.581179831910084 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2203337570364798 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 90/10000 [00:00<01:05, 150.40 it/sec, obj=[53.29991518 11.70138187]]
INFO - 10:40:03: 1%| | 91/10000 [00:00<01:05, 151.47 it/sec, obj=[50.08202514 12.48509203]]
INFO - 10:40:03: 1%| | 92/10000 [00:00<01:04, 152.54 it/sec, obj=[50.06545915 12.48458197]]
INFO - 10:40:03: 1%| | 93/10000 [00:00<01:04, 153.62 it/sec, obj=[50.06356185 12.48411964]]
INFO - 10:40:03: 1%| | 94/10000 [00:00<01:04, 154.64 it/sec, obj=[50.06328232 12.48418943]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.2237243648301845
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.61119331e-08 4.02797085e-09]
INFO - 10:40:03: wrapped_ineq1 = [-16.25632384]
INFO - 10:40:03: wrapped_ineq2 = [-28.03418671]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.50158157567626 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.50158153967253 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2237243648301845 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [77.71428571 23.71428571]
INFO - 10:40:03: Optimization problem:
INFO - 10:40:03: minimize -t_extraction
INFO - 10:40:03: with respect to t, x, y
INFO - 10:40:03: subject to constraints:
INFO - 10:40:03: wrapped_ineq1 <= 0.0
INFO - 10:40:03: wrapped_ineq2 <= 0.0
INFO - 10:40:03: beta_sub_optim_constraint <= 0.0
INFO - 10:40:03: over the design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.50158157567626 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.50158153967253 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2237243648301845 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: 1%| | 95/10000 [00:00<01:05, 150.97 it/sec, obj=[50.06328232 12.48418943]]
INFO - 10:40:03: 1%| | 96/10000 [00:00<01:05, 151.96 it/sec, obj=[46.93193733 13.29712012]]
INFO - 10:40:03: 1%| | 97/10000 [00:00<01:04, 152.98 it/sec, obj=[46.91539824 13.29666496]]
INFO - 10:40:03: 1%| | 98/10000 [00:00<01:04, 154.01 it/sec, obj=[46.91349492 13.2962073 ]]
INFO - 10:40:03: 1%| | 99/10000 [00:00<01:03, 154.98 it/sec, obj=[46.91322377 13.29627947]]
INFO - 10:40:03: Optimization result:
INFO - 10:40:03: Optimizer info:
INFO - 10:40:03: Status: None
INFO - 10:40:03: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:03: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:03: Solution:
INFO - 10:40:03: The solution is feasible.
INFO - 10:40:03: Objective: -0.22647839676414017
INFO - 10:40:03: Standardized constraints:
INFO - 10:40:03: beta_sub_optim_constraint = [1.60554023e-08 4.01385591e-09]
INFO - 10:40:03: wrapped_ineq1 = [-15.93025908]
INFO - 10:40:03: wrapped_ineq2 = [-28.84013285]
INFO - 10:40:03: Design space:
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: | x | 0 | 2.42160134571469 | 5 | float |
INFO - 10:40:03: | y | 0 | 2.421601301832548 | 3 | float |
INFO - 10:40:03: | t | -inf | 0.2264783967641402 | inf | float |
INFO - 10:40:03: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:03: Solving mNBI sub-problem for phi_beta = [74.93877551 24.65306122]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.42160134571469 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.421601301832548 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2264783967641402 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%| | 100/10000 [00:00<01:05, 151.66 it/sec, obj=[46.91322377 13.29627947]]
INFO - 10:40:04: 1%| | 101/10000 [00:00<01:04, 152.52 it/sec, obj=[46.91322377 13.29627947]]
INFO - 10:40:04: 1%| | 102/10000 [00:00<01:04, 153.51 it/sec, obj=[43.86967855 14.13885408]]
INFO - 10:40:04: 1%| | 103/10000 [00:00<01:04, 154.48 it/sec, obj=[43.8531693 14.13845566]]
INFO - 10:40:04: 1%| | 104/10000 [00:00<01:03, 155.46 it/sec, obj=[43.85126056 14.13800305]]
INFO - 10:40:04: 1%| | 105/10000 [00:00<01:03, 156.38 it/sec, obj=[43.85099789 14.13807762]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.22858660027257494
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.59653339e-08 3.99134237e-09]
INFO - 10:40:04: wrapped_ineq1 = [-15.58972776]
INFO - 10:40:04: wrapped_ineq2 = [-29.66287052]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.341233618774564 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.341233566860503 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2285866002725749 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [72.16326531 25.59183673]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.341233618774564 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.341233566860503 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2285866002725749 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%| | 106/10000 [00:00<01:04, 153.02 it/sec, obj=[43.85099789 14.13807762]]
INFO - 10:40:04: 1%| | 107/10000 [00:00<01:04, 153.93 it/sec, obj=[40.89653773 15.01072985]]
INFO - 10:40:04: 1%| | 108/10000 [00:00<01:03, 154.88 it/sec, obj=[40.88006141 15.01039006]]
INFO - 10:40:04: 1%| | 109/10000 [00:00<01:03, 155.81 it/sec, obj=[40.87814788 15.00994293]]
INFO - 10:40:04: 1%| | 110/10000 [00:00<01:03, 156.67 it/sec, obj=[40.8778938 15.0100199]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.23003949645199573
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.58414615e-08 3.96036448e-09]
INFO - 10:40:04: wrapped_ineq1 = [-15.23451757]
INFO - 10:40:04: wrapped_ineq2 = [-30.50264619]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.260472707502382 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.260472647425706 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2300394964519957 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [69.3877551 26.53061224]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.260472707502382 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.260472647425706 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2300394964519957 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%| | 111/10000 [00:00<01:04, 153.89 it/sec, obj=[40.8778938 15.0100199]]
INFO - 10:40:04: 1%| | 112/10000 [00:00<01:03, 154.67 it/sec, obj=[40.8778938 15.0100199]]
INFO - 10:40:04: 1%| | 113/10000 [00:00<01:03, 155.56 it/sec, obj=[38.01383568 15.91319413]]
INFO - 10:40:04: 1%| | 114/10000 [00:00<01:03, 156.45 it/sec, obj=[37.99739551 15.91291495]]
INFO - 10:40:04: 1%| | 115/10000 [00:00<01:02, 157.34 it/sec, obj=[37.99547788 15.91247373]]
INFO - 10:40:04: 1%| | 116/10000 [00:00<01:02, 158.18 it/sec, obj=[37.99523249 15.91255313]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.2308273722331007
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.56837210e-08 3.92095956e-09]
INFO - 10:40:04: wrapped_ineq1 = [-14.86441091]
INFO - 10:40:04: wrapped_ineq2 = [-31.35971238]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.17931278404434 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.179312715698789 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2308273722331007 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [66.6122449 27.46938776]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.17931278404434 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.179312715698789 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2308273722331007 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%| | 117/10000 [00:00<01:03, 155.07 it/sec, obj=[37.99523249 15.91255313]]
INFO - 10:40:04: 1%| | 118/10000 [00:00<01:03, 155.91 it/sec, obj=[35.22292611 16.84670478]]
INFO - 10:40:04: 1%| | 119/10000 [00:00<01:03, 156.75 it/sec, obj=[35.20652549 16.84648825]]
INFO - 10:40:04: 1%| | 120/10000 [00:00<01:02, 157.61 it/sec, obj=[35.20460449 16.84605342]]
INFO - 10:40:04: 1%| | 121/10000 [00:00<01:02, 158.41 it/sec, obj=[35.20436789 16.84613524]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.23094027220111932
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.54923541e-08 3.87308496e-09]
INFO - 10:40:04: wrapped_ineq1 = [-14.47918481]
INFO - 10:40:04: wrapped_ineq2 = [-32.23432794]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.097747875138651 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.09774779844383 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2309402722011193 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [63.83673469 28.40816327]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.097747875138651 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.09774779844383 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2309402722011193 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%| | 122/10000 [00:00<01:03, 155.40 it/sec, obj=[35.20436789 16.84613524]]
INFO - 10:40:04: 1%| | 123/10000 [00:00<01:03, 156.19 it/sec, obj=[32.52519687 17.81173117]]
INFO - 10:40:04: 1%| | 124/10000 [00:00<01:02, 157.01 it/sec, obj=[32.50883935 17.81157943]]
INFO - 10:40:04: 1%|▏ | 125/10000 [00:00<01:02, 157.81 it/sec, obj=[32.50691576 17.81115147]]
INFO - 10:40:04: 1%|▏ | 126/10000 [00:00<01:02, 158.58 it/sec, obj=[32.50668806 17.81123573]]
INFO - 10:40:04: 1%|▏ | 127/10000 [00:00<01:01, 159.35 it/sec, obj=[32.50668806 17.81123573]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.23036799004509081
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.52680286e-08 3.81700360e-09]
INFO - 10:40:04: wrapped_ineq1 = [-14.07861062]
INFO - 10:40:04: wrapped_ineq2 = [-33.12675824]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.015771856985211 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.015771771887767 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2303679900450908 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [61.06122449 29.34693878]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 2.015771856985211 | 5 | float |
INFO - 10:40:04: | y | 0 | 2.015771771887767 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2303679900450908 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%|▏ | 128/10000 [00:00<01:03, 156.39 it/sec, obj=[32.50668806 17.81123573]]
INFO - 10:40:04: 1%|▏ | 129/10000 [00:00<01:02, 157.15 it/sec, obj=[29.92207113 18.80875463]]
INFO - 10:40:04: 1%|▏ | 130/10000 [00:00<01:02, 157.93 it/sec, obj=[29.90576043 18.8086699 ]]
INFO - 10:40:04: 1%|▏ | 131/10000 [00:00<01:02, 158.71 it/sec, obj=[29.90383509 18.80824933]]
INFO - 10:40:04: 1%|▏ | 132/10000 [00:00<01:01, 159.45 it/sec, obj=[29.9036164 18.80833604]]
INFO - 10:40:04: 1%|▏ | 133/10000 [00:00<01:01, 160.18 it/sec, obj=[29.9036164 18.80833604]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.22910005961528188
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.50115262e-08 3.75288423e-09]
INFO - 10:40:04: wrapped_ineq1 = [-13.66245391]
INFO - 10:40:04: wrapped_ineq2 = [-34.03727539]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.933378449880203 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.933378356354741 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2291000596152819 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [58.28571429 30.28571429]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.933378449880203 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.933378356354741 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2291000596152819 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%|▏ | 134/10000 [00:00<01:02, 157.64 it/sec, obj=[29.9036164 18.80833604]]
INFO - 10:40:04: 1%|▏ | 135/10000 [00:00<01:02, 158.29 it/sec, obj=[29.9036164 18.80833604]]
INFO - 10:40:04: 1%|▏ | 136/10000 [00:00<01:02, 159.03 it/sec, obj=[27.41500866 19.83826888]]
INFO - 10:40:04: 1%|▏ | 137/10000 [00:00<01:01, 159.78 it/sec, obj=[27.39874869 19.83825346]]
INFO - 10:40:04: 1%|▏ | 138/10000 [00:00<01:01, 160.52 it/sec, obj=[27.39682249 19.83784083]]
INFO - 10:40:04: 1%|▏ | 139/10000 [00:00<01:01, 161.22 it/sec, obj=[27.3966129 19.83792999]]
INFO - 10:40:04: 1%|▏ | 140/10000 [00:00<01:00, 161.93 it/sec, obj=[27.3966129 19.83792999]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.22712574556659043
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.47245132e-08 3.68113184e-09]
INFO - 10:40:04: wrapped_ineq1 = [-13.23047421]
INFO - 10:40:04: wrapped_ineq2 = [-34.96615851]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.850561212602271 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.850561110651984 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2271257455665904 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [55.51020408 31.2244898 ]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.850561212602271 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.850561110651984 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2271257455665904 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%|▏ | 141/10000 [00:00<01:01, 159.49 it/sec, obj=[27.3966129 19.83792999]]
INFO - 10:40:04: 1%|▏ | 142/10000 [00:00<01:01, 160.10 it/sec, obj=[27.3966129 19.83792999]]
INFO - 10:40:04: 1%|▏ | 143/10000 [00:00<01:01, 160.82 it/sec, obj=[25.00550716 20.90078046]]
INFO - 10:40:04: 1%|▏ | 144/10000 [00:00<01:01, 161.53 it/sec, obj=[24.98930201 20.90083674]]
INFO - 10:40:04: 1%|▏ | 145/10000 [00:00<01:00, 162.24 it/sec, obj=[24.98737592 20.90043263]]
INFO - 10:40:04: 1%|▏ | 146/10000 [00:00<01:00, 162.91 it/sec, obj=[24.98717553 20.90052426]]
INFO - 10:40:04: 1%|▏ | 147/10000 [00:00<01:00, 163.59 it/sec, obj=[24.98717553 20.90052426]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.22443403356423705
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.44083181e-08 3.60208396e-09]
INFO - 10:40:04: wrapped_ineq1 = [-12.7824248]
INFO - 10:40:04: wrapped_ineq2 = [-35.91369398]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.767313536536093 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.767313426193087 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2244340335642371 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [52.73469388 32.16326531]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.767313536536093 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.767313426193087 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2244340335642371 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 1%|▏ | 148/10000 [00:00<01:01, 161.26 it/sec, obj=[24.98717553 20.90052426]]
INFO - 10:40:04: 1%|▏ | 149/10000 [00:00<01:00, 161.86 it/sec, obj=[24.98717553 20.90052426]]
INFO - 10:40:04: 2%|▏ | 150/10000 [00:00<01:00, 162.54 it/sec, obj=[22.69510365 21.99680919]]
INFO - 10:40:04: 2%|▏ | 151/10000 [00:00<01:00, 163.22 it/sec, obj=[22.67895762 21.99693969]]
INFO - 10:40:04: 2%|▏ | 152/10000 [00:00<01:00, 163.90 it/sec, obj=[22.67703266 21.9965447 ]]
INFO - 10:40:04: 2%|▏ | 153/10000 [00:00<00:59, 164.55 it/sec, obj=[22.67684157 21.99663879]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.2210136200264868
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.40650585e-08 3.51626817e-09]
INFO - 10:40:04: wrapped_ineq1 = [-12.31805248]
INFO - 10:40:04: wrapped_ineq2 = [-36.88017568]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.683628639518165 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.683628520843471 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2210136200264868 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [49.95918367 33.10204082]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.683628639518165 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.683628520843471 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2210136200264868 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 154/10000 [00:00<01:00, 162.24 it/sec, obj=[22.67684157 21.99663879]]
INFO - 10:40:04: 2%|▏ | 155/10000 [00:00<01:00, 162.80 it/sec, obj=[22.67684157 21.99663879]]
INFO - 10:40:04: 2%|▏ | 156/10000 [00:00<01:00, 163.44 it/sec, obj=[20.48537592 23.12688872]]
INFO - 10:40:04: 2%|▏ | 157/10000 [00:00<00:59, 164.09 it/sec, obj=[20.46929353 23.12709601]]
INFO - 10:40:04: 2%|▏ | 158/10000 [00:00<00:59, 164.74 it/sec, obj=[20.46737079 23.1267108 ]]
INFO - 10:40:04: 2%|▏ | 159/10000 [00:00<00:59, 165.30 it/sec, obj=[20.4671891 23.12680736]]
INFO - 10:40:04: 2%|▏ | 160/10000 [00:00<00:59, 165.92 it/sec, obj=[20.4671891 23.12680736]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.21685290137733731
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.36968410e-08 3.42421558e-09]
INFO - 10:40:04: wrapped_ineq1 = [-11.83709732]
INFO - 10:40:04: wrapped_ineq2 = [-37.86590532]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.59949955938857 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.599499432471848 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2168529013773373 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [47.18367347 34.04081633]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.59949955938857 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.599499432471848 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2168529013773373 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 161/10000 [00:00<01:00, 163.25 it/sec, obj=[20.4671891 23.12680736]]
INFO - 10:40:04: 2%|▏ | 162/10000 [00:00<01:00, 163.88 it/sec, obj=[18.37794411 24.29156696]]
INFO - 10:40:04: 2%|▏ | 163/10000 [00:00<00:59, 164.46 it/sec, obj=[18.36193008 24.29185378]]
INFO - 10:40:04: 2%|▏ | 164/10000 [00:00<00:59, 165.10 it/sec, obj=[18.36001075 24.29147902]]
INFO - 10:40:04: 2%|▏ | 165/10000 [00:00<00:59, 165.69 it/sec, obj=[18.35983854 24.29157804]]
INFO - 10:40:04: 2%|▏ | 166/10000 [00:00<00:59, 166.24 it/sec, obj=[18.35983854 24.29157804]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.21193996278027655
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.33061242e-08 3.32655148e-09]
INFO - 10:40:04: wrapped_ineq1 = [-11.33929244]
INFO - 10:40:04: wrapped_ineq2 = [-38.87119268]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.514919147231388 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.514919012190323 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2119399627802765 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [44.40816327 34.97959184]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.514919147231388 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.514919012190323 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2119399627802765 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 167/10000 [00:01<00:59, 164.00 it/sec, obj=[18.35983854 24.29157804]]
INFO - 10:40:04: 2%|▏ | 168/10000 [00:01<00:59, 164.51 it/sec, obj=[18.35983854 24.29157804]]
INFO - 10:40:04: 2%|▏ | 169/10000 [00:01<00:59, 165.08 it/sec, obj=[16.37447224 25.49140673]]
INFO - 10:40:04: 2%|▏ | 170/10000 [00:01<00:59, 165.69 it/sec, obj=[16.35853156 25.49177588]]
INFO - 10:40:04: 2%|▏ | 171/10000 [00:01<00:59, 166.30 it/sec, obj=[16.35661688 25.4914123 ]]
INFO - 10:40:04: 2%|▏ | 172/10000 [00:01<00:58, 166.83 it/sec, obj=[16.35645426 25.49151379]]
INFO - 10:40:04: 2%|▏ | 173/10000 [00:01<00:58, 167.41 it/sec, obj=[16.35645426 25.49151379]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.20626256632213952
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.28955868e-08 3.22389582e-09]
INFO - 10:40:04: wrapped_ineq1 = [-10.8243637]
INFO - 10:40:04: wrapped_ineq2 = [-39.89635594]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.429880060285158 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.429879917264565 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2062625663221395 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [41.63265306 35.91836735]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.429880060285158 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.429879917264565 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.2062625663221395 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 174/10000 [00:01<00:59, 164.61 it/sec, obj=[16.35645426 25.49151379]]
INFO - 10:40:04: 2%|▏ | 175/10000 [00:01<00:59, 165.18 it/sec, obj=[14.47666997 26.72698623]]
INFO - 10:40:04: 2%|▏ | 176/10000 [00:01<00:59, 165.76 it/sec, obj=[14.46080786 26.72744066]]
INFO - 10:40:04: 2%|▏ | 177/10000 [00:01<00:59, 166.31 it/sec, obj=[14.45889916 26.72708902]]
INFO - 10:40:04: 2%|▏ | 178/10000 [00:01<00:58, 166.85 it/sec, obj=[14.45874622 26.72719297]]
INFO - 10:40:04: 2%|▏ | 179/10000 [00:01<00:58, 167.35 it/sec, obj=[14.45874622 26.72719297]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.19980813861390445
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.24679609e-08 3.11699822e-09]
INFO - 10:40:04: wrapped_ineq1 = [-10.29202946]
INFO - 10:40:04: wrapped_ineq2 = [-40.94172201]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.344374754503482 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.344374603674137 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1998081386139045 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [38.85714286 36.85714286]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.344374754503482 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.344374603674137 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1998081386139045 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 180/10000 [00:01<00:59, 165.11 it/sec, obj=[14.45874622 26.72719297]]
INFO - 10:40:04: 2%|▏ | 181/10000 [00:01<00:59, 165.59 it/sec, obj=[14.45874622 26.72719297]]
INFO - 10:40:04: 2%|▏ | 182/10000 [00:01<00:59, 166.15 it/sec, obj=[12.68629432 27.99889971]]
INFO - 10:40:04: 2%|▏ | 183/10000 [00:01<00:58, 166.71 it/sec, obj=[12.67051626 27.99944249]]
INFO - 10:40:04: 2%|▏ | 184/10000 [00:01<00:58, 167.28 it/sec, obj=[12.66861499 27.99910359]]
INFO - 10:40:04: 2%|▏ | 185/10000 [00:01<00:58, 167.81 it/sec, obj=[12.66847181 27.99921 ]]
INFO - 10:40:04: 2%|▏ | 186/10000 [00:01<00:58, 168.36 it/sec, obj=[12.66847181 27.99921 ]]
INFO - 10:40:04: 2%|▏ | 187/10000 [00:01<00:58, 168.89 it/sec, obj=[12.66847181 27.99921 ]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 10
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.19256375777287138
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.20262058e-08 3.00653724e-09]
INFO - 10:40:04: wrapped_ineq1 = [-9.74200027]
INFO - 10:40:04: wrapped_ineq2 = [-42.00762687]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.258395476744471 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.2583953183017 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1925637577728714 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [36.08163265 37.79591837]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.258395476744471 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.2583953183017 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1925637577728714 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 188/10000 [00:01<00:58, 166.85 it/sec, obj=[12.66847181 27.99921 ]]
INFO - 10:40:04: 2%|▏ | 189/10000 [00:01<00:58, 167.27 it/sec, obj=[12.66847181 27.99921 ]]
INFO - 10:40:04: 2%|▏ | 190/10000 [00:01<00:58, 167.81 it/sec, obj=[11.00515155 29.30775805]]
INFO - 10:40:04: 2%|▏ | 191/10000 [00:01<00:58, 168.35 it/sec, obj=[10.98946331 29.30839237]]
INFO - 10:40:04: 2%|▏ | 192/10000 [00:01<00:58, 168.87 it/sec, obj=[10.98757099 29.30806707]]
INFO - 10:40:04: 2%|▏ | 193/10000 [00:01<00:57, 169.38 it/sec, obj=[10.98743766 29.30817594]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 7
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.18451613974807748
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.15729684e-08 2.89323232e-09]
INFO - 10:40:04: wrapped_ineq1 = [-9.17397857]
INFO - 10:40:04: wrapped_ineq2 = [-43.09441589]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.171934256566079 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.171934090728132 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1845161397480775 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [33.30612245 38.73469388]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.171934256566079 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.171934090728132 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1845161397480775 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 194/10000 [00:01<00:58, 167.08 it/sec, obj=[10.98743766 29.30817594]]
INFO - 10:40:04: 2%|▏ | 195/10000 [00:01<00:58, 167.61 it/sec, obj=[ 9.43509909 30.65418945]]
INFO - 10:40:04: 2%|▏ | 196/10000 [00:01<00:58, 168.13 it/sec, obj=[ 9.41950673 30.65491867]]
INFO - 10:40:04: 2%|▏ | 197/10000 [00:01<00:58, 168.63 it/sec, obj=[ 9.41762501 30.65460787]]
INFO - 10:40:04: 2%|▏ | 198/10000 [00:01<00:57, 169.12 it/sec, obj=[ 9.4175016 30.65471918]]
INFO - 10:40:04: 2%|▏ | 199/10000 [00:01<00:57, 169.63 it/sec, obj=[ 9.4175016 30.65471918]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.17565162394799738
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.11113430e-08 2.77787038e-09]
INFO - 10:40:04: wrapped_ineq1 = [-8.58765836]
INFO - 10:40:04: wrapped_ineq2 = [-44.20244425]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.084982897602836 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.084982724609076 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1756516239479974 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [30.53061224 39.67346939]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 1.084982897602836 | 5 | float |
INFO - 10:40:04: | y | 0 | 1.084982724609076 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1756516239479974 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 200/10000 [00:01<00:58, 167.29 it/sec, obj=[ 9.4175016 30.65471918]]
INFO - 10:40:04: 2%|▏ | 201/10000 [00:01<00:58, 167.78 it/sec, obj=[ 7.97804763 32.03884013]]
INFO - 10:40:04: 2%|▏ | 202/10000 [00:01<00:58, 168.29 it/sec, obj=[ 7.96255753 32.03966771]]
INFO - 10:40:04: 2%|▏ | 203/10000 [00:01<00:58, 168.81 it/sec, obj=[ 7.96068815 32.03937238]]
INFO - 10:40:04: 2%|▏ | 204/10000 [00:01<00:57, 169.29 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 10:40:04: 2%|▏ | 205/10000 [00:01<00:57, 169.78 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 10:40:04: 2%|▏ | 206/10000 [00:01<00:57, 170.27 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.16595615812650483
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.06442890e-08 2.66105626e-09]
INFO - 10:40:04: wrapped_ineq1 = [-7.98272487]
INFO - 10:40:04: wrapped_ineq2 = [-45.33207732]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.997532968497486 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.9975327886064471 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1659561581265048 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [27.75510204 40.6122449 ]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.997532968497486 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.9975327886064471 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1659561581265048 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 207/10000 [00:01<00:58, 168.41 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 10:40:04: 2%|▏ | 208/10000 [00:01<00:57, 168.83 it/sec, obj=[ 7.96057475 32.03948612]]
INFO - 10:40:04: 2%|▏ | 209/10000 [00:01<00:57, 169.32 it/sec, obj=[ 6.63596326 33.46237501]]
INFO - 10:40:04: 2%|▏ | 210/10000 [00:01<00:57, 169.82 it/sec, obj=[ 6.62058212 33.46330461]]
INFO - 10:40:04: 2%|▏ | 211/10000 [00:01<00:57, 170.32 it/sec, obj=[ 6.61872696 33.46302576]]
INFO - 10:40:04: 2%|▏ | 212/10000 [00:01<00:57, 170.79 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 10:40:04: 2%|▏ | 213/10000 [00:01<00:57, 171.26 it/sec, obj=[ 6.61862363 33.46314191]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.15541528247977868
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [1.01743218e-08 2.54355026e-09]
INFO - 10:40:04: wrapped_ineq1 = [-7.3588542]
INFO - 10:40:04: wrapped_ineq2 = [-46.48369104]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.9095757933591742 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.9095756068465055 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1554152824797787 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [24.97959184 41.55102041]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.9095757933591742 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.9095756068465055 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1554152824797787 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 214/10000 [00:01<00:57, 169.21 it/sec, obj=[ 6.61862363 33.46314191]]
INFO - 10:40:04: 2%|▏ | 215/10000 [00:01<00:57, 169.68 it/sec, obj=[ 5.41086973 34.92547856]]
INFO - 10:40:04: 2%|▏ | 216/10000 [00:01<00:57, 170.16 it/sec, obj=[ 5.39560462 34.92651398]]
INFO - 10:40:04: 2%|▏ | 217/10000 [00:01<00:57, 170.64 it/sec, obj=[ 5.39376567 34.92625268]]
INFO - 10:40:04: 2%|▏ | 218/10000 [00:01<00:57, 171.09 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 10:40:04: 2%|▏ | 219/10000 [00:01<00:57, 171.56 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 10:40:04: 2%|▏ | 220/10000 [00:01<00:56, 172.02 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 10
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.14401411290320923
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [9.70413438e-09 2.42602471e-09]
INFO - 10:40:04: wrapped_ineq1 = [-6.71571295]
INFO - 10:40:04: wrapped_ineq2 = [-47.6576724]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.8211024417176317 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.8211022488739839 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1440141129032092 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [22.20408163 42.48979592]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.8211024417176317 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.8211022488739839 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1440141129032092 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 221/10000 [00:01<00:57, 169.97 it/sec, obj=[ 5.39367249 34.92637122]]
INFO - 10:40:04: 2%|▏ | 222/10000 [00:01<00:57, 170.42 it/sec, obj=[ 4.30485089 36.42885552]]
INFO - 10:40:04: 2%|▏ | 223/10000 [00:01<00:57, 170.89 it/sec, obj=[ 4.28970922 36.43000074]]
INFO - 10:40:04: 2%|▏ | 224/10000 [00:01<00:57, 171.36 it/sec, obj=[ 4.28788864 36.42975814]]
INFO - 10:40:04: 2%|▏ | 225/10000 [00:01<00:56, 171.80 it/sec, obj=[ 4.28780567 36.42987905]]
INFO - 10:40:04: 2%|▏ | 226/10000 [00:01<00:56, 172.25 it/sec, obj=[ 4.28780567 36.42987905]]
INFO - 10:40:04: 2%|▏ | 227/10000 [00:01<00:56, 172.70 it/sec, obj=[ 4.28780567 36.42987905]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 10
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.131737323353427
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [9.23636279e-09 2.30908981e-09]
INFO - 10:40:04: wrapped_ineq1 = [-6.05295781]
INFO - 10:40:04: wrapped_ineq2 = [-48.85441989]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.732103717940411 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.7321035190692963 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.131737323353427 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [19.42857143 43.42857143]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.732103717940411 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.7321035190692963 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.131737323353427 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 228/10000 [00:01<00:57, 170.73 it/sec, obj=[ 4.28780567 36.42987905]]
INFO - 10:40:04: 2%|▏ | 229/10000 [00:01<00:57, 171.18 it/sec, obj=[ 3.32005317 37.97323185]]
INFO - 10:40:04: 2%|▏ | 230/10000 [00:01<00:56, 171.64 it/sec, obj=[ 3.30504278 37.97449102]]
INFO - 10:40:04: 2%|▏ | 231/10000 [00:01<00:56, 172.08 it/sec, obj=[ 3.30324284 37.97426832]]
INFO - 10:40:04: 2%|▏ | 232/10000 [00:01<00:56, 172.52 it/sec, obj=[ 3.30317013 37.97439158]]
INFO - 10:40:04: 2%|▏ | 233/10000 [00:01<00:56, 172.94 it/sec, obj=[ 3.30317013 37.97439158]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.11856912725629419
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [8.77312090e-09 2.19328733e-09]
INFO - 10:40:04: wrapped_ineq1 = [-5.37023516]
INFO - 10:40:04: wrapped_ineq2 = [-50.07434394]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.6425701500776614 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.6425699454933228 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1185691272562942 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [16.65306122 44.36734694]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.6425701500776614 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.6425699454933228 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1185691272562942 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 234/10000 [00:01<00:57, 171.04 it/sec, obj=[ 3.30317013 37.97439158]]
INFO - 10:40:04: 2%|▏ | 235/10000 [00:01<00:56, 171.47 it/sec, obj=[ 2.45868826 39.55935553]]
INFO - 10:40:04: 2%|▏ | 236/10000 [00:01<00:56, 171.91 it/sec, obj=[ 2.44381739 39.56073302]]
INFO - 10:40:04: 2%|▏ | 237/10000 [00:01<00:56, 172.35 it/sec, obj=[ 2.44204056 39.56053151]]
INFO - 10:40:04: 2%|▏ | 238/10000 [00:01<00:56, 172.77 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 10:40:04: 2%|▏ | 239/10000 [00:01<00:56, 173.19 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 10:40:04: 2%|▏ | 240/10000 [00:01<00:56, 173.61 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.1044932578970259
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [8.31648883e-09 2.07914130e-09]
INFO - 10:40:04: wrapped_ineq1 = [-4.66718063]
INFO - 10:40:04: wrapped_ineq2 = [-51.3178675]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.5524919780961888 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.5524917681215171 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1044932578970259 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [13.87755102 45.30612245]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.5524919780961888 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.5524917681215171 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.1044932578970259 | inf | float |
INFO - 10:40:04: +------+-------------+--------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 241/10000 [00:01<00:56, 171.69 it/sec, obj=[ 2.44197816 39.56065708]]
INFO - 10:40:04: 2%|▏ | 242/10000 [00:01<00:56, 172.10 it/sec, obj=[ 1.72303596 41.18799761]]
INFO - 10:40:04: 2%|▏ | 243/10000 [00:01<00:56, 172.53 it/sec, obj=[ 1.70831329 41.18949797]]
INFO - 10:40:04: 2%|▏ | 244/10000 [00:01<00:56, 172.96 it/sec, obj=[ 1.70656219 41.18931901]]
INFO - 10:40:04: 2%|▏ | 245/10000 [00:01<00:56, 173.37 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 10:40:04: 2%|▏ | 246/10000 [00:01<00:56, 173.77 it/sec, obj=[ 1.70651014 41.18944686]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.08949294772348224
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [7.86845433e-09 1.96711625e-09]
INFO - 10:40:04: wrapped_ineq1 = [-3.94341862]
INFO - 10:40:04: wrapped_ineq2 = [-52.58542653]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.4618591414613932 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.4618589264259151 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.08949294772348224 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [11.10204082 46.24489796]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.4618591414613932 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.4618589264259151 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.08949294772348224 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: 2%|▏ | 247/10000 [00:01<00:56, 172.10 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 10:40:04: 2%|▏ | 248/10000 [00:01<00:56, 172.45 it/sec, obj=[ 1.70651014 41.18944686]]
INFO - 10:40:04: 2%|▏ | 249/10000 [00:01<00:56, 172.87 it/sec, obj=[ 1.11544708 42.85995314]]
INFO - 10:40:04: 2%|▎ | 250/10000 [00:01<00:56, 173.29 it/sec, obj=[ 1.10088177 42.86158114]]
INFO - 10:40:04: 3%|▎ | 251/10000 [00:01<00:56, 173.71 it/sec, obj=[ 1.09915921 42.86142618]]
INFO - 10:40:04: 3%|▎ | 252/10000 [00:01<00:55, 174.10 it/sec, obj=[ 1.09911754 42.86155626]]
INFO - 10:40:04: 3%|▎ | 253/10000 [00:01<00:55, 174.50 it/sec, obj=[ 1.09911754 42.86155626]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 9
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.07355090648809494
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [7.43084527e-09 1.85768645e-09]
INFO - 10:40:04: wrapped_ineq1 = [-3.19856184]
INFO - 10:40:04: wrapped_ineq2 = [-53.87747056]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.3706612660223688 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.3706610462603387 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.07355090648809494 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [ 8.32653061 47.18367347]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.3706612660223688 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.3706610462603387 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.07355090648809494 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: 3%|▎ | 254/10000 [00:01<00:56, 172.48 it/sec, obj=[ 1.09911754 42.86155626]]
INFO - 10:40:04: 3%|▎ | 255/10000 [00:01<00:56, 172.88 it/sec, obj=[ 0.63834664 44.57604226]]
INFO - 10:40:04: 3%|▎ | 256/10000 [00:01<00:56, 173.29 it/sec, obj=[ 0.62394835 44.57780289]]
INFO - 10:40:04: 3%|▎ | 257/10000 [00:01<00:56, 173.70 it/sec, obj=[ 0.62225734 44.57767347]]
INFO - 10:40:04: 3%|▎ | 258/10000 [00:01<00:55, 174.08 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 10:40:04: 3%|▎ | 259/10000 [00:01<00:55, 174.48 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 10:40:04: 3%|▎ | 260/10000 [00:01<00:55, 174.87 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 10
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.05664929814776543
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [7.00444502e-09 1.75109349e-09]
INFO - 10:40:04: wrapped_ineq1 = [-2.43221075]
INFO - 10:40:04: wrapped_ineq2 = [-55.19446334]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.2788876501517835 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.2788874260003971 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.05664929814776543 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [ 5.55102041 48.12244898]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.2788876501517835 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.2788874260003971 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.05664929814776543 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: 3%|▎ | 261/10000 [00:01<00:56, 173.07 it/sec, obj=[ 0.62222607 44.57780576]]
INFO - 10:40:04: 3%|▎ | 262/10000 [00:01<00:56, 173.46 it/sec, obj=[ 0.29423714 46.33711132]]
INFO - 10:40:04: 3%|▎ | 263/10000 [00:01<00:56, 173.86 it/sec, obj=[ 0.28001607 46.33900984]]
INFO - 10:40:04: 3%|▎ | 264/10000 [00:01<00:55, 174.26 it/sec, obj=[ 0.27835982 46.3389076 ]]
INFO - 10:40:04: 3%|▎ | 265/10000 [00:01<00:55, 174.64 it/sec, obj=[ 0.27833898 46.33904203]]
INFO - 10:40:04: 3%|▎ | 266/10000 [00:01<00:55, 175.02 it/sec, obj=[ 0.27833898 46.33904203]]
ERROR - 10:40:04: 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.9/site-packages/gemseo/algos/opt/lib_nlopt.py", line 471, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/develop/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 8
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.03876971643439605
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [6.59082389e-09 1.64768976e-09]
INFO - 10:40:04: wrapped_ineq1 = [-1.64395306]
INFO - 10:40:04: wrapped_ineq2 = [-56.53688344]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.1865272500881318 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.186527021885882 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.03876971643439605 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: Solving mNBI sub-problem for phi_beta = [ 2.7755102 49.06122449]
INFO - 10:40:04: Optimization problem:
INFO - 10:40:04: minimize -t_extraction
INFO - 10:40:04: with respect to t, x, y
INFO - 10:40:04: subject to constraints:
INFO - 10:40:04: wrapped_ineq1 <= 0.0
INFO - 10:40:04: wrapped_ineq2 <= 0.0
INFO - 10:40:04: beta_sub_optim_constraint <= 0.0
INFO - 10:40:04: over the design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.1865272500881318 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.186527021885882 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.03876971643439605 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: 3%|▎ | 267/10000 [00:01<00:56, 173.45 it/sec, obj=[ 0.27833898 46.33904203]]
INFO - 10:40:04: 3%|▎ | 268/10000 [00:01<00:56, 173.77 it/sec, obj=[ 0.27833898 46.33904203]]
INFO - 10:40:04: 3%|▎ | 269/10000 [00:01<00:55, 174.15 it/sec, obj=[ 0.08570212 48.14403408]]
INFO - 10:40:04: 3%|▎ | 270/10000 [00:01<00:55, 174.53 it/sec, obj=[ 0.07166903 48.14607599]]
INFO - 10:40:04: 3%|▎ | 271/10000 [00:01<00:55, 174.92 it/sec, obj=[ 0.07005101 48.14600266]]
INFO - 10:40:04: 3%|▎ | 272/10000 [00:01<00:55, 175.28 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 10:40:04: 3%|▎ | 273/10000 [00:01<00:55, 175.65 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 10:40:04: 3%|▎ | 274/10000 [00:01<00:55, 176.02 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 10
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: -0.019893159001294167
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: beta_sub_optim_constraint = [6.19088514e-09 1.54772239e-09]
INFO - 10:40:04: wrapped_ineq1 = [-0.83336312]
INFO - 10:40:04: wrapped_ineq2 = [-57.9052249]
INFO - 10:40:04: Design space:
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | Name | Lower bound | Value | Upper bound | Type |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: | x | 0 | 0.09356866442349428 | 5 | float |
INFO - 10:40:04: | y | 0 | 0.09356843250868149 | 3 | float |
INFO - 10:40:04: | t | -inf | 0.01989315900129417 | inf | float |
INFO - 10:40:04: +------+-------------+---------------------+-------------+-------+
INFO - 10:40:04: 3%|▎ | 275/10000 [00:01<00:57, 170.44 it/sec, obj=[ 0.07004059 48.14613918]]
INFO - 10:40:04: Optimization result:
INFO - 10:40:04: Optimizer info:
INFO - 10:40:04: Status: None
INFO - 10:40:04: Message: None
INFO - 10:40:04: Number of calls to the objective function by the optimizer: 402
INFO - 10:40:04: Solution:
INFO - 10:40:04: The solution is feasible.
INFO - 10:40:04: Objective: 30.287800698724556
INFO - 10:40:04: Standardized constraints:
INFO - 10:40:04: ineq1 = [-10.82436329]
INFO - 10:40:04: ineq2 = [-39.89635666]
INFO - 10:40:04: Pareto efficient solutions:
INFO - 10:40:04: Pareto optimal points : 203 / 275
INFO - 10:40:04: Utopia point : [2.61376505e-22 4.00000000e+00]
INFO - 10:40:04: Compromise solution (closest to utopia) : [[14.45889916 26.72708902]]
INFO - 10:40:04: Distance from utopia : 26.936598526002314
INFO - 10:40:04: Objective values:
INFO - 10:40:04: +------+-------------+-----+---------+
INFO - 10:40:04: | name | 1 | 2 | 3 |
INFO - 10:40:04: +------+-------------+-----+---------+
INFO - 10:40:04: | 1 | 2.61377e-22 | 136 | 14.4589 |
INFO - 10:40:04: | 2 | 50 | 4 | 26.7271 |
INFO - 10:40:04: +------+-------------+-----+---------+
INFO - 10:40:04: Constraint values:
INFO - 10:40:04: +------+-------------+------+----------+
INFO - 10:40:04: | name | 1 | 2 | 3 |
INFO - 10:40:04: +------+-------------+------+----------+
INFO - 10:40:04: | 1 | 3.61311e-12 | -22 | -10.2918 |
INFO - 10:40:04: | 2 | -59.3 | -7.3 | -40.9422 |
INFO - 10:40:04: +------+-------------+------+----------+
INFO - 10:40:04: Design space:
INFO - 10:40:04: +-------+-------------+-------------+---+---------+-------------+-------+
INFO - 10:40:04: | name | lower_bound | 1 | 2 | 3 | upper_bound | type |
INFO - 10:40:04: +-------+-------------+-------------+---+---------+-------------+-------+
INFO - 10:40:04: | x (1) | 0 | 4.45921e-13 | 5 | 1.34434 | 5 | float |
INFO - 10:40:04: | y (1) | 0 | 8.07126e-12 | 3 | 1.34442 | 3 | float |
INFO - 10:40:04: +-------+-------------+-------------+---+---------+-------------+-------+
Display the Pareto front¶
GEMSEO detects the Pareto optimal points and the dominated ones.
execute_post(problem, "ParetoFront", save=False, show=True)
<gemseo.post.pareto_front.ParetoFront object at 0x7f13201ff700>
Total running time of the script: (0 minutes 1.903 seconds)