.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/formulations/plot_doe_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_doe_sobieski_mdf_example.py>`
        to download the full example code

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

.. _sphx_glr_examples_formulations_plot_doe_sobieski_mdf_example.py:


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

.. GENERATED FROM PYTHON SOURCE LINES 24-36

.. code-block:: default

    from __future__ import annotations

    from os import name as os_name

    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 37-44

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 44-53

.. code-block:: default

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








.. GENERATED FROM PYTHON SOURCE LINES 54-57

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 57-61

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




.. GENERATED FROM PYTHON SOURCE LINES 62-66

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 66-68

.. code-block:: default

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




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





.. GENERATED FROM PYTHON SOURCE LINES 69-77

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:`.BiLevel` formulation. We tell the scenario to minimize -y_4
instead of minimizing y_4 (range), which is the default option.

We need to define the design space.

.. GENERATED FROM PYTHON SOURCE LINES 77-80

.. code-block:: default

    design_space = SobieskiProblem().design_space
    print(design_space)





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

 .. code-block:: none

    Design space: 
    +-------------+-------------+--------------------+-------------+-------+
    | name        | lower_bound |       value        | upper_bound | type  |
    +-------------+-------------+--------------------+-------------+-------+
    | x_shared[0] |     0.01    |        0.05        |     0.09    | float |
    | x_shared[1] |    30000    |       45000        |    60000    | float |
    | x_shared[2] |     1.4     |        1.6         |     1.8     | float |
    | x_shared[3] |     2.5     |        5.5         |     8.5     | float |
    | x_shared[4] |      40     |         55         |      70     | float |
    | x_shared[5] |     500     |        1000        |     1500    | float |
    | x_1[0]      |     0.1     |        0.25        |     0.4     | float |
    | x_1[1]      |     0.75    |         1          |     1.25    | float |
    | x_2         |     0.75    |         1          |     1.25    | float |
    | x_3         |     0.1     |        0.5         |      1      | float |
    | y_14[0]     |    24850    |   50606.9741711    |    77100    | float |
    | y_14[1]     |    -7700    |   7306.20262124    |    45000    | float |
    | y_32        |    0.235    | 0.5027962499999999 |    0.795    | float |
    | y_31        |     2960    |   6354.32430691    |    10185    | float |
    | y_24        |     0.44    |     4.15006276     |    11.13    | float |
    | y_34        |     0.44    |     1.10754577     |     1.98    | float |
    | y_23        |     3365    |   12194.2671934    |    26400    | float |
    | y_21        |    24850    |   50606.9741711    |    77250    | float |
    | y_12[0]     |    24850    |     50606.9742     |    77250    | float |
    | y_12[1]     |     0.45    |        0.95        |     1.5     | float |
    +-------------+-------------+--------------------+-------------+-------+




.. GENERATED FROM PYTHON SOURCE LINES 81-83

Instantiate the scenario
^^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 83-92

.. code-block:: default

    scenario = create_scenario(
        disciplines,
        formulation="MDF",
        objective_name="y_4",
        design_space=design_space,
        maximize_objective=True,
        scenario_type="DOE",
    )








.. GENERATED FROM PYTHON SOURCE LINES 93-95

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

.. GENERATED FROM PYTHON SOURCE LINES 95-98

.. code-block:: default

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








.. GENERATED FROM PYTHON SOURCE LINES 99-106

Visualize the XDSM
^^^^^^^^^^^^^^^^^^
Generate the XDSM 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 106-108

.. 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_2, x_1, x_shared, x_3"
          },
          {
            "from": "Dis1",
            "to": "Opt",
            "name": "g_2, g_1, g_3, y_4"
          },
          {
            "from": "Opt",
            "to": "Dis2",
            "name": "x_2, x_1, x_shared, x_3"
          },
          {
            "from": "Dis2",
            "to": "Opt",
            "name": "g_2, g_1, g_3"
          },
          {
            "from": "Opt",
            "to": "Dis3",
            "name": "x_shared, x_3"
          },
          {
            "from": "Dis3",
            "to": "Opt",
            "name": "g_3"
          },
          {
            "from": "Opt",
            "to": "Dis4",
            "name": "x_2, x_shared"
          },
          {
            "from": "Dis4",
            "to": "Opt",
            "name": "g_2"
          },
          {
            "from": "Opt",
            "to": "Dis5",
            "name": "x_shared, x_1"
          },
          {
            "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 109-112

Execute the scenario
^^^^^^^^^^^^^^^^^^^^
Use provided analytic derivatives

.. GENERATED FROM PYTHON SOURCE LINES 112-114

.. code-block:: default

    scenario.set_differentiation_method()








.. GENERATED FROM PYTHON SOURCE LINES 115-120

Multiprocessing
^^^^^^^^^^^^^^^
It is possible to run a DOE in parallel using multiprocessing, in order to do
this, we specify the number of processes to be used for the computation of
the samples.

.. GENERATED FROM PYTHON SOURCE LINES 122-129

.. warning::
   The multiprocessing option has some limitations on Windows.
   Due to problems with sphinx, we disable it in this example.
   The features :class:`.MemoryFullCache` and :class:`.HDF5Cache` are not
   available for multiprocessing on Windows.
   As an alternative, we recommend the method
   :meth:`.DOEScenario.set_optimization_history_backup`.

.. GENERATED FROM PYTHON SOURCE LINES 131-133

We define the algorithm options. Here the criterion = center option of pyDOE
centers the points within the sampling intervals.

.. GENERATED FROM PYTHON SOURCE LINES 133-144

.. code-block:: default

    algo_options = {
        "criterion": "center",
        # Evaluate gradient of the MDA
        # with coupled adjoint
        "eval_jac": True,
        # Run in parallel on 1 or 4 processors
        "n_processes": 1 if os_name == "nt" else 4,
    }

    scenario.execute({"n_samples": 30, "algo": "lhs", "algo_options": algo_options})





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

 .. code-block:: none

        INFO - 16:26:32:  
        INFO - 16:26:32: *** Start DOEScenario execution ***
        INFO - 16:26:32: DOEScenario
        INFO - 16:26:32:    Disciplines: SobieskiAerodynamics SobieskiMission SobieskiPropulsion SobieskiStructure
        INFO - 16:26:32:    MDO formulation: MDF
        INFO - 16:26:32: Optimization problem:
        INFO - 16:26:32:    minimize -y_4(x_shared, x_1, x_2, x_3)
        INFO - 16:26:32:    with respect to x_1, x_2, x_3, x_shared
        INFO - 16:26:32:    subject to constraints:
        INFO - 16:26:32:       g_1(x_shared, x_1, x_2, x_3) <= 0.0
        INFO - 16:26:32:       g_2(x_shared, x_1, x_2, x_3) <= 0.0
        INFO - 16:26:32:       g_3(x_shared, x_1, x_2, x_3) <= 0.0
        INFO - 16:26:32:    over the design space:
        INFO - 16:26:32:    +-------------+-------------+-------+-------------+-------+
        INFO - 16:26:32:    | name        | lower_bound | value | upper_bound | type  |
        INFO - 16:26:32:    +-------------+-------------+-------+-------------+-------+
        INFO - 16:26:32:    | x_shared[0] |     0.01    |  0.05 |     0.09    | float |
        INFO - 16:26:32:    | x_shared[1] |    30000    | 45000 |    60000    | float |
        INFO - 16:26:32:    | x_shared[2] |     1.4     |  1.6  |     1.8     | float |
        INFO - 16:26:32:    | x_shared[3] |     2.5     |  5.5  |     8.5     | float |
        INFO - 16:26:32:    | x_shared[4] |      40     |   55  |      70     | float |
        INFO - 16:26:32:    | x_shared[5] |     500     |  1000 |     1500    | float |
        INFO - 16:26:32:    | x_1[0]      |     0.1     |  0.25 |     0.4     | float |
        INFO - 16:26:32:    | x_1[1]      |     0.75    |   1   |     1.25    | float |
        INFO - 16:26:32:    | x_2         |     0.75    |   1   |     1.25    | float |
        INFO - 16:26:32:    | x_3         |     0.1     |  0.5  |      1      | float |
        INFO - 16:26:32:    +-------------+-------------+-------+-------------+-------+
        INFO - 16:26:32: Solving optimization problem with algorithm lhs:
        INFO - 16:26:32: ...   0%|          | 0/30 [00:00<?, ?it]
        INFO - 16:26:32: Running DOE in parallel on n_processes = 4
        INFO - 16:26:33: ...   3%|▎         | 1/30 [00:00<00:26,  1.11 it/sec, obj=-285]
        INFO - 16:26:33: ...   7%|▋         | 2/30 [00:00<00:13,  2.14 it/sec, obj=-344]
        INFO - 16:26:33: ...  10%|█         | 3/30 [00:00<00:08,  3.18 it/sec, obj=-567]
        INFO - 16:26:33: ...  13%|█▎        | 4/30 [00:00<00:06,  4.22 it/sec, obj=-222]
        INFO - 16:26:33: ...  17%|█▋        | 5/30 [00:01<00:06,  3.73 it/sec, obj=-481]
        INFO - 16:26:33: ...  20%|██        | 6/30 [00:01<00:05,  4.02 it/sec, obj=-429]
        INFO - 16:26:33: ...  23%|██▎       | 7/30 [00:01<00:05,  4.57 it/sec, obj=-519]
        INFO - 16:26:33: ...  27%|██▋       | 8/30 [00:01<00:04,  5.14 it/sec, obj=-235]
        INFO - 16:26:34: ...  30%|███       | 9/30 [00:02<00:04,  4.49 it/sec, obj=-538]
        INFO - 16:26:34: ...  33%|███▎      | 10/30 [00:02<00:04,  4.85 it/sec, obj=-306]
        INFO - 16:26:34: ...  37%|███▋      | 11/30 [00:02<00:03,  5.33 it/sec, obj=-556]
        INFO - 16:26:34: ...  40%|████      | 12/30 [00:02<00:03,  5.75 it/sec, obj=-419]
        INFO - 16:26:34: ...  43%|████▎     | 13/30 [00:02<00:03,  5.07 it/sec, obj=-495]
        INFO - 16:26:34: ...  47%|████▋     | 14/30 [00:02<00:02,  5.43 it/sec, obj=-433]
        INFO - 16:26:34: ...  50%|█████     | 15/30 [00:02<00:02,  5.78 it/sec, obj=-621]
        INFO - 16:26:34: ...  53%|█████▎    | 16/30 [00:02<00:02,  6.09 it/sec, obj=-574]
        INFO - 16:26:35: ...  57%|█████▋    | 17/30 [00:03<00:02,  5.65 it/sec, obj=-481]
        INFO - 16:26:35: ...  60%|██████    | 18/30 [00:03<00:02,  5.77 it/sec, obj=-1.07e+3]
        INFO - 16:26:35: ...  63%|██████▎   | 19/30 [00:03<00:01,  6.02 it/sec, obj=-287]
        INFO - 16:26:35: ...  67%|██████▋   | 20/30 [00:03<00:01,  6.32 it/sec, obj=-602]
        INFO - 16:26:35: ...  70%|███████   | 21/30 [00:03<00:01,  5.80 it/sec, obj=-247]
        INFO - 16:26:35: ...  73%|███████▎  | 22/30 [00:03<00:01,  5.99 it/sec, obj=-273]
        INFO - 16:26:35: ...  77%|███████▋  | 23/30 [00:03<00:01,  6.23 it/sec, obj=-414]
        INFO - 16:26:35: ...  80%|████████  | 24/30 [00:03<00:00,  6.45 it/sec, obj=-1.18e+3]
        INFO - 16:26:36: ...  83%|████████▎ | 25/30 [00:04<00:00,  6.08 it/sec, obj=-383]
        INFO - 16:26:36: ...  87%|████████▋ | 26/30 [00:04<00:00,  6.16 it/sec, obj=-624]
        INFO - 16:26:36: ...  90%|█████████ | 27/30 [00:04<00:00,  6.32 it/sec, obj=-606]
        INFO - 16:26:36: ...  93%|█████████▎| 28/30 [00:04<00:00,  6.54 it/sec, obj=-351]
        INFO - 16:26:36: ...  97%|█████████▋| 29/30 [00:04<00:00,  6.42 it/sec, obj=-485]
        INFO - 16:26:36: ... 100%|██████████| 30/30 [00:04<00:00,  6.63 it/sec, obj=-405]
        INFO - 16:26:36: Optimization result:
        INFO - 16:26:36:    Optimizer info:
        INFO - 16:26:36:       Status: None
        INFO - 16:26:36:       Message: None
        INFO - 16:26:36:       Number of calls to the objective function by the optimizer: 30
        INFO - 16:26:36:    Solution:
        INFO - 16:26:36:       The solution is feasible.
        INFO - 16:26:36:       Objective: -485.49220045924955
        INFO - 16:26:36:       Standardized constraints:
        INFO - 16:26:36:          g_1 = [-0.11350951 -0.10812292 -0.1045109  -0.10204971 -0.10028641 -0.01838903
        INFO - 16:26:36:  -0.22161097]
        INFO - 16:26:36:          g_2 = -0.02400000000000002
        INFO - 16:26:36:          g_3 = [-0.33063157 -0.66936843 -0.73821755 -0.07789536]
        INFO - 16:26:36:       Design space: 
        INFO - 16:26:36:       +-------------+-------------+---------------------+-------------+-------+
        INFO - 16:26:36:       | name        | lower_bound |        value        | upper_bound | type  |
        INFO - 16:26:36:       +-------------+-------------+---------------------+-------------+-------+
        INFO - 16:26:36:       | x_shared[0] |     0.01    | 0.05400000000000001 |     0.09    | float |
        INFO - 16:26:36:       | x_shared[1] |    30000    |        46500        |    60000    | float |
        INFO - 16:26:36:       | x_shared[2] |     1.4     |  1.686666666666667  |     1.8     | float |
        INFO - 16:26:36:       | x_shared[3] |     2.5     |         5.2         |     8.5     | float |
        INFO - 16:26:36:       | x_shared[4] |      40     |         66.5        |      70     | float |
        INFO - 16:26:36:       | x_shared[5] |     500     |  583.3333333333334  |     1500    | float |
        INFO - 16:26:36:       | x_1[0]      |     0.1     |        0.185        |     0.4     | float |
        INFO - 16:26:36:       | x_1[1]      |     0.75    |  0.9416666666666667 |     1.25    | float |
        INFO - 16:26:36:       | x_2         |     0.75    |        0.775        |     1.25    | float |
        INFO - 16:26:36:       | x_3         |     0.1     |        0.115        |      1      | float |
        INFO - 16:26:36:       +-------------+-------------+---------------------+-------------+-------+
        INFO - 16:26:36: *** End DOEScenario execution (time: 0:00:04.601650) ***

    {'eval_jac': False, 'algo_options': {'criterion': 'center', 'eval_jac': True, 'n_processes': 4}, 'n_samples': 30, 'algo': 'lhs'}



.. GENERATED FROM PYTHON SOURCE LINES 145-151

.. warning::
   On Windows, the progress bar may show duplicated instances during the
   initialization of each subprocess. In some cases it may also print the
   conclusion of an iteration ahead of another one that was concluded first.
   This is a consequence of the pickling process and does not affect the
   computations of the scenario.

.. GENERATED FROM PYTHON SOURCE LINES 153-158

Exporting the problem data.
^^^^^^^^^^^^^^^^^^^^^^^^^^^
After the execution of the scenario, you may want to export your data to use it
elsewhere. The method :meth:`.Scenario.to_dataset` will allow you to export
your results to a :class:`.Dataset`, the basic |g| class to store data.

.. GENERATED FROM PYTHON SOURCE LINES 158-160

.. code-block:: default

    dataset = scenario.to_dataset("a_name_for_my_dataset")








.. GENERATED FROM PYTHON SOURCE LINES 161-163

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

.. GENERATED FROM PYTHON SOURCE LINES 163-165

.. code-block:: default

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




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


    *

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

    *

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

    *

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

    *

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

    *

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


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

 .. code-block:: none

    /home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.0.0/lib/python3.9/site-packages/genson/schema/strategies/base.py:42: UserWarning: Schema incompatible. Keyword 'description' has conflicting values ('The width and height of the figure in inches, e.g. ``(w, h)``.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.' vs. 'The width and height of the figure in inches, e.g. `(w, h)`.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.'). Using 'The width and height of the figure in inches, e.g. ``(w, h)``.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.'
      warn(('Schema incompatible. Keyword {0!r} has conflicting '

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



.. GENERATED FROM PYTHON SOURCE LINES 166-173

.. tip::

   Each post-processing method requires different inputs and offers a variety
   of customization options. Use the high-level function
   :func:`.get_post_processing_options_schema` to print a table with
   the attributes for any post-processing algo. Or refer to our dedicated page:
   :ref:`gen_post_algos`.

.. GENERATED FROM PYTHON SOURCE LINES 175-177

Plot the scatter matrix
^^^^^^^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 177-184

.. code-block:: default

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




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


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

 .. code-block:: none


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



.. GENERATED FROM PYTHON SOURCE LINES 185-187

Plot correlations
^^^^^^^^^^^^^^^^^

.. GENERATED FROM PYTHON SOURCE LINES 187-188

.. code-block:: default

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



.. image-sg:: /examples/formulations/images/sphx_glr_plot_doe_sobieski_mdf_example_008.png
   :alt: R=0.98987, R=0.97511, R=0.99671, R=0.96235, R=0.99119, R=0.99866, R=0.95205, R=0.98584, R=0.99618, R=0.99936
   :srcset: /examples/formulations/images/sphx_glr_plot_doe_sobieski_mdf_example_008.png
   :class: sphx-glr-single-img


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

 .. code-block:: none

    /home/docs/checkouts/readthedocs.org/user_builds/gemseo/envs/5.0.0/lib/python3.9/site-packages/genson/schema/strategies/base.py:42: UserWarning: Schema incompatible. Keyword 'description' has conflicting values ('The width and height of the figure in inches, e.g. ``(w, h)``.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.' vs. 'The width and height of the figure in inches, e.g. `(w, h)`.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.'). Using 'The width and height of the figure in inches, e.g. ``(w, h)``.\nIf ``None``, use the :attr:`.OptPostProcessor.DEFAULT_FIG_SIZE`\nof the post-processor.'
      warn(('Schema incompatible. Keyword {0!r} has conflicting '
        INFO - 16:26:40: Detected 10 correlations > 0.95

    <gemseo.post.correlations.Correlations object at 0x7fea43210eb0>




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

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


.. _sphx_glr_download_examples_formulations_plot_doe_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_doe_sobieski_mdf_example.py <plot_doe_sobieski_mdf_example.py>`

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

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


.. only:: html

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

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