.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/post_process/plot_basic_history.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_post_process_plot_basic_history.py: Basic history ============= In this example, we illustrate the use of the :class:`.BasicHistory` plot on the Sobieski's SSBJ problem. .. GENERATED FROM PYTHON SOURCE LINES 28-34 .. code-block:: default from gemseo.api import configure_logger from gemseo.api import create_discipline from gemseo.api import create_scenario from gemseo.problems.sobieski.core.problem import SobieskiProblem from matplotlib import pyplot as plt .. GENERATED FROM PYTHON SOURCE LINES 35-39 Import ------ The first step is to import some functions from the API and a method to get the design space. .. GENERATED FROM PYTHON SOURCE LINES 39-42 .. code-block:: default configure_logger() .. rst-class:: sphx-glr-script-out Out: .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 43-48 Description ----------- The :class:`~gemseo.post.basic_history.BasicHistory` post-processing plots any of the constraint or objective functions w.r.t. the optimization iterations or sampling snapshots. .. GENERATED FROM PYTHON SOURCE LINES 50-54 Create disciplines ------------------ At this point, we instantiate the disciplines of Sobieski's SSBJ problem: Propulsion, Aerodynamics, Structure and Mission .. GENERATED FROM PYTHON SOURCE LINES 54-63 .. code-block:: default disciplines = create_discipline( [ "SobieskiPropulsion", "SobieskiAerodynamics", "SobieskiStructure", "SobieskiMission", ] ) .. GENERATED FROM PYTHON SOURCE LINES 64-67 Create design space ------------------- We also read the design space from the :class:`.SobieskiProblem`. .. GENERATED FROM PYTHON SOURCE LINES 67-69 .. code-block:: default design_space = SobieskiProblem().design_space .. GENERATED FROM PYTHON SOURCE LINES 70-77 Create and execute scenario --------------------------- The next step is to build an MDO scenario in order to maximize the range, encoded 'y_4', with respect to the design parameters, while satisfying the inequality constraints 'g_1', 'g_2' and 'g_3'. We can use the MDF formulation, the SLSQP optimization algorithm and a maximum number of iterations equal to 100. .. GENERATED FROM PYTHON SOURCE LINES 77-89 .. code-block:: default scenario = create_scenario( disciplines, formulation="MDF", objective_name="y_4", maximize_objective=True, design_space=design_space, ) scenario.set_differentiation_method("user") for constraint in ["g_1", "g_2", "g_3"]: scenario.add_constraint(constraint, "ineq") scenario.execute({"algo": "SLSQP", "max_iter": 10}) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none INFO - 07:14:21: INFO - 07:14:21: *** Start MDOScenario execution *** INFO - 07:14:21: MDOScenario INFO - 07:14:21: Disciplines: SobieskiPropulsion SobieskiAerodynamics SobieskiStructure SobieskiMission INFO - 07:14:21: MDO formulation: MDF INFO - 07:14:21: Optimization problem: INFO - 07:14:21: minimize -y_4(x_shared, x_1, x_2, x_3) INFO - 07:14:21: with respect to x_1, x_2, x_3, x_shared INFO - 07:14:21: subject to constraints: INFO - 07:14:21: g_1(x_shared, x_1, x_2, x_3) <= 0.0 INFO - 07:14:21: g_2(x_shared, x_1, x_2, x_3) <= 0.0 INFO - 07:14:21: g_3(x_shared, x_1, x_2, x_3) <= 0.0 INFO - 07:14:21: over the design space: INFO - 07:14:21: +----------+-------------+-------+-------------+-------+ INFO - 07:14:21: | name | lower_bound | value | upper_bound | type | INFO - 07:14:21: +----------+-------------+-------+-------------+-------+ INFO - 07:14:21: | x_shared | 0.01 | 0.05 | 0.09 | float | INFO - 07:14:21: | x_shared | 30000 | 45000 | 60000 | float | INFO - 07:14:21: | x_shared | 1.4 | 1.6 | 1.8 | float | INFO - 07:14:21: | x_shared | 2.5 | 5.5 | 8.5 | float | INFO - 07:14:21: | x_shared | 40 | 55 | 70 | float | INFO - 07:14:21: | x_shared | 500 | 1000 | 1500 | float | INFO - 07:14:21: | x_1 | 0.1 | 0.25 | 0.4 | float | INFO - 07:14:21: | x_1 | 0.75 | 1 | 1.25 | float | INFO - 07:14:21: | x_2 | 0.75 | 1 | 1.25 | float | INFO - 07:14:21: | x_3 | 0.1 | 0.5 | 1 | float | INFO - 07:14:21: +----------+-------------+-------+-------------+-------+ INFO - 07:14:21: Solving optimization problem with algorithm SLSQP: INFO - 07:14:21: ... 0%| | 0/10 [00:00 .. GENERATED FROM PYTHON SOURCE LINES 110-115 .. note:: Set the boolean instance attribute :attr:`.OptimizationProblem.use_standardized_objective` to ``False`` to plot the objective to maximize as a performance function. .. GENERATED FROM PYTHON SOURCE LINES 115-119 .. code-block:: default scenario.use_standardized_objective = False scenario.post_process("BasicHistory", variable_names=["y_4"], save=False, show=False) # Workaround for HTML rendering, instead of ``show=True`` plt.show() .. image-sg:: /examples/post_process/images/sphx_glr_plot_basic_history_003.png :alt: History plot :srcset: /examples/post_process/images/sphx_glr_plot_basic_history_003.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 1.964 seconds) .. _sphx_glr_download_examples_post_process_plot_basic_history.py: .. only :: html .. container:: sphx-glr-footer :class: sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_basic_history.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_basic_history.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_