.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/formulations/plot_sobieski_idf_example.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_formulations_plot_sobieski_idf_example.py: IDF-based MDO on the Sobieski SSBJ test case ============================================ .. GENERATED FROM PYTHON SOURCE LINES 26-36 .. code-block:: default from __future__ import absolute_import, division, print_function, unicode_literals from future import standard_library from gemseo.api import configure_logger, create_discipline, create_scenario from gemseo.problems.sobieski.core import SobieskiProblem standard_library.install_aliases() configure_logger() .. rst-class:: sphx-glr-script-out Out: .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 37-44 Instantiate the disciplines ---------------------------- First, we instantiate the four disciplines of the use case: :class:`~gemseo.problems.sobieski.wrappers.SobieskiPropulsion`, :class:`~gemseo.problems.sobieski.wrappers.SobieskiAerodynamics`, :class:`~gemseo.problems.sobieski.wrappers.SobieskiMission` and :class:`~gemseo.problems.sobieski.wrappers.SobieskiStructure`. .. GENERATED FROM PYTHON SOURCE LINES 44-53 .. code-block:: default disciplines = create_discipline( [ "SobieskiPropulsion", "SobieskiAerodynamics", "SobieskiMission", "SobieskiStructure", ] ) .. GENERATED FROM PYTHON SOURCE LINES 54-63 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:`.IDF` formulation. We tell the scenario to minimize -y_4 instead of minimizing y_4 (range), which is the default option. Instantiate the scenario ^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 63-72 .. code-block:: default design_space = SobieskiProblem().read_design_space() scenario = create_scenario( disciplines, "IDF", objective_name="y_4", design_space=design_space, maximize_objective=True, ) .. GENERATED FROM PYTHON SOURCE LINES 73-75 Set the design constraints ^^^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 75-78 .. code-block:: default for c_name in ["g_1", "g_2", "g_3"]: scenario.add_constraint(c_name, "ineq") .. GENERATED FROM PYTHON SOURCE LINES 79-83 Define the algorithm inputs ^^^^^^^^^^^^^^^^^^^^^^^^^^^ We set the maximum number of iterations, the optimizer and the optimizer options .. GENERATED FROM PYTHON SOURCE LINES 83-91 .. code-block:: default algo_options = { "ftol_rel": 1e-10, "ineq_tolerance": 1e-3, "eq_tolerance": 1e-3, "normalize_design_space": True, } scn_inputs = {"max_iter": 20, "algo": "SLSQP", "algo_options": algo_options} .. GENERATED FROM PYTHON SOURCE LINES 92-94 Execute the scenario ^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 94-96 .. code-block:: default scenario.execute(scn_inputs) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none {'max_iter': 20, 'algo': 'SLSQP', 'algo_options': {'ftol_rel': 1e-10, 'ineq_tolerance': 0.001, 'eq_tolerance': 0.001, 'normalize_design_space': True}} .. GENERATED FROM PYTHON SOURCE LINES 97-101 Save the optimization history ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ We can save the whole optimization problem and its history for further post processing: .. GENERATED FROM PYTHON SOURCE LINES 101-102 .. code-block:: default scenario.save_optimization_history("idf_history.h5", file_format="hdf5") .. GENERATED FROM PYTHON SOURCE LINES 103-104 We can also save only calls to functions and design variables history: .. GENERATED FROM PYTHON SOURCE LINES 104-106 .. code-block:: default scenario.save_optimization_history("idf_history.xml", file_format="ggobi") .. GENERATED FROM PYTHON SOURCE LINES 107-109 Print optimization metrics ^^^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 109-111 .. code-block:: default scenario.print_execution_metrics() .. GENERATED FROM PYTHON SOURCE LINES 112-114 Plot the optimization history view ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 114-116 .. code-block:: default scenario.post_process("OptHistoryView", save=False, show=True) .. rst-class:: sphx-glr-horizontal * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_001.png :alt: Evolution of the optimization variables :class: sphx-glr-multi-img * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_002.png :alt: Evolution of the objective value :class: sphx-glr-multi-img * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_003.png :alt: Distance to the optimum :class: sphx-glr-multi-img * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_004.png :alt: Hessian diagonal approximation :class: sphx-glr-multi-img * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_005.png :alt: Evolution of the inequality constraints :class: sphx-glr-multi-img * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_006.png :alt: Evolution of the equality constraints :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out Out: .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/gemseo/conda/3.0.3/lib/python3.8/site-packages/gemseo/post/opt_history_view.py:312: UserWarning: FixedFormatter should only be used together with FixedLocator ax1.set_yticklabels(y_labels) /home/docs/checkouts/readthedocs.org/user_builds/gemseo/conda/3.0.3/lib/python3.8/site-packages/gemseo/post/opt_history_view.py:716: MatplotlibDeprecationWarning: default base will change from np.e to 10 in 3.4. To suppress this warning specify the base keyword argument. norm=SymLogNorm(linthresh=linthresh, vmin=-vmax, vmax=vmax), /home/docs/checkouts/readthedocs.org/user_builds/gemseo/conda/3.0.3/lib/python3.8/site-packages/gemseo/post/opt_history_view.py:626: MatplotlibDeprecationWarning: default base will change from np.e to 10 in 3.4. To suppress this warning specify the base keyword argument. norm=SymLogNorm(linthresh=1.0, vmin=-vmax, vmax=vmax), /home/docs/checkouts/readthedocs.org/user_builds/gemseo/conda/3.0.3/lib/python3.8/site-packages/gemseo/post/opt_history_view.py:619: MatplotlibDeprecationWarning: Passing parameters norm and vmin/vmax simultaneously is deprecated since 3.3 and will become an error two minor releases later. Please pass vmin/vmax directly to the norm when creating it. im1 = ax1.imshow( .. GENERATED FROM PYTHON SOURCE LINES 117-119 Plot the quadratic approximation of the objective ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 119-120 .. code-block:: default scenario.post_process("QuadApprox", function="-y_4", save=False, show=True) .. rst-class:: sphx-glr-horizontal * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_007.png :alt: Hessian matrix SR1 approximation of -y_4 :class: sphx-glr-multi-img * .. image:: /examples/formulations/images/sphx_glr_plot_sobieski_idf_example_008.png :alt: plot sobieski idf example :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out Out: .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/gemseo/conda/3.0.3/lib/python3.8/site-packages/gemseo/post/quad_approx.py:151: MatplotlibDeprecationWarning: default base will change from np.e to 10 in 3.4. To suppress this warning specify the base keyword argument. norm=SymLogNorm(linthresh=linthresh, vmin=-vmax, vmax=vmax), .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 2.922 seconds) .. _sphx_glr_download_examples_formulations_plot_sobieski_idf_example.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_sobieski_idf_example.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_sobieski_idf_example.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_