Note
Go to the end to download the full example code
Multistart optimization¶
Runs simple optimization problem with multiple starting points
Nests a MDOScenario
in a DOEScenario
using a MDOScenarioAdapter
.
from __future__ import annotations
from gemseo import configure_logger
from gemseo import create_design_space
from gemseo import create_discipline
from gemseo import create_scenario
from gemseo.disciplines.scenario_adapters.mdo_scenario_adapter import MDOScenarioAdapter
configure_logger()
<RootLogger root (INFO)>
Create the disciplines¶
objective = create_discipline("AnalyticDiscipline", expressions={"obj": "x**3-x+1"})
constraint = create_discipline(
"AnalyticDiscipline", expressions={"cstr": "x**2+obj**2-1.5"}
)
Create the design space¶
design_space = create_design_space()
design_space.add_variable("x", l_b=-1.5, u_b=1.5, value=1.5)
Create the MDO scenario¶
scenario = create_scenario(
[objective, constraint],
formulation="DisciplinaryOpt",
objective_name="obj",
design_space=design_space,
)
scenario.default_inputs = {"algo": "SLSQP", "max_iter": 10}
scenario.add_constraint("cstr", "ineq")
Create the scenario adapter¶
dv_names = scenario.formulation.opt_problem.design_space.variable_names
adapter = MDOScenarioAdapter(
scenario, dv_names, ["obj", "cstr"], set_x0_before_opt=True
)
Create the DOE scenario¶
scenario_doe = create_scenario(
adapter,
formulation="DisciplinaryOpt",
objective_name="obj",
design_space=design_space,
scenario_type="DOE",
)
scenario_doe.add_constraint("cstr", "ineq")
run_inputs = {"n_samples": 10, "algo": "fullfact"}
scenario_doe.execute(run_inputs)
INFO - 08:37:35:
INFO - 08:37:35: *** Start DOEScenario execution ***
INFO - 08:37:35: DOEScenario
INFO - 08:37:35: Disciplines: MDOScenario_adapter
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 1.5 | 1.5 | float |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm fullfact:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -1.5 | 1.5 | float |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 137.73 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 172.99 it/sec, obj=-.267]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 224.43 it/sec, obj=-.809]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 264.30 it/sec, obj=-.868]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 296.43 it/sec, obj=-.872]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 304.00 it/sec, obj=-.265]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 309.99 it/sec, obj=0.136]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 314.89 it/sec, obj=0.579]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 329.54 it/sec, obj=0.289]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 331.83 it/sec, obj=1.25]
INFO - 08:37:35:
WARNING - 08:37:35: Optimization found no feasible point ! The least infeasible point is selected.
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:35: Solution:
WARNING - 08:37:35: The solution is not feasible.
INFO - 08:37:35: Objective: 0.2888129873625884
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = 0.1513713660745195
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -1.252181466244097 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.041569) ***
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 20.82 it/sec, obj=0.289]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -1.166666666666667 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 919.60 it/sec, obj=0.579]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 434.35 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 399.06 it/sec, obj=-.267]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 431.50 it/sec, obj=-.809]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 455.31 it/sec, obj=-.868]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 470.64 it/sec, obj=-.874]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 449.88 it/sec, obj=-.266]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 435.23 it/sec, obj=0.135]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 424.28 it/sec, obj=0.577]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 436.26 it/sec, obj=0.288]
INFO - 08:37:35:
WARNING - 08:37:35: Optimization found no feasible point ! The least infeasible point is selected.
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:35: Solution:
WARNING - 08:37:35: The solution is not feasible.
INFO - 08:37:35: Objective: 0.28811314244670916
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = 0.15144075010172386
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -1.252370379421043 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.034791) ***
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 23.72 it/sec, obj=0.288]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -0.8333333333333334 | 1.5 | float |
INFO - 08:37:35: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 929.79 it/sec, obj=1.25]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 427.44 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 464.25 it/sec, obj=1.01]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 428.53 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 451.54 it/sec, obj=0.819]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 424.29 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 440.84 it/sec, obj=0.693]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 427.24 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 439.52 it/sec, obj=0.584]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 429.54 it/sec, obj=-.875]
INFO - 08:37:35:
WARNING - 08:37:35: Optimization found no feasible point ! The least infeasible point is selected.
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:35: Solution:
WARNING - 08:37:35: The solution is not feasible.
INFO - 08:37:35: Objective: 0.5837798748832244
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = 0.19806414473664136
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -1.165017254128868 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.034288) ***
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 24.83 it/sec, obj=0.584]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -0.5 | 1.5 | float |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 925.89 it/sec, obj=1.38]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 413.88 it/sec, obj=2.88]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 468.74 it/sec, obj=1.23]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 431.15 it/sec, obj=2.87]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 454.13 it/sec, obj=0.848]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 426.75 it/sec, obj=0.658]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 402.51 it/sec, obj=0.643]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 387.80 it/sec, obj=0.616]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 377.57 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 369.71 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 13
INFO - 08:37:35: Solution:
INFO - 08:37:35: The solution is feasible.
INFO - 08:37:35: Objective: 0.6150998219543254
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = -0.788285881879476
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.5773788419679762 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.038752) ***
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 24.84 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -0.1666666666666667 | 1.5 | float |
INFO - 08:37:35: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 960.67 it/sec, obj=1.16]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 434.82 it/sec, obj=2.87]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 487.03 it/sec, obj=0.785]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 423.32 it/sec, obj=2.88]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 430.79 it/sec, obj=0.644]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 405.20 it/sec, obj=0.624]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 387.73 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 376.22 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 365.71 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 358.77 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:35: Solution:
INFO - 08:37:35: The solution is feasible.
INFO - 08:37:35: Objective: 0.6150998205402495
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = -0.7883188793606977
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.5773502675245377 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.039658) ***
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 24.74 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.1666666666666667 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 967.32 it/sec, obj=0.838]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 436.57 it/sec, obj=2.88]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 486.39 it/sec, obj=0.656]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 425.75 it/sec, obj=0.639]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 394.97 it/sec, obj=0.616]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 377.66 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 366.80 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 359.48 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 353.65 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 10
INFO - 08:37:35: Solution:
INFO - 08:37:35: The solution is feasible.
INFO - 08:37:35: Objective: 0.6150998205402508
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = -0.7883189092934597
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.5773502416020033 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.035369) ***
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 25.03 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.5 | 1.5 | float |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 960.45 it/sec, obj=0.625]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 433.41 it/sec, obj=2.87]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 480.47 it/sec, obj=0.616]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 422.64 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 393.23 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 375.13 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 364.60 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 8
INFO - 08:37:35: Solution:
INFO - 08:37:35: The solution is feasible.
INFO - 08:37:35: Objective: 0.6150998205402526
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = -0.7883188285028959
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.5773503115686811 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.029191) ***
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 25.82 it/sec, obj=0.615]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.8333333333333335 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 946.15 it/sec, obj=0.745]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 435.18 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 403.57 it/sec, obj=-.267]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 418.33 it/sec, obj=-.809]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 442.19 it/sec, obj=-.868]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 459.98 it/sec, obj=-.874]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 440.64 it/sec, obj=-.266]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 427.63 it/sec, obj=0.135]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 418.42 it/sec, obj=0.577]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 429.95 it/sec, obj=0.288]
INFO - 08:37:35:
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:35: Solution:
INFO - 08:37:35: The solution is feasible.
INFO - 08:37:35: Objective: 0.7453703703703706
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = -0.2499785665294918
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 0.8333333333333335 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.034541) ***
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 25.98 it/sec, obj=0.745]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+-------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+-------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 1.166666666666667 | 1.5 | float |
INFO - 08:37:35: +------+-------------+-------------------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 931.65 it/sec, obj=1.42]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 438.14 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:35: ... 30%|███ | 3/10 [00:00<00:00, 409.28 it/sec, obj=-.267]
INFO - 08:37:35:
INFO - 08:37:35: ... 40%|████ | 4/10 [00:00<00:00, 441.23 it/sec, obj=-.809]
INFO - 08:37:35:
INFO - 08:37:35: ... 50%|█████ | 5/10 [00:00<00:00, 464.69 it/sec, obj=-.868]
INFO - 08:37:35:
INFO - 08:37:35: ... 60%|██████ | 6/10 [00:00<00:00, 481.22 it/sec, obj=-.874]
INFO - 08:37:35:
INFO - 08:37:35: ... 70%|███████ | 7/10 [00:00<00:00, 457.79 it/sec, obj=-.266]
INFO - 08:37:35:
INFO - 08:37:35: ... 80%|████████ | 8/10 [00:00<00:00, 434.79 it/sec, obj=0.135]
INFO - 08:37:35:
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 424.16 it/sec, obj=0.577]
INFO - 08:37:35:
INFO - 08:37:35: ... 100%|██████████| 10/10 [00:00<00:00, 436.35 it/sec, obj=0.288]
INFO - 08:37:35:
WARNING - 08:37:35: Optimization found no feasible point ! The least infeasible point is selected.
INFO - 08:37:35:
INFO - 08:37:35: Optimization result:
INFO - 08:37:35: Optimizer info:
INFO - 08:37:35: Status: None
INFO - 08:37:35: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:35: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:35: Solution:
WARNING - 08:37:35: The solution is not feasible.
INFO - 08:37:35: Objective: 0.28811314244698893
INFO - 08:37:35: Standardized constraints:
INFO - 08:37:35: cstr = 0.1514407501016959
INFO - 08:37:35: Design space:
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | -1.252370379420967 | 1.5 | float |
INFO - 08:37:35: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:35: *** End MDOScenario execution (time: 0:00:00.035059) ***
INFO - 08:37:35: ... 90%|█████████ | 9/10 [00:00<00:00, 26.14 it/sec, obj=0.288]
INFO - 08:37:35:
INFO - 08:37:35: *** Start MDOScenario execution ***
INFO - 08:37:35: MDOScenario
INFO - 08:37:35: Disciplines: AnalyticDiscipline AnalyticDiscipline
INFO - 08:37:35: MDO formulation: DisciplinaryOpt
INFO - 08:37:35: Optimization problem:
INFO - 08:37:35: minimize obj(x)
INFO - 08:37:35: with respect to x
INFO - 08:37:35: subject to constraints:
INFO - 08:37:35: cstr(x) <= 0.0
INFO - 08:37:35: over the design space:
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: | x | -1.5 | 1.5 | 1.5 | float |
INFO - 08:37:35: +------+-------------+-------+-------------+-------+
INFO - 08:37:35: Solving optimization problem with algorithm SLSQP:
INFO - 08:37:35: ... 0%| | 0/10 [00:00<?, ?it]
INFO - 08:37:35:
INFO - 08:37:35: ... 10%|█ | 1/10 [00:00<00:00, 965.10 it/sec, obj=2.88]
INFO - 08:37:35:
INFO - 08:37:35: ... 20%|██ | 2/10 [00:00<00:00, 437.16 it/sec, obj=-.875]
INFO - 08:37:35:
INFO - 08:37:36: ... 30%|███ | 3/10 [00:00<00:00, 412.69 it/sec, obj=-.267]
INFO - 08:37:36:
INFO - 08:37:36: ... 40%|████ | 4/10 [00:00<00:00, 442.22 it/sec, obj=-.809]
INFO - 08:37:36:
INFO - 08:37:36: ... 50%|█████ | 5/10 [00:00<00:00, 464.39 it/sec, obj=-.868]
INFO - 08:37:36:
INFO - 08:37:36: ... 60%|██████ | 6/10 [00:00<00:00, 479.33 it/sec, obj=-.874]
INFO - 08:37:36:
INFO - 08:37:36: ... 70%|███████ | 7/10 [00:00<00:00, 456.91 it/sec, obj=-.266]
INFO - 08:37:36:
INFO - 08:37:36: ... 80%|████████ | 8/10 [00:00<00:00, 440.53 it/sec, obj=0.135]
INFO - 08:37:36:
INFO - 08:37:36: ... 90%|█████████ | 9/10 [00:00<00:00, 428.83 it/sec, obj=0.577]
INFO - 08:37:36:
INFO - 08:37:36: ... 100%|██████████| 10/10 [00:00<00:00, 440.16 it/sec, obj=0.288]
INFO - 08:37:36:
WARNING - 08:37:36: Optimization found no feasible point ! The least infeasible point is selected.
INFO - 08:37:36:
INFO - 08:37:36: Optimization result:
INFO - 08:37:36: Optimizer info:
INFO - 08:37:36: Status: None
INFO - 08:37:36: Message: Maximum number of iterations reached. GEMSEO Stopped the driver
INFO - 08:37:36: Number of calls to the objective function by the optimizer: 12
INFO - 08:37:36: Solution:
WARNING - 08:37:36: The solution is not feasible.
INFO - 08:37:36: Objective: 0.28811314244587316
INFO - 08:37:36: Standardized constraints:
INFO - 08:37:36: cstr = 0.15144075010180735
INFO - 08:37:36: Design space:
INFO - 08:37:36: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:36: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:36: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:36: | x | -1.5 | -1.252370379421268 | 1.5 | float |
INFO - 08:37:36: +------+-------------+--------------------+-------------+-------+
INFO - 08:37:36: *** End MDOScenario execution (time: 0:00:00.035634) ***
INFO - 08:37:36: ... 100%|██████████| 10/10 [00:00<00:00, 26.22 it/sec, obj=0.288]
INFO - 08:37:36: Optimization result:
INFO - 08:37:36: Optimizer info:
INFO - 08:37:36: Status: None
INFO - 08:37:36: Message: None
INFO - 08:37:36: Number of calls to the objective function by the optimizer: 10
INFO - 08:37:36: Solution:
INFO - 08:37:36: The solution is feasible.
INFO - 08:37:36: Objective: 0.6150998205402495
INFO - 08:37:36: Standardized constraints:
INFO - 08:37:36: cstr = -0.7883188793606977
INFO - 08:37:36: Design space:
INFO - 08:37:36: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:36: | name | lower_bound | value | upper_bound | type |
INFO - 08:37:36: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:36: | x | -1.5 | -0.1666666666666667 | 1.5 | float |
INFO - 08:37:36: +------+-------------+---------------------+-------------+-------+
INFO - 08:37:36: *** End DOEScenario execution (time: 0:00:00.392954) ***
{'eval_jac': False, 'n_samples': 10, 'algo': 'fullfact'}
Plot the optimum objective for different x0¶
scenario_doe.post_process("BasicHistory", variable_names=["obj"], save=False, show=True)
/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.0.1/lib/python3.9/site-packages/genson/schema/strategies/base.py:42: UserWarning: Schema incompatible. Keyword 'description' has conflicting values ('The width and height of the figure in inches, e.g. ``(w, h)``.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.' vs. 'The width and height of the figure in inches, e.g. `(w, h)`.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.'). Using 'The width and height of the figure in inches, e.g. ``(w, h)``.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.'
warn(('Schema incompatible. Keyword {0!r} has conflicting '
<gemseo.post.basic_history.BasicHistory object at 0x7f873adaf4f0>
Total running time of the script: ( 0 minutes 0.603 seconds)