Multistart optimization

Runs simple optimization problem with multiple starting points Nests an MDOScenario in a DOEScenario using an 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],
    "DisciplinaryOpt",
    "obj",
    design_space,
)
scenario.default_inputs = {"algo": "SLSQP", "max_iter": 10}
scenario.add_constraint("cstr", constraint_type="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,
    "DisciplinaryOpt",
    "obj",
    design_space,
    scenario_type="DOE",
)
scenario_doe.add_constraint("cstr", constraint_type="ineq")
run_inputs = {"n_samples": 10, "algo": "fullfact"}
scenario_doe.execute(run_inputs)
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start DOEScenario execution ***
    INFO - 08:56:30: DOEScenario
    INFO - 08:56:30:    Disciplines: MDOScenario_adapter
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound | Value | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    |  1.5  |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm fullfact:
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound | Value | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    |  -1.5 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 171.52 it/sec, obj=-0.875]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 211.79 it/sec, obj=-0.267]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 280.46 it/sec, obj=-0.809]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 330.95 it/sec, obj=-0.868]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 369.93 it/sec, obj=-0.872]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 381.79 it/sec, obj=-0.265]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 388.97 it/sec, obj=0.136]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 396.48 it/sec, obj=0.579]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 416.66 it/sec, obj=0.289]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 419.82 it/sec, obj=1.25]
 WARNING - 08:56:30: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
 WARNING - 08:56:30:       The solution is not feasible.
    INFO - 08:56:30:       Objective: 0.2888129873625884
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = 0.1513713660745195
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | -1.252181466244097 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.036041) ***
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 24.44 it/sec, obj=0.289]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    | -1.166666666666667 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1112.55 it/sec, obj=0.579]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 529.42 it/sec, obj=-0.875]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 512.35 it/sec, obj=-0.267]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 555.87 it/sec, obj=-0.809]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 580.37 it/sec, obj=-0.868]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 600.46 it/sec, obj=-0.874]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 571.51 it/sec, obj=-0.266]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 553.00 it/sec, obj=0.135]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 537.20 it/sec, obj=0.577]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 550.41 it/sec, obj=0.288]
 WARNING - 08:56:30: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
 WARNING - 08:56:30:       The solution is not feasible.
    INFO - 08:56:30:       Objective: 0.28811314244670916
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = 0.15144075010172386
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | -1.252370379421043 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.029255) ***
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 28.04 it/sec, obj=0.288]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound |        Value        | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    | -0.8333333333333334 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1122.07 it/sec, obj=1.25]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 531.53 it/sec, obj=-0.875]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 601.36 it/sec, obj=1.01]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 553.32 it/sec, obj=-0.875]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 582.61 it/sec, obj=0.819]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 553.19 it/sec, obj=-0.875]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 563.80 it/sec, obj=0.693]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 532.16 it/sec, obj=-0.875]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 548.21 it/sec, obj=0.584]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 534.97 it/sec, obj=-0.875]
 WARNING - 08:56:30: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
 WARNING - 08:56:30:       The solution is not feasible.
    INFO - 08:56:30:       Objective: 0.5837798748832244
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = 0.19806414473664136
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | -1.165017254128868 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.028935) ***
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 29.39 it/sec, obj=0.584]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound | Value | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    |  -0.5 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1116.40 it/sec, obj=1.38]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 534.54 it/sec, obj=2.88]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 604.45 it/sec, obj=1.23]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 553.69 it/sec, obj=2.87]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 584.21 it/sec, obj=0.848]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 548.49 it/sec, obj=0.658]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 517.38 it/sec, obj=0.643]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 496.80 it/sec, obj=0.616]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 484.51 it/sec, obj=0.615]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 474.25 it/sec, obj=0.615]
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 13
    INFO - 08:56:30:    Solution:
    INFO - 08:56:30:       The solution is feasible.
    INFO - 08:56:30:       Objective: 0.6150998219543254
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = -0.788285881879476
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | 0.5773788419679762 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.031886) ***
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 29.64 it/sec, obj=0.615]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound |        Value        | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    | -0.1666666666666667 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1140.07 it/sec, obj=1.16]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 534.58 it/sec, obj=2.87]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 603.15 it/sec, obj=0.785]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 531.13 it/sec, obj=2.88]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 540.80 it/sec, obj=0.644]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 508.91 it/sec, obj=0.624]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 486.53 it/sec, obj=0.615]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 473.99 it/sec, obj=0.615]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 463.65 it/sec, obj=0.615]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 455.87 it/sec, obj=0.615]
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
    INFO - 08:56:30:       The solution is feasible.
    INFO - 08:56:30:       Objective: 0.6150998205402495
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = -0.7883188793606977
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | 0.5773502675245377 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.032869) ***
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 29.61 it/sec, obj=0.615]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    | 0.1666666666666667 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1106.38 it/sec, obj=0.838]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 540.82 it/sec, obj=2.88]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 601.33 it/sec, obj=0.656]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 534.77 it/sec, obj=0.639]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 495.71 it/sec, obj=0.616]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 475.14 it/sec, obj=0.615]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 459.20 it/sec, obj=0.615]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 449.14 it/sec, obj=0.615]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 441.35 it/sec, obj=0.615]
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 10
    INFO - 08:56:30:    Solution:
    INFO - 08:56:30:       The solution is feasible.
    INFO - 08:56:30:       Objective: 0.6150998205402508
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = -0.7883189092934597
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | 0.5773502416020033 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.029859) ***
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 29.93 it/sec, obj=0.615]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound | Value | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    |  0.5  |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1039.48 it/sec, obj=0.625]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 533.46 it/sec, obj=2.87]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 594.71 it/sec, obj=0.616]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 531.06 it/sec, obj=0.615]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 494.24 it/sec, obj=0.615]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 474.79 it/sec, obj=0.615]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 463.26 it/sec, obj=0.615]
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 8
    INFO - 08:56:30:    Solution:
    INFO - 08:56:30:       The solution is feasible.
    INFO - 08:56:30:       Objective: 0.6150998205402526
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = -0.7883188285028959
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | 0.5773503115686811 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.024517) ***
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 30.87 it/sec, obj=0.615]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    | 0.8333333333333335 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1111.66 it/sec, obj=0.745]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 535.71 it/sec, obj=-0.875]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 511.40 it/sec, obj=-0.267]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 528.68 it/sec, obj=-0.809]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 560.15 it/sec, obj=-0.868]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 580.29 it/sec, obj=-0.874]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 558.50 it/sec, obj=-0.266]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 540.11 it/sec, obj=0.135]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 529.05 it/sec, obj=0.577]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 542.96 it/sec, obj=0.288]
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
    INFO - 08:56:30:       The solution is feasible.
    INFO - 08:56:30:       Objective: 0.7453703703703706
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = -0.2499785665294918
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | 0.8333333333333335 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.029043) ***
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 31.05 it/sec, obj=0.745]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+-------------------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound |       Value       | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+-------------------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    | 1.166666666666667 |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+-------------------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1111.37 it/sec, obj=1.42]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 534.03 it/sec, obj=-0.875]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 517.20 it/sec, obj=-0.267]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 555.52 it/sec, obj=-0.809]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 584.15 it/sec, obj=-0.868]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 605.17 it/sec, obj=-0.874]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 575.75 it/sec, obj=-0.266]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 514.69 it/sec, obj=0.135]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 506.97 it/sec, obj=0.577]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 522.78 it/sec, obj=0.288]
 WARNING - 08:56:30: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
 WARNING - 08:56:30:       The solution is not feasible.
    INFO - 08:56:30:       Objective: 0.28811314244698893
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = 0.1514407501016959
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | -1.252370379420967 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.030023) ***
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 31.18 it/sec, obj=0.288]
    INFO - 08:56:30:
    INFO - 08:56:30: *** Start MDOScenario execution ***
    INFO - 08:56:30: MDOScenario
    INFO - 08:56:30:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 08:56:30:    MDO formulation: DisciplinaryOpt
    INFO - 08:56:30: Optimization problem:
    INFO - 08:56:30:    minimize obj(x)
    INFO - 08:56:30:    with respect to x
    INFO - 08:56:30:    subject to constraints:
    INFO - 08:56:30:       cstr(x) <= 0.0
    INFO - 08:56:30:    over the design space:
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | Name | Lower bound | Value | Upper bound | Type  |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30:       | x    |     -1.5    |  1.5  |     1.5     | float |
    INFO - 08:56:30:       +------+-------------+-------+-------------+-------+
    INFO - 08:56:30: Solving optimization problem with algorithm SLSQP:
    INFO - 08:56:30:     10%|█         | 1/10 [00:00<00:00, 1043.88 it/sec, obj=2.88]
    INFO - 08:56:30:     20%|██        | 2/10 [00:00<00:00, 533.90 it/sec, obj=-0.875]
    INFO - 08:56:30:     30%|███       | 3/10 [00:00<00:00, 510.82 it/sec, obj=-0.267]
    INFO - 08:56:30:     40%|████      | 4/10 [00:00<00:00, 549.80 it/sec, obj=-0.809]
    INFO - 08:56:30:     50%|█████     | 5/10 [00:00<00:00, 573.51 it/sec, obj=-0.868]
    INFO - 08:56:30:     60%|██████    | 6/10 [00:00<00:00, 595.84 it/sec, obj=-0.874]
    INFO - 08:56:30:     70%|███████   | 7/10 [00:00<00:00, 569.07 it/sec, obj=-0.266]
    INFO - 08:56:30:     80%|████████  | 8/10 [00:00<00:00, 551.46 it/sec, obj=0.135]
    INFO - 08:56:30:     90%|█████████ | 9/10 [00:00<00:00, 536.02 it/sec, obj=0.577]
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 549.66 it/sec, obj=0.288]
 WARNING - 08:56:30: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 12
    INFO - 08:56:30:    Solution:
 WARNING - 08:56:30:       The solution is not feasible.
    INFO - 08:56:30:       Objective: 0.28811314244587316
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = 0.15144075010180735
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |       Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | -1.252370379421268 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+--------------------+-------------+-------+
    INFO - 08:56:30: *** End MDOScenario execution (time: 0:00:00.028951) ***
    INFO - 08:56:30:    100%|██████████| 10/10 [00:00<00:00, 31.38 it/sec, obj=0.288]
    INFO - 08:56:30: Optimization result:
    INFO - 08:56:30:    Optimizer info:
    INFO - 08:56:30:       Status: None
    INFO - 08:56:30:       Message: None
    INFO - 08:56:30:       Number of calls to the objective function by the optimizer: 10
    INFO - 08:56:30:    Solution:
    INFO - 08:56:30:       The solution is feasible.
    INFO - 08:56:30:       Objective: 0.6150998205402495
    INFO - 08:56:30:       Standardized constraints:
    INFO - 08:56:30:          cstr = -0.7883188793606977
    INFO - 08:56:30:       Design space:
    INFO - 08:56:30:          +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30:          | Name | Lower bound |        Value        | Upper bound | Type  |
    INFO - 08:56:30:          +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30:          | x    |     -1.5    | -0.1666666666666667 |     1.5     | float |
    INFO - 08:56:30:          +------+-------------+---------------------+-------------+-------+
    INFO - 08:56:30: *** End DOEScenario execution (time: 0:00:00.329708) ***

{'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)
History plot
<gemseo.post.basic_history.BasicHistory object at 0x7f1db15ce2b0>

Total running time of the script: (0 minutes 0.520 seconds)

Gallery generated by Sphinx-Gallery