.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/mda/plot_quasinewton_sobieski.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_mda_plot_quasinewton_sobieski.py: Quasi-Newton MDA ================ .. GENERATED FROM PYTHON SOURCE LINES 24-34 .. code-block:: Python from __future__ import annotations from gemseo import configure_logger from gemseo import create_discipline from gemseo import create_mda configure_logger() .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 35-37 Define a way to display results ------------------------------- .. GENERATED FROM PYTHON SOURCE LINES 37-67 .. code-block:: Python def display_result(res, mda_name): """Display coupling and output variables in logger. Args: res: result of MDA mda_name: name of the current MDA """ # names of the coupling variables coupling_names = [ "y_11", "y_12", "y_14", "y_21", "y_23", "y_24", "y_31", "y_32", "y_34", ] for coupling_var in coupling_names: print(f"{mda_name}, coupling variable {coupling_var}: {res[coupling_var]}") # names of the output variables output_names = ["y_1", "y_2", "y_3", "y_4", "g_1", "g_2", "g_3"] for output_name in output_names: print( f"{mda_name}, output variable {output_name}: {res[output_name]}", ) .. GENERATED FROM PYTHON SOURCE LINES 68-73 Create, execute and post-process MDA ------------------------------------ We do not need to specify the inputs, the default inputs of the MDA will be used and computed from the Default inputs of the disciplines .. GENERATED FROM PYTHON SOURCE LINES 73-90 .. code-block:: Python disciplines = create_discipline([ "SobieskiStructure", "SobieskiPropulsion", "SobieskiAerodynamics", "SobieskiMission", ]) mda = create_mda("MDAQuasiNewton", disciplines, method="broyden1") res = mda.execute() display_result(res, mda.name) mda.plot_residual_history( n_iterations=10, logscale=[1e-8, 10.0], save=False, show=True, fig_size=(10, 2), ) .. image-sg:: /examples/mda/images/sphx_glr_plot_quasinewton_sobieski_001.png :alt: MDAQuasiNewton: residual plot :srcset: /examples/mda/images/sphx_glr_plot_quasinewton_sobieski_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none MDAQuasiNewton, coupling variable y_11: [0.15591864] MDAQuasiNewton, coupling variable y_12: [5.06070636e+04 9.50000000e-01] MDAQuasiNewton, coupling variable y_14: [50607.0636194 7306.20262124] MDAQuasiNewton, coupling variable y_21: [50607.06361948] MDAQuasiNewton, coupling variable y_23: [12194.40918574] MDAQuasiNewton, coupling variable y_24: [4.15002177] MDAQuasiNewton, coupling variable y_31: [6354.40199386] MDAQuasiNewton, coupling variable y_32: [0.50275822] MDAQuasiNewton, coupling variable y_34: [1.10754577] MDAQuasiNewton, output variable y_1: [5.06070636e+04 7.30620262e+03 9.50000000e-01] MDAQuasiNewton, output variable y_2: [5.06070636e+04 1.21944092e+04 4.15002177e+00] MDAQuasiNewton, output variable y_3: [1.10754577e+00 6.35440200e+03 5.02802106e-01] MDAQuasiNewton, output variable y_4: [535.78213201] MDAQuasiNewton, output variable g_1: [ 0.035 -0.00666667 -0.0275 -0.04 -0.04833333 -0.09 -0.15 ] MDAQuasiNewton, output variable g_2: [-0.04] MDAQuasiNewton, output variable g_3: [-0.99719789 -0.00280211 0.16206032 -0.02 ] INFO - 10:56:14: Requested 10 iterations but the residual history contains only 8, plotting all the residual history.
.. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 0.283 seconds) .. _sphx_glr_download_examples_mda_plot_quasinewton_sobieski.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_quasinewton_sobieski.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_quasinewton_sobieski.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_