.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/formulations/plot_sobieski_mdf_example.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_formulations_plot_sobieski_mdf_example.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_formulations_plot_sobieski_mdf_example.py:


MDF-based MDO on the Sobieski SSBJ test case
============================================

.. GENERATED FROM PYTHON SOURCE LINES 24-34

.. code-block:: default

    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.problems.sobieski.core.problem import SobieskiProblem

    configure_logger()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <RootLogger root (INFO)>



.. GENERATED FROM PYTHON SOURCE LINES 35-42

Instantiate the disciplines
---------------------------
First, we instantiate the four disciplines of the use case:
:class:`.SobieskiPropulsion`,
:class:`.SobieskiAerodynamics`,
:class:`.SobieskiMission`
and :class:`.SobieskiStructure`.

.. GENERATED FROM PYTHON SOURCE LINES 42-51

.. code-block:: default

    disciplines = create_discipline(
        [
            "SobieskiPropulsion",
            "SobieskiAerodynamics",
            "SobieskiMission",
            "SobieskiStructure",
        ]
    )








.. GENERATED FROM PYTHON SOURCE LINES 52-55

We can quickly access the most relevant information of any discipline (name, inputs,
and outputs) with Python's ``print()`` function. Moreover, we can get the default
input values of a discipline with the attribute :attr:`.MDODiscipline.default_inputs`

.. GENERATED FROM PYTHON SOURCE LINES 55-59

.. code-block:: default

    for discipline in disciplines:
        print(discipline)
        print(f"Default inputs: {discipline.default_inputs}")





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    SobieskiPropulsion
    Default inputs: {'x_3': array([0.5]), 'c_3': array([4360.]), 'y_23': array([12562.01206488]), 'x_shared': array([5.0e-02, 4.5e+04, 1.6e+00, 5.5e+00, 5.5e+01, 1.0e+03])}
    SobieskiAerodynamics
    Default inputs: {'y_32': array([0.50279625]), 'x_2': array([1.]), 'y_12': array([5.06069742e+04, 9.50000000e-01]), 'c_4': array([0.01375]), 'x_shared': array([5.0e-02, 4.5e+04, 1.6e+00, 5.5e+00, 5.5e+01, 1.0e+03])}
    SobieskiMission
    Default inputs: {'y_24': array([4.15006276]), 'x_shared': array([5.0e-02, 4.5e+04, 1.6e+00, 5.5e+00, 5.5e+01, 1.0e+03]), 'y_34': array([1.10754577]), 'y_14': array([50606.9741711 ,  7306.20262124])}
    SobieskiStructure
    Default inputs: {'c_1': array([25000.]), 'c_0': array([2000.]), 'c_2': array([6.]), 'y_31': array([6354.32430691]), 'y_21': array([50606.9741711]), 'x_1': array([0.25, 1.  ]), 'x_shared': array([5.0e-02, 4.5e+04, 1.6e+00, 5.5e+00, 5.5e+01, 1.0e+03])}




.. GENERATED FROM PYTHON SOURCE LINES 60-64

You may also be interested in plotting the couplings of your disciplines.
A quick way of getting this information is the high-level function
:func:`.generate_n2_plot`. A much more detailed explanation of coupling
visualization is available :ref:`here <coupling_visualization>`.

.. GENERATED FROM PYTHON SOURCE LINES 64-66

.. code-block:: default

    generate_n2_plot(disciplines, save=False, show=True)




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_001.png
   :alt: plot sobieski mdf example
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_001.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 67-78

Build, execute and post-process the scenario
--------------------------------------------
Then, we build the scenario which links the disciplines
with the formulation and the optimization algorithm. Here, we use the
:class:`.MDF` formulation. We tell the scenario to minimize -y_4 instead of
minimizing y_4 (range), which is the default option.

Instantiate the scenario
^^^^^^^^^^^^^^^^^^^^^^^^
During the instantiation of the scenario, we provide some options for the
MDF formulations:

.. GENERATED FROM PYTHON SOURCE LINES 78-86

.. code-block:: default

    formulation_options = {
        "tolerance": 1e-14,
        "max_mda_iter": 50,
        "warm_start": True,
        "use_lu_fact": False,
        "linear_solver_tolerance": 1e-14,
    }








.. GENERATED FROM PYTHON SOURCE LINES 87-96

- ``'warm_start``: warm starts MDA,
- ``'warm_start``: optimize the adjoints resolution by storing
  the Jacobian matrix LU factorization for the multiple RHS
  (objective + constraints). This saves CPU time if you can pay for
  the memory and have the full Jacobians available, not just matrix vector
  products.
- ``'linear_solver_tolerance'``: set the linear solver tolerance,
  idem we need full convergence


.. GENERATED FROM PYTHON SOURCE LINES 97-100

.. code-block:: default

    design_space = SobieskiProblem().design_space
    design_space






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div style='background-color: #fafae2; padding: 10px;'>Design space:<br/><table>
        <tr>
            <th style='text-align: left;'>name</th>
            <th style='text-align: left;'>lower_bound</th>
            <th style='text-align: left;'>value</th>
            <th style='text-align: left;'>upper_bound</th>
            <th style='text-align: left;'>type</th>
        </tr>
        <tr>
            <td>x_shared[0]</td>
            <td>0.01</td>
            <td>0.05</td>
            <td>0.09</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_shared[1]</td>
            <td>30000</td>
            <td>45000</td>
            <td>60000</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_shared[2]</td>
            <td>1.4</td>
            <td>1.6</td>
            <td>1.8</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_shared[3]</td>
            <td>2.5</td>
            <td>5.5</td>
            <td>8.5</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_shared[4]</td>
            <td>40</td>
            <td>55</td>
            <td>70</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_shared[5]</td>
            <td>500</td>
            <td>1000</td>
            <td>1500</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_1[0]</td>
            <td>0.1</td>
            <td>0.25</td>
            <td>0.4</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_1[1]</td>
            <td>0.75</td>
            <td>1</td>
            <td>1.25</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_2</td>
            <td>0.75</td>
            <td>1</td>
            <td>1.25</td>
            <td>float</td>
        </tr>
        <tr>
            <td>x_3</td>
            <td>0.1</td>
            <td>0.5</td>
            <td>1</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_14[0]</td>
            <td>24850</td>
            <td>50606.9741711</td>
            <td>77100</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_14[1]</td>
            <td>-7700</td>
            <td>7306.20262124</td>
            <td>45000</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_32</td>
            <td>0.235</td>
            <td>0.5027962499999999</td>
            <td>0.795</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_31</td>
            <td>2960</td>
            <td>6354.32430691</td>
            <td>10185</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_24</td>
            <td>0.44</td>
            <td>4.15006276</td>
            <td>11.13</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_34</td>
            <td>0.44</td>
            <td>1.10754577</td>
            <td>1.98</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_23</td>
            <td>3365</td>
            <td>12194.2671934</td>
            <td>26400</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_21</td>
            <td>24850</td>
            <td>50606.9741711</td>
            <td>77250</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_12[0]</td>
            <td>24850</td>
            <td>50606.9742</td>
            <td>77250</td>
            <td>float</td>
        </tr>
        <tr>
            <td>y_12[1]</td>
            <td>0.45</td>
            <td>0.95</td>
            <td>1.5</td>
            <td>float</td>
        </tr>
    </table></div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 101-110

.. code-block:: default

    scenario = create_scenario(
        disciplines,
        "MDF",
        objective_name="y_4",
        design_space=design_space,
        maximize_objective=True,
        **formulation_options,
    )








.. GENERATED FROM PYTHON SOURCE LINES 111-113

Set the design constraints
^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 113-116

.. code-block:: default

    for c_name in ["g_1", "g_2", "g_3"]:
        scenario.add_constraint(c_name, "ineq")








.. GENERATED FROM PYTHON SOURCE LINES 117-124

XDSMIZE the scenario
^^^^^^^^^^^^^^^^^^^^
Generate the XDSM file on the fly:

- ``log_workflow_status=True`` will log the status of the workflow  in the console,
- ``save_html`` (default ``True``) will generate a self-contained HTML file,
  that can be automatically opened using ``show_html=True``.

.. GENERATED FROM PYTHON SOURCE LINES 124-126

.. code-block:: default

    scenario.xdsmize(save_html=False)






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

        <script type="text/javascript">
            document.addEventListener('DOMContentLoaded', () => {
              const mdo = {
      "root": {
        "nodes": [
          {
            "id": "Opt",
            "name": "Optimizer",
            "type": "optimization"
          },
          {
            "id": "Dis1",
            "name": "MDAChain",
            "type": "mda"
          },
          {
            "id": "Dis2",
            "name": "MDAJacobi",
            "type": "mda"
          },
          {
            "id": "Dis3",
            "name": "SobieskiPropulsion",
            "type": "analysis"
          },
          {
            "id": "Dis4",
            "name": "SobieskiAerodynamics",
            "type": "analysis"
          },
          {
            "id": "Dis5",
            "name": "SobieskiStructure",
            "type": "analysis"
          },
          {
            "id": "Dis6",
            "name": "SobieskiMission",
            "type": "analysis"
          }
        ],
        "edges": [
          {
            "from": "_U_",
            "to": "Opt",
            "name": "x_shared^(0), x_1^(0), x_2^(0), x_3^(0)"
          },
          {
            "from": "Opt",
            "to": "_U_",
            "name": "-y_4^*, g_1^*, g_2^*, g_3^*"
          },
          {
            "from": "Opt",
            "to": "Dis1",
            "name": "x_3, x_shared, x_1, x_2"
          },
          {
            "from": "Dis1",
            "to": "Opt",
            "name": "y_4, g_3, g_2, g_1"
          },
          {
            "from": "Opt",
            "to": "Dis2",
            "name": "x_3, x_shared, x_1, x_2"
          },
          {
            "from": "Dis2",
            "to": "Opt",
            "name": "g_3, g_2, g_1"
          },
          {
            "from": "Opt",
            "to": "Dis3",
            "name": "x_3, x_shared"
          },
          {
            "from": "Dis3",
            "to": "Opt",
            "name": "g_3"
          },
          {
            "from": "Opt",
            "to": "Dis4",
            "name": "x_shared, x_2"
          },
          {
            "from": "Dis4",
            "to": "Opt",
            "name": "g_2"
          },
          {
            "from": "Opt",
            "to": "Dis5",
            "name": "x_1, x_shared"
          },
          {
            "from": "Dis5",
            "to": "Opt",
            "name": "g_1"
          },
          {
            "from": "Opt",
            "to": "Dis6",
            "name": "x_shared"
          },
          {
            "from": "Dis6",
            "to": "Opt",
            "name": "y_4"
          },
          {
            "from": "Dis3",
            "to": "_U_",
            "name": "y_32^*, y_31^*, y_3^*, y_34^*"
          },
          {
            "from": "Dis3",
            "to": "Opt",
            "name": "g_3"
          },
          {
            "from": "Dis4",
            "to": "_U_",
            "name": "y_21^*, y_23^*, y_24^*, y_2^*"
          },
          {
            "from": "Dis4",
            "to": "Opt",
            "name": "g_2"
          },
          {
            "from": "Dis5",
            "to": "_U_",
            "name": "y_1^*, y_11^*, y_14^*, y_12^*"
          },
          {
            "from": "Dis5",
            "to": "Opt",
            "name": "g_1"
          },
          {
            "from": "Dis6",
            "to": "Opt",
            "name": "y_4"
          },
          {
            "from": "Dis2",
            "to": "Dis6",
            "name": "y_14, y_24, y_34"
          },
          {
            "from": "Dis2",
            "to": "Dis3",
            "name": "y_23"
          },
          {
            "from": "Dis2",
            "to": "Dis4",
            "name": "y_12, y_32"
          },
          {
            "from": "Dis2",
            "to": "Dis5",
            "name": "y_21, y_31"
          },
          {
            "from": "Dis3",
            "to": "Dis2",
            "name": "y_31, y_32"
          },
          {
            "from": "Dis3",
            "to": "Dis4",
            "name": "y_32"
          },
          {
            "from": "Dis3",
            "to": "Dis5",
            "name": "y_31"
          },
          {
            "from": "Dis4",
            "to": "Dis2",
            "name": "y_21, y_23"
          },
          {
            "from": "Dis4",
            "to": "Dis3",
            "name": "y_23"
          },
          {
            "from": "Dis4",
            "to": "Dis5",
            "name": "y_21"
          },
          {
            "from": "Dis5",
            "to": "Dis2",
            "name": "y_12"
          },
          {
            "from": "Dis5",
            "to": "Dis4",
            "name": "y_12"
          }
        ],
        "workflow": [
          "_U_",
          [
            "Opt",
            [
              "Dis1",
              "Dis2",
              [
                {
                  "parallel": [
                    "Dis3",
                    "Dis4",
                    "Dis5"
                  ]
                }
              ],
              "Dis6"
            ]
          ]
        ],
        "optpb": "Optimization problem:\n   minimize -y_4(x_shared, x_1, x_2, x_3)\n   with respect to x_1, x_2, x_3, x_shared\n   subject to constraints:\n      g_1(x_shared, x_1, x_2, x_3) <= 0.0\n      g_2(x_shared, x_1, x_2, x_3) <= 0.0\n      g_3(x_shared, x_1, x_2, x_3) <= 0.0"
      }
    };
              const config = {
                labelizer: {
                    ellipsis: 5,
                    subSupScript: false,
                    showLinkNbOnly: false,
                },
                layout: {
                    origin: { x: 50, y: 20 },
                    cellsize: { w: 150, h: 50 },
                    padding: 10,
                },
                withDefaultDriver: true,
                withTitleTooltip: true,
              };
              xdsmjs.XDSMjs(config).createXdsm(mdo);
            });
        </script>
        <div class='xdsm-toolbar'></div><div class='xdsm2'></div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 127-138

Define the algorithm inputs
^^^^^^^^^^^^^^^^^^^^^^^^^^^
We set the maximum number of iterations, the optimizer
and the optimizer options. Algorithm specific options are passed there.
Use the high-level function :func:`.get_algorithm_options_schema`
for more information or read the documentation.

Here ftol_rel option is a stop criteria based on the relative difference
in the objective between two iterates ineq_tolerance the tolerance
determination of the optimum; this is specific to the |g| wrapping and not
in the solver.

.. GENERATED FROM PYTHON SOURCE LINES 138-145

.. code-block:: default

    algo_options = {
        "ftol_rel": 1e-10,
        "ineq_tolerance": 2e-3,
        "normalize_design_space": True,
    }
    scn_inputs = {"max_iter": 15, "algo": "SLSQP", "algo_options": algo_options}








.. GENERATED FROM PYTHON SOURCE LINES 146-160

.. seealso::

   We can also generate a backup file for the optimization,
   as well as plots on the fly of the optimization history if option
   ``generate_opt_plot`` is ``True``.
   This slows down a lot the process, here since SSBJ is very light

   .. code::

    scenario.set_optimization_history_backup(file_path="mdf_backup.h5",
                                             each_new_iter=True,
                                             each_store=False, erase=True,
                                             pre_load=False,
                                             generate_opt_plot=True)

.. GENERATED FROM PYTHON SOURCE LINES 162-164

Execute the scenario
^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 164-166

.. code-block:: default

    scenario.execute(scn_inputs)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

        INFO - 08:24:45:  
        INFO - 08:24:45: *** Start MDOScenario execution ***
        INFO - 08:24:45: MDOScenario
        INFO - 08:24:45:    Disciplines: SobieskiAerodynamics SobieskiMission SobieskiPropulsion SobieskiStructure
        INFO - 08:24:45:    MDO formulation: MDF
        INFO - 08:24:45: Optimization problem:
        INFO - 08:24:45:    minimize -y_4(x_shared, x_1, x_2, x_3)
        INFO - 08:24:45:    with respect to x_1, x_2, x_3, x_shared
        INFO - 08:24:45:    subject to constraints:
        INFO - 08:24:45:       g_1(x_shared, x_1, x_2, x_3) <= 0.0
        INFO - 08:24:45:       g_2(x_shared, x_1, x_2, x_3) <= 0.0
        INFO - 08:24:45:       g_3(x_shared, x_1, x_2, x_3) <= 0.0
        INFO - 08:24:45:    over the design space:
        INFO - 08:24:45:    +-------------+-------------+-------+-------------+-------+
        INFO - 08:24:45:    | name        | lower_bound | value | upper_bound | type  |
        INFO - 08:24:45:    +-------------+-------------+-------+-------------+-------+
        INFO - 08:24:45:    | x_shared[0] |     0.01    |  0.05 |     0.09    | float |
        INFO - 08:24:45:    | x_shared[1] |    30000    | 45000 |    60000    | float |
        INFO - 08:24:45:    | x_shared[2] |     1.4     |  1.6  |     1.8     | float |
        INFO - 08:24:45:    | x_shared[3] |     2.5     |  5.5  |     8.5     | float |
        INFO - 08:24:45:    | x_shared[4] |      40     |   55  |      70     | float |
        INFO - 08:24:45:    | x_shared[5] |     500     |  1000 |     1500    | float |
        INFO - 08:24:45:    | x_1[0]      |     0.1     |  0.25 |     0.4     | float |
        INFO - 08:24:45:    | x_1[1]      |     0.75    |   1   |     1.25    | float |
        INFO - 08:24:45:    | x_2         |     0.75    |   1   |     1.25    | float |
        INFO - 08:24:45:    | x_3         |     0.1     |  0.5  |      1      | float |
        INFO - 08:24:45:    +-------------+-------------+-------+-------------+-------+
        INFO - 08:24:45: Solving optimization problem with algorithm SLSQP:
        INFO - 08:24:45: ...   0%|          | 0/15 [00:00<?, ?it]
        INFO - 08:24:45: ...   7%|▋         | 1/15 [00:00<00:01,  8.57 it/sec, obj=-536]
        INFO - 08:24:45: ...  13%|█▎        | 2/15 [00:00<00:01,  6.54 it/sec, obj=-2.12e+3]
        INFO - 08:24:46: ...  20%|██        | 3/15 [00:00<00:02,  5.33 it/sec, obj=-3.72e+3]
        INFO - 08:24:46: ...  27%|██▋       | 4/15 [00:00<00:02,  5.35 it/sec, obj=-3.97e+3]
        INFO - 08:24:46: Optimization result:
        INFO - 08:24:46:    Optimizer info:
        INFO - 08:24:46:       Status: 8
        INFO - 08:24:46:       Message: Positive directional derivative for linesearch
        INFO - 08:24:46:       Number of calls to the objective function by the optimizer: 5
        INFO - 08:24:46:    Solution:
        INFO - 08:24:46:       The solution is feasible.
        INFO - 08:24:46:       Objective: -3716.555963095829
        INFO - 08:24:46:       Standardized constraints:
        INFO - 08:24:46:          g_1 = [-0.01608807 -0.03194613 -0.04316738 -0.05095364 -0.05658344 -0.1380806
        INFO - 08:24:46:  -0.1019194 ]
        INFO - 08:24:46:          g_2 = -0.0005956359157315294
        INFO - 08:24:46:          g_3 = [-0.67076432 -0.32923568 -0.10429595 -0.183255  ]
        INFO - 08:24:46:       Design space:
        INFO - 08:24:46:       +-------------+-------------+---------------------+-------------+-------+
        INFO - 08:24:46:       | name        | lower_bound |        value        | upper_bound | type  |
        INFO - 08:24:46:       +-------------+-------------+---------------------+-------------+-------+
        INFO - 08:24:46:       | x_shared[0] |     0.01    | 0.05985109102106711 |     0.09    | float |
        INFO - 08:24:46:       | x_shared[1] |    30000    |   59785.5639131558  |    60000    | float |
        INFO - 08:24:46:       | x_shared[2] |     1.4     |         1.4         |     1.8     | float |
        INFO - 08:24:46:       | x_shared[3] |     2.5     |  2.540129655171779  |     8.5     | float |
        INFO - 08:24:46:       | x_shared[4] |      40     |  69.80684214977607  |      70     | float |
        INFO - 08:24:46:       | x_shared[5] |     500     |  1493.746505655324  |     1500    | float |
        INFO - 08:24:46:       | x_1[0]      |     0.1     |         0.4         |     0.4     | float |
        INFO - 08:24:46:       | x_1[1]      |     0.75    |  0.7530970468499044 |     1.25    | float |
        INFO - 08:24:46:       | x_2         |     0.75    |  0.7530625218813826 |     1.25    | float |
        INFO - 08:24:46:       | x_3         |     0.1     |  0.1411034879427379 |      1      | float |
        INFO - 08:24:46:       +-------------+-------------+---------------------+-------------+-------+
        INFO - 08:24:46: *** End MDOScenario execution (time: 0:00:00.860172) ***

    {'max_iter': 15, 'algo_options': {'ftol_rel': 1e-10, 'ineq_tolerance': 0.002, 'normalize_design_space': True}, 'algo': 'SLSQP'}



.. GENERATED FROM PYTHON SOURCE LINES 167-171

Save the optimization history
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
We can save the whole optimization problem and its history for further post
processing:

.. GENERATED FROM PYTHON SOURCE LINES 171-173

.. code-block:: default

    scenario.save_optimization_history("mdf_history.h5", file_format="hdf5")





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

        INFO - 08:24:46: Export optimization problem to file: mdf_history.h5




.. GENERATED FROM PYTHON SOURCE LINES 174-175

We can also save only calls to functions and design variables history:

.. GENERATED FROM PYTHON SOURCE LINES 175-177

.. code-block:: default

    scenario.save_optimization_history("mdf_history.xml", file_format="ggobi")





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

        INFO - 08:24:46: Export to ggobi for functions: ['-y_4', 'g_1', 'g_2', 'g_3']
        INFO - 08:24:46: Export to ggobi file: mdf_history.xml




.. GENERATED FROM PYTHON SOURCE LINES 178-180

Print optimization metrics
^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 180-182

.. code-block:: default

    scenario.print_execution_metrics()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

        INFO - 08:24:46: Scenario Execution Statistics
        INFO - 08:24:46:    Discipline: SobieskiPropulsion
        INFO - 08:24:46:       Executions number: 67
        INFO - 08:24:46:       Execution time: 0.026514842003052763 s
        INFO - 08:24:46:       Linearizations number: 4
        INFO - 08:24:46:    Discipline: SobieskiAerodynamics
        INFO - 08:24:46:       Executions number: 76
        INFO - 08:24:46:       Execution time: 0.038043101001676405 s
        INFO - 08:24:46:       Linearizations number: 4
        INFO - 08:24:46:    Discipline: SobieskiMission
        INFO - 08:24:46:       Executions number: 4
        INFO - 08:24:46:       Execution time: 0.00036684400129161077 s
        INFO - 08:24:46:       Linearizations number: 4
        INFO - 08:24:46:    Discipline: SobieskiStructure
        INFO - 08:24:46:       Executions number: 75
        INFO - 08:24:46:       Execution time: 0.20256956800221815 s
        INFO - 08:24:46:       Linearizations number: 4
        INFO - 08:24:46:    Total number of executions calls: 222
        INFO - 08:24:46:    Total number of linearizations: 16




.. GENERATED FROM PYTHON SOURCE LINES 183-188

Post-process the results
------------------------

Plot the optimization history view
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 188-190

.. code-block:: default

    scenario.post_process("OptHistoryView", save=False, show=True)




.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_002.png
         :alt: Evolution of the optimization variables
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_002.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_003.png
         :alt: Evolution of the objective value
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_003.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_004.png
         :alt: Distance to the optimum
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_004.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_005.png
         :alt: Hessian diagonal approximation
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_005.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_006.png
         :alt: Evolution of the inequality constraints
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_006.png
         :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.opt_history_view.OptHistoryView object at 0x7f0cbc969eb0>



.. GENERATED FROM PYTHON SOURCE LINES 191-193

Plot the basic history view
^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 193-197

.. code-block:: default

    scenario.post_process(
        "BasicHistory", variable_names=["x_shared"], save=False, show=True
    )




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_007.png
   :alt: History plot
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_007.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.basic_history.BasicHistory object at 0x7f0cb8542040>



.. GENERATED FROM PYTHON SOURCE LINES 198-200

Plot the constraints and objective history
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 200-202

.. code-block:: default

    scenario.post_process("ObjConstrHist", save=False, show=True)




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_008.png
   :alt: Evolution of the objective and maximum constraint
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_008.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.obj_constr_hist.ObjConstrHist object at 0x7f0cb811be50>



.. GENERATED FROM PYTHON SOURCE LINES 203-205

Plot the constraints history
^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 205-212

.. code-block:: default

    scenario.post_process(
        "ConstraintsHistory",
        constraint_names=["g_1", "g_2", "g_3"],
        save=False,
        show=True,
    )




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_009.png
   :alt: Evolution of the constraints w.r.t. iterations, g_1[0] (inequality), g_1[1] (inequality), g_1[2] (inequality), g_1[3] (inequality), g_1[4] (inequality), g_1[5] (inequality), g_1[6] (inequality), g_2 (inequality), g_3[0] (inequality), g_3[1] (inequality), g_3[2] (inequality), g_3[3] (inequality)
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_009.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.constraints_history.ConstraintsHistory object at 0x7f0cb7f365e0>



.. GENERATED FROM PYTHON SOURCE LINES 213-215

Plot the constraints history using a radar chart
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 215-222

.. code-block:: default

    scenario.post_process(
        "RadarChart",
        constraint_names=["g_1", "g_2", "g_3"],
        save=False,
        show=True,
    )




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_010.png
   :alt: Constraints at iteration 2 (optimum)
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_010.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.radar_chart.RadarChart object at 0x7f0cb8388f10>



.. GENERATED FROM PYTHON SOURCE LINES 223-225

Plot the quadratic approximation of the objective
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 225-227

.. code-block:: default

    scenario.post_process("QuadApprox", function="-y_4", save=False, show=True)




.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_011.png
         :alt: Hessian matrix SR1 approximation of -y_4
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_011.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_012.png
         :alt: plot sobieski mdf example
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_012.png
         :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.quad_approx.QuadApprox object at 0x7f0cb79c1340>



.. GENERATED FROM PYTHON SOURCE LINES 228-230

Plot the functions using a SOM
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 230-232

.. code-block:: default

    scenario.post_process("SOM", save=False, show=True)




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_013.png
   :alt: Self Organizing Maps of the design space, -y_4, g_1_0, g_1_1, g_1_2, g_1_3, g_1_4, g_1_5, g_1_6, g_2, g_3_0, g_3_1, g_3_2, g_3_3
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_013.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

        INFO - 08:24:49: Building Self Organizing Map from optimization history:
        INFO - 08:24:49:     Number of neurons in x direction = 4
        INFO - 08:24:49:     Number of neurons in y direction = 4

    <gemseo.post.som.SOM object at 0x7f0cb79a3e80>



.. GENERATED FROM PYTHON SOURCE LINES 233-235

Plot the scatter matrix of variables of interest
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 235-243

.. code-block:: default

    scenario.post_process(
        "ScatterPlotMatrix",
        variable_names=["-y_4", "g_1"],
        save=False,
        show=True,
        fig_size=(14, 14),
    )




.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_014.png
   :alt: plot sobieski mdf example
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_014.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.scatter_mat.ScatterPlotMatrix object at 0x7f0cb6d5dbe0>



.. GENERATED FROM PYTHON SOURCE LINES 244-246

Plot the variables using the parallel coordinates
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 246-248

.. code-block:: default

    scenario.post_process("ParallelCoordinates", save=False, show=True)




.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_015.png
         :alt: Design variables history colored by '-y_4' value
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_015.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_016.png
         :alt: Objective function and constraints history colored by '-y_4' value.
         :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_016.png
         :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.para_coord.ParallelCoordinates object at 0x7f0cb6e3c160>



.. GENERATED FROM PYTHON SOURCE LINES 249-251

Plot the robustness of the solution
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 251-253

.. code-block:: default

    scenario.post_process("Robustness", save=True, show=True)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gemseo.post.robustness.Robustness object at 0x7f0cb58c27f0>



.. GENERATED FROM PYTHON SOURCE LINES 254-256

Plot the influence of the design variables
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 256-257

.. code-block:: default

    scenario.post_process("VariableInfluence", fig_size=(14, 14), save=False, show=True)



.. image-sg:: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_017.png
   :alt: Partial variation of the functions wrt design variables, 9 variables required to explain 99% of -y_4 variations, 5 variables required to explain 99% of g_1[0] variations, 5 variables required to explain 99% of g_1[1] variations, 5 variables required to explain 99% of g_1[2] variations, 5 variables required to explain 99% of g_1[3] variations, 5 variables required to explain 99% of g_1[4] variations, 4 variables required to explain 99% of g_1[5] variations, 4 variables required to explain 99% of g_1[6] variations, 1 variables required to explain 99% of g_2 variations, 7 variables required to explain 99% of g_3[0] variations, 7 variables required to explain 99% of g_3[1] variations, 3 variables required to explain 99% of g_3[2] variations, 3 variables required to explain 99% of g_3[3] variations
   :srcset: /examples/formulations/images/sphx_glr_plot_sobieski_mdf_example_017.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

        INFO - 08:24:55: Output name; most influential variables to explain 0.99% of the output variation 
        INFO - 08:24:55:    -y_4; x_1[1], x_2, x_3, x_shared[0], x_shared[1], x_shared[2], x_shared[3], x_shared[4], x_shared[5]
    /home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.1.1/lib/python3.9/site-packages/gemseo/post/variable_influence.py:238: UserWarning: set_ticklabels() should only be used with a fixed number of ticks, i.e. after set_ticks() or using a FixedLocator.
      axe.set_xticklabels(x_labels, fontsize=font_size, rotation=rotation)
        INFO - 08:24:55:    g_1[0]; x_1[0], x_1[1], x_shared[0], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_1[1]; x_1[0], x_1[1], x_shared[0], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_1[2]; x_1[0], x_1[1], x_shared[0], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_1[3]; x_1[0], x_1[1], x_shared[0], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_1[4]; x_1[0], x_1[1], x_shared[0], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_1[5]; x_1[0], x_1[1], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_1[6]; x_1[0], x_1[1], x_shared[3], x_shared[5]
        INFO - 08:24:55:    g_2; x_shared[0]
        INFO - 08:24:55:    g_3[0]; x_2, x_3, x_shared[0], x_shared[1], x_shared[2], x_shared[4], x_shared[5]
        INFO - 08:24:55:    g_3[1]; x_2, x_3, x_shared[0], x_shared[1], x_shared[2], x_shared[4], x_shared[5]
        INFO - 08:24:55:    g_3[2]; x_3, x_shared[1], x_shared[2]
        INFO - 08:24:55:    g_3[3]; x_3, x_shared[1], x_shared[2]

    <gemseo.post.variable_influence.VariableInfluence object at 0x7f0cb59e0d00>




.. rst-class:: sphx-glr-timing

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


.. _sphx_glr_download_examples_formulations_plot_sobieski_mdf_example.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example




    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: plot_sobieski_mdf_example.py <plot_sobieski_mdf_example.py>`

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: plot_sobieski_mdf_example.ipynb <plot_sobieski_mdf_example.ipynb>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_