Multistart optimization

Runs simple optimization problem with multiple starting points Nests a MDOScenario in a DOEScenario using a MDOScenarioAdapter.

from __future__ import division, unicode_literals

from matplotlib import pyplot as plt

from gemseo.api import (
    configure_logger,
    create_design_space,
    create_discipline,
    create_scenario,
)
from gemseo.core.mdo_scenario import MDOScenarioAdapter

configure_logger()

Out:

<RootLogger root (INFO)>

Create the disciplines

objective = create_discipline(
    "AnalyticDiscipline", expressions_dict={"obj": "x**3-x+1"}
)
constraint = create_discipline(
    "AnalyticDiscipline", expressions_dict={"cstr": "x**2+obj**2-1.5"}
)

Create the design space

design_space = create_design_space()
design_space.add_variable("x", 1, 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.variables_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)

Out:

    INFO - 12:57:13:
    INFO - 12:57:13: *** Start DOE Scenario execution ***
    INFO - 12:57:13: DOEScenario
    INFO - 12:57:13:    Disciplines: MDOScenario_adapter
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: fullfact
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Full factorial design required. Number of samples along each direction for a design vector of size 1 with 10 samples: 10
    INFO - 12:57:13: Final number of samples for DOE = 10 vs 10 requested
    INFO - 12:57:13: DOE sampling:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound | value | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    |  -1.5 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
 WARNING - 12:57:13: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 521.73 it/sec, obj=1.25]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.2888129873625884
    INFO - 12:57:13: The result is not feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = 0.1513713660745195
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -1.252181466244097 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.024347 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -1.166666666666667 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
 WARNING - 12:57:13: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 623.58 it/sec, obj=0.288]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.28811314244670916
    INFO - 12:57:13: The result is not feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = 0.15144075010172386
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -1.252370379421043 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.020701 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |        value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -0.8333333333333334 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
 WARNING - 12:57:13: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 670.45 it/sec, obj=-.875]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.5837798751488654
    INFO - 12:57:13: The result is not feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = 0.19806414484529755
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -1.165017254042391 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.019606 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound | value | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    |  -0.5 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 513.38 it/sec, obj=0.615]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.6150998219543254
    INFO - 12:57:13: The result is feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 13
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = -0.788285881879476
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.5773788419679762 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.024194 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |        value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -0.1666666666666667 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 498.57 it/sec, obj=0.615]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.6150998205402495
    INFO - 12:57:13: The result is feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = -0.7883188793606977
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.5773502675245377 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.024801 ***
    INFO - 12:57:13: DOE sampling:  50%|█████     | 5/10 [00:00<00:00, 82.09 it/sec]
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.1666666666666667 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
    INFO - 12:57:13: Optimization:  90%|█████████ | 9/10 [00:00<00:00, 578.41 it/sec, obj=0.615]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.6150998205402508
    INFO - 12:57:13: The result is feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 10
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = -0.7883189092934597
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.5773502416020033 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.021990 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound | value | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    |  0.5  |     1.5     | float |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
    INFO - 12:57:13: Optimization:  70%|███████   | 7/10 [00:00<00:00, 726.26 it/sec, obj=0.615]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.6150998205402526
    INFO - 12:57:13: The result is feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 8
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = -0.7883188285028959
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.5773503115686811 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.018634 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.8333333333333335 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 627.36 it/sec, obj=0.288]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.7453703703703706
    INFO - 12:57:13: The result is feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = -0.2499785665294918
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 0.8333333333333335 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.020701 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+-------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value       | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+-------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | 1.166666666666667 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+-------------------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
 WARNING - 12:57:13: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 625.99 it/sec, obj=0.288]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.28811314244698893
    INFO - 12:57:13: The result is not feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = 0.1514407501016959
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -1.252370379420967 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.020660 ***
    INFO - 12:57:13:
    INFO - 12:57:13: *** Start MDO Scenario execution ***
    INFO - 12:57:13: MDOScenario
    INFO - 12:57:13:    Disciplines: AnalyticDiscipline AnalyticDiscipline
    INFO - 12:57:13:    MDOFormulation: DisciplinaryOpt
    INFO - 12:57:13:    Algorithm: SLSQP
    INFO - 12:57:13: Optimization problem:
    INFO - 12:57:13:    Minimize: obj(x)
    INFO - 12:57:13:    With respect to: x
    INFO - 12:57:13:    Subject to constraints:
    INFO - 12:57:13:       cstr(x) <= 0.0
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound | value | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    |  1.5  |     1.5     | float |
    INFO - 12:57:13: +------+-------------+-------+-------------+-------+
    INFO - 12:57:13: Optimization:   0%|          | 0/10 [00:00<?, ?it]
    INFO - 12:57:13: 
 WARNING - 12:57:13: Optimization found no feasible point !  The least infeasible point is selected.
    INFO - 12:57:13: Optimization: 100%|██████████| 10/10 [00:00<00:00, 609.57 it/sec, obj=0.288]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.28811314244587316
    INFO - 12:57:13: The result is not feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: Maximum number of iterations reached. GEMSEO Stopped the driver
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 12
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = 0.15144075010180735
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |       value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -1.252370379421268 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+--------------------+-------------+-------+
    INFO - 12:57:13: *** MDO Scenario run terminated in 0:00:00.021354 ***
    INFO - 12:57:13: DOE sampling: 100%|██████████| 10/10 [00:00<00:00, 43.25 it/sec]
    INFO - 12:57:13: DOE sampling: 100%|██████████| 10/10 [00:00<00:00, 43.13 it/sec]
    INFO - 12:57:13: Optimization result:
    INFO - 12:57:13: Objective value = 0.6150998205402495
    INFO - 12:57:13: The result is feasible.
    INFO - 12:57:13: Status: None
    INFO - 12:57:13: Optimizer message: None
    INFO - 12:57:13: Number of calls to the objective function by the optimizer: 10
    INFO - 12:57:13: Constraints values:
    INFO - 12:57:13:    cstr = -0.7883188793606977
    INFO - 12:57:13: Design space:
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: | name | lower_bound |        value        | upper_bound | type  |
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: | x    |     -1.5    | -0.1666666666666667 |     1.5     | float |
    INFO - 12:57:13: +------+-------------+---------------------+-------------+-------+
    INFO - 12:57:13: *** DOE Scenario run terminated ***

{'eval_jac': False, 'algo': 'fullfact', 'n_samples': 10}

Plot the optimum objective for different x0

scenario_doe.post_process("BasicHistory", data_list=["obj"], save=False, show=False)
# Workaround for HTML rendering, instead of ``show=True``
plt.show()
History plot

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

Gallery generated by Sphinx-Gallery