Note
Go to the end to download the full example code.
MDO formulations for a toy example in aerostructure#
from __future__ import annotations
from gemseo import configure_logger
from gemseo import create_discipline
from gemseo import create_scenario
from gemseo import generate_n2_plot
from gemseo.algos.opt.nlopt.settings.nlopt_cobyla_settings import NLOPT_COBYLA_Settings
from gemseo.algos.opt.nlopt.settings.nlopt_slsqp_settings import NLOPT_SLSQP_Settings
from gemseo.problems.mdo.aerostructure.aerostructure_design_space import (
AerostructureDesignSpace,
)
configure_logger()
# Passed to algo settings
cobyla_settings = NLOPT_COBYLA_Settings(
max_iter=7,
xtol_rel=1e-8,
xtol_abs=1e-8,
ftol_rel=1e-8,
ftol_abs=1e-8,
ineq_tolerance=1e-5,
eq_tolerance=1e-3,
)
slsqp_settings = NLOPT_SLSQP_Settings(
max_iter=10,
xtol_rel=1e-8,
xtol_abs=1e-8,
ftol_rel=1e-8,
ftol_abs=1e-8,
ineq_tolerance=1e-5,
eq_tolerance=1e-3,
)
Create discipline#
First, we create disciplines (aero, structure, mission) with dummy formulas
using the AnalyticDiscipline class.
aero_formulas = {
"drag": "0.1*((sweep/360)**2 + 200 + thick_airfoils**2-thick_airfoils -4*displ)",
"forces": "10*sweep + 0.2*thick_airfoils-0.2*displ",
"lift": "(sweep + 0.2*thick_airfoils-2.*displ)/3000.",
}
aerodynamics = create_discipline(
"AnalyticDiscipline", name="Aerodynamics", expressions=aero_formulas
)
struc_formulas = {
"mass": "4000*(sweep/360)**3 + 200000 + 100*thick_panels +200.0*forces",
"reserve_fact": "-3*sweep -6*thick_panels+0.1*forces+55",
"displ": "2*sweep + 3*thick_panels-2.*forces",
}
structure = create_discipline(
"AnalyticDiscipline", name="Structure", expressions=struc_formulas
)
mission_formulas = {"range": "8e11*lift/(mass*drag)"}
mission = create_discipline(
"AnalyticDiscipline", name="Mission", expressions=mission_formulas
)
disciplines = [aerodynamics, structure, mission]
We can see that structure and aerodynamics are strongly coupled:
generate_n2_plot(disciplines, save=False, show=True)
Create an MDO scenario with MDF formulation#
Then, we create an MDO scenario based on the MDF formulation
design_space = AerostructureDesignSpace()
scenario = create_scenario(
disciplines,
"range",
design_space,
maximize_objective=True,
formulation_name="MDF",
)
scenario.add_constraint("reserve_fact", constraint_type="ineq", value=0.5)
scenario.add_constraint("lift", value=0.5)
scenario.execute(slsqp_settings)
scenario.post_process(post_name="OptHistoryView", save=False, show=True)
/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/design_space.py:424: ComplexWarning: Casting complex values to real discards the imaginary part
self.__current_value[name] = array_value.astype(
INFO - 08:36:05: Variable reserve_fact was removed from the Design Space, it is not an input of any discipline.
INFO - 08:36:05:
INFO - 08:36:05: *** Start MDOScenario execution ***
INFO - 08:36:05: MDOScenario
INFO - 08:36:05: Disciplines: Aerodynamics Mission Structure
INFO - 08:36:05: MDO formulation: MDF
INFO - 08:36:05: Optimization problem:
INFO - 08:36:05: minimize -range(thick_airfoils, thick_panels, sweep)
INFO - 08:36:05: with respect to sweep, thick_airfoils, thick_panels
INFO - 08:36:05: subject to constraints:
INFO - 08:36:05: reserve_fact(thick_airfoils, thick_panels, sweep) <= 0.5
INFO - 08:36:05: lift(thick_airfoils, thick_panels, sweep) == 0.5
INFO - 08:36:05: over the design space:
INFO - 08:36:05: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:05: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:05: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:05: | thick_airfoils | 5 | 15 | 25 | float |
INFO - 08:36:05: | thick_panels | 1 | 3 | 20 | float |
INFO - 08:36:05: | sweep | 10 | 25 | 35 | float |
INFO - 08:36:05: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:05: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:05: 10%|█ | 1/10 [00:00<00:00, 80.23 it/sec, obj=-4.25e+3]
INFO - 08:36:05: 20%|██ | 2/10 [00:00<00:00, 28.71 it/sec, obj=-4.51e+3]
INFO - 08:36:05: 30%|███ | 3/10 [00:00<00:00, 33.14 it/sec, obj=-4.51e+3]
INFO - 08:36:05: 40%|████ | 4/10 [00:00<00:00, 43.82 it/sec, obj=Not evaluated]
INFO - 08:36:05: Optimization result:
INFO - 08:36:05: Optimizer info:
INFO - 08:36:05: Status: None
INFO - 08:36:05: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 08:36:05: Number of calls to the objective function by the optimizer: 4
INFO - 08:36:05: Solution:
INFO - 08:36:05: The solution is feasible.
INFO - 08:36:05: Objective: -4509.505446934416
INFO - 08:36:05: Standardized constraints:
INFO - 08:36:05: [lift-0.5] = 8.881784197001252e-15
INFO - 08:36:05: [reserve_fact-0.5] = 6.464347279688809e-09
INFO - 08:36:05: Design space:
INFO - 08:36:05: +----------------+-------------+-------------------+-------------+-------+
INFO - 08:36:05: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:05: +----------------+-------------+-------------------+-------------+-------+
INFO - 08:36:05: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:05: | thick_panels | 1 | 3.225589224567829 | 20 | float |
INFO - 08:36:05: | sweep | 10 | 24.99326599309897 | 35 | float |
INFO - 08:36:05: +----------------+-------------+-------------------+-------------+-------+
INFO - 08:36:05: *** End MDOScenario execution (time: 0:00:00.098212) ***
<gemseo.post.opt_history_view.OptHistoryView object at 0x7f6dc8c69fd0>
Create an MDO scenario with bilevel formulation#
Then, we create an MDO scenario based on the bilevel formulation
design_space_ref = AerostructureDesignSpace()
/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/design_space.py:424: ComplexWarning: Casting complex values to real discards the imaginary part
self.__current_value[name] = array_value.astype(
Create the aeronautics sub-scenario#
For this purpose, we create a first sub-scenario to maximize the range with respect to the thick airfoils, based on the aerodynamics discipline.
aero_scenario = create_scenario(
[aerodynamics, mission],
"range",
design_space_ref.filter(["thick_airfoils"], copy=True),
maximize_objective=True,
formulation_name="DisciplinaryOpt",
)
aero_scenario.set_algorithm(slsqp_settings)
Create the structure sub-scenario#
We create a second sub-scenario to maximize the range with respect to the thick panels, based on the structure discipline.
struct_scenario = create_scenario(
[structure, mission],
"range",
design_space_ref.filter(["thick_panels"], copy=True),
maximize_objective=True,
formulation_name="DisciplinaryOpt",
)
struct_scenario.set_algorithm(slsqp_settings)
Create the system scenario#
Lastly, we build a system scenario to maximize the range with respect to the sweep, which is a shared variable, based on the previous sub-scenarios.
design_space_system = design_space_ref.filter(["sweep"], copy=True)
system_scenario = create_scenario(
[aero_scenario, struct_scenario, mission],
"range",
design_space_system,
formulation_name="BiLevel",
maximize_objective=True,
main_mda_settings={"inner_mda_name": "MDAJacobi", "tolerance": 1e-8},
)
system_scenario.add_constraint("reserve_fact", constraint_type="ineq", value=0.5)
system_scenario.add_constraint("lift", value=0.5)
system_scenario.execute(cobyla_settings)
system_scenario.post_process(post_name="OptHistoryView", save=False, show=True)
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: MDOScenario MDOScenario Mission
INFO - 08:36:06: MDO formulation: BiLevel
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(sweep)
INFO - 08:36:06: with respect to sweep
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(sweep) <= 0.5
INFO - 08:36:06: lift(sweep) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +-------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +-------+-------------+-------+-------------+-------+
INFO - 08:36:06: | sweep | 10 | 25 | 35 | float |
INFO - 08:36:06: +-------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_COBYLA:
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 15 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1239.45 it/sec, obj=-4.27e+3]
WARNING - 08:36:06: Optimization found no feasible point; the least infeasible point is selected.
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 326.49 it/sec, obj=-4.51e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: 5
INFO - 08:36:06: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 16
INFO - 08:36:06: Solution:
WARNING - 08:36:06: The solution is not feasible.
INFO - 08:36:06: Objective: -4513.429203824652
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = 0.008666666666666822
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.009325) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1291.35 it/sec, obj=-4.51e+3]
ERROR - 08:36:06: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 427.88 it/sec, obj=-4.5e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 3
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4504.955637332531
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = 1.062133492268913e-09
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.266666666489645 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.007811) ***
INFO - 08:36:06: 14%|█▍ | 1/7 [00:00<00:00, 35.09 it/sec, obj=-4.51e+3]
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1693.30 it/sec, obj=-4.27e+3]
WARNING - 08:36:06: Optimization found no feasible point; the least infeasible point is selected.
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: 5
INFO - 08:36:06: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 16
INFO - 08:36:06: Solution:
WARNING - 08:36:06: The solution is not feasible.
INFO - 08:36:06: Objective: -4267.910320106269
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = 0.12708333333392363
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.006288) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.266666666489645 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 2088.80 it/sec, obj=-4.27e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 438.67 it/sec, obj=-4.27e+3]
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 482.55 it/sec, obj=-4.27e+3]
INFO - 08:36:06: 40%|████ | 4/10 [00:00<00:00, 621.93 it/sec, obj=-4.27e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 5
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4270.047650046897
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = 1.5540514652911952e-08
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 1.764999997412866 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.009533) ***
INFO - 08:36:06: 29%|██▊ | 2/7 [00:00<00:00, 36.02 it/sec, obj=-4.27e+3]
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1760.83 it/sec, obj=-4.53e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 498.14 it/sec, obj=-3.83e+3]
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 707.26 it/sec, obj=-4.43e+3]
INFO - 08:36:06: 40%|████ | 4/10 [00:00<00:00, 757.47 it/sec, obj=-4.5e+3]
INFO - 08:36:06: 50%|█████ | 5/10 [00:00<00:00, 778.71 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 60%|██████ | 6/10 [00:00<00:00, 798.15 it/sec, obj=-4.52e+3]
INFO - 08:36:06: 70%|███████ | 7/10 [00:00<00:00, 814.52 it/sec, obj=-4.52e+3]
INFO - 08:36:06: 80%|████████ | 8/10 [00:00<00:00, 821.99 it/sec, obj=-4.52e+3]
INFO - 08:36:06: 90%|█████████ | 9/10 [00:00<00:00, 832.37 it/sec, obj=-4.53e+3]
INFO - 08:36:06: 100%|██████████| 10/10 [00:00<00:00, 840.04 it/sec, obj=-4.53e+3]
WARNING - 08:36:06: Optimization found no feasible point; the least infeasible point is selected.
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: Maximum number of iterations reached. GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 12
INFO - 08:36:06: Solution:
WARNING - 08:36:06: The solution is not feasible.
INFO - 08:36:06: Objective: -3831.807218275781
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = -0.0023374064594913757
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 25 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.015612) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 1.764999997412866 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1476.87 it/sec, obj=-3.82e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 442.62 it/sec, obj=-3.82e+3]
ERROR - 08:36:06: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 493.72 it/sec, obj=-3.82e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 4
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -3818.5902856031416
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = 2.8144597763457568e-11
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.414591822896584 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.008765) ***
INFO - 08:36:06: 43%|████▎ | 3/7 [00:00<00:00, 32.81 it/sec, obj=-3.82e+3]
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 25 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1763.79 it/sec, obj=-3.82e+3]
WARNING - 08:36:06: Optimization found no feasible point; the least infeasible point is selected.
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 386.39 it/sec, obj=-4.49e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: 5
INFO - 08:36:06: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 16
INFO - 08:36:06: Solution:
WARNING - 08:36:06: The solution is not feasible.
INFO - 08:36:06: Objective: -4490.54017449428
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = 0.00918103190482844
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.007924) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.414591822896584 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1467.57 it/sec, obj=-4.5e+3]
ERROR - 08:36:06: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 489.39 it/sec, obj=-4.5e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 3
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4503.594576633829
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = 1.1323209037072957e-10
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.256301558340416 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.006853) ***
INFO - 08:36:06: 57%|█████▋ | 4/7 [00:00<00:00, 34.27 it/sec, obj=-4.51e+3]
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1510.37 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 459.17 it/sec, obj=-3.85e+3]
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 654.64 it/sec, obj=-4.48e+3]
INFO - 08:36:06: 40%|████ | 4/10 [00:00<00:00, 649.88 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 50%|█████ | 5/10 [00:00<00:00, 663.15 it/sec, obj=-4.51e+3]
ERROR - 08:36:06: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 08:36:06: 60%|██████ | 6/10 [00:00<00:00, 543.89 it/sec, obj=-3.85e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 8
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -3854.0819929120544
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = 0.0
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+------------------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 24.3567703351527 | 25 | float |
INFO - 08:36:06: +----------------+-------------+------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.013700) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.256301558340416 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1470.14 it/sec, obj=-3.84e+3]
ERROR - 08:36:06: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 491.80 it/sec, obj=-3.84e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 3
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -3843.382025541499
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = -1.1290524071227992e-10
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.303233590674794 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.006748) ***
INFO - 08:36:06: 71%|███████▏ | 5/7 [00:00<00:00, 33.84 it/sec, obj=-3.85e+3]
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+------------------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 24.3567703351527 | 25 | float |
INFO - 08:36:06: +----------------+-------------+------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1766.77 it/sec, obj=-3.85e+3]
WARNING - 08:36:06: Optimization found no feasible point; the least infeasible point is selected.
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 366.51 it/sec, obj=-4.49e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: 5
INFO - 08:36:06: Message: NLOPT_MAXEVAL_REACHED: Optimization stopped because maxeval (above) was reached
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 16
INFO - 08:36:06: Solution:
WARNING - 08:36:06: The solution is not feasible.
INFO - 08:36:06: Objective: -4492.38920700084
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = 0.00858587320790527
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.008150) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.303233590674794 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1372.93 it/sec, obj=-4.5e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 430.85 it/sec, obj=-4.5e+3]
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 480.26 it/sec, obj=-4.5e+3]
INFO - 08:36:06: 40%|████ | 4/10 [00:00<00:00, 618.04 it/sec, obj=-4.5e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: Successive iterates of the objective function are closer than ftol_rel or ftol_abs. GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 5
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4504.850158478924
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = 5.3804996014150674e-09
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.264536834765667 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.009656) ***
INFO - 08:36:06: 86%|████████▌ | 6/7 [00:00<00:00, 33.82 it/sec, obj=-4.51e+3]
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Aerodynamics Mission
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_airfoils)
INFO - 08:36:06: with respect to thick_airfoils
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: lift(thick_airfoils) == 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1751.28 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 471.85 it/sec, obj=-4.48e+3]
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 673.24 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 40%|████ | 4/10 [00:00<00:00, 689.40 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 50%|█████ | 5/10 [00:00<00:00, 695.85 it/sec, obj=-4.51e+3]
ERROR - 08:36:06: NLopt run failed: NLopt roundoff-limited, RoundoffLimited
Traceback (most recent call last):
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/gemseo/algos/opt/nlopt/nlopt.py", line 384, in _run
nlopt_problem.optimize(x_0.real)
File "/home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/6.0.0/lib/python3.9/site-packages/nlopt/nlopt.py", line 335, in optimize
return _nlopt.opt_optimize(self, *args)
nlopt.RoundoffLimited: NLopt roundoff-limited
INFO - 08:36:06: 60%|██████ | 6/10 [00:00<00:00, 539.18 it/sec, obj=-4.48e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 8
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4509.519641910682
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = -0.0001369125207841826
INFO - 08:36:06: Design space:
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: | thick_airfoils | 5 | 5 | 25 | float |
INFO - 08:36:06: +----------------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.013876) ***
INFO - 08:36:06:
INFO - 08:36:06: *** Start MDOScenario execution ***
INFO - 08:36:06: MDOScenario
INFO - 08:36:06: Disciplines: Mission Structure
INFO - 08:36:06: MDO formulation: DisciplinaryOpt
INFO - 08:36:06: Optimization problem:
INFO - 08:36:06: minimize -range(thick_panels)
INFO - 08:36:06: with respect to thick_panels
INFO - 08:36:06: subject to constraints:
INFO - 08:36:06: reserve_fact(thick_panels) <= 0.5
INFO - 08:36:06: over the design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.264536834765667 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: Solving optimization problem with algorithm NLOPT_SLSQP:
INFO - 08:36:06: 10%|█ | 1/10 [00:00<00:00, 1702.92 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 20%|██ | 2/10 [00:00<00:00, 452.88 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 30%|███ | 3/10 [00:00<00:00, 466.86 it/sec, obj=-4.51e+3]
INFO - 08:36:06: 40%|████ | 4/10 [00:00<00:00, 566.78 it/sec, obj=Not evaluated]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: Successive iterates of the design variables are closer than xtol_rel or xtol_abs. GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 4
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4509.583287345491
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [reserve_fact-0.5] = -2.105693397425057e-10
INFO - 08:36:06: Design space:
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: | thick_panels | 1 | 3.225021426834291 | 20 | float |
INFO - 08:36:06: +--------------+-------------+-------------------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.009759) ***
INFO - 08:36:06: 100%|██████████| 7/7 [00:00<00:00, 32.87 it/sec, obj=-4.51e+3]
INFO - 08:36:06: Optimization result:
INFO - 08:36:06: Optimizer info:
INFO - 08:36:06: Status: None
INFO - 08:36:06: Message: Maximum number of iterations reached. GEMSEO stopped the driver.
INFO - 08:36:06: Number of calls to the objective function by the optimizer: 9
INFO - 08:36:06: Solution:
INFO - 08:36:06: The solution is feasible.
INFO - 08:36:06: Objective: -4509.243708600623
INFO - 08:36:06: Standardized constraints:
INFO - 08:36:06: [lift-0.5] = 5.903055821931957e-13
INFO - 08:36:06: [reserve_fact-0.5] = -0.2599999989201578
INFO - 08:36:06: Design space:
INFO - 08:36:06: +-------+-------------+-------+-------------+-------+
INFO - 08:36:06: | Name | Lower bound | Value | Upper bound | Type |
INFO - 08:36:06: +-------+-------------+-------+-------------+-------+
INFO - 08:36:06: | sweep | 10 | 25 | 35 | float |
INFO - 08:36:06: +-------+-------------+-------+-------------+-------+
INFO - 08:36:06: *** End MDOScenario execution (time: 0:00:00.217560) ***
<gemseo.post.opt_history_view.OptHistoryView object at 0x7f6df7e94850>
Total running time of the script: (0 minutes 2.887 seconds)