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

.. only:: html

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

        Click :ref:`here <sphx_glr_download_examples_post_process_algorithms_plot_pareto_front_binhkorn.py>`
        to download the full example code

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

.. _sphx_glr_examples_post_process_algorithms_plot_pareto_front_binhkorn.py:


Pareto front on Binh and Korn problem
=====================================

In this example, we illustrate the use of the :class:`.ParetoFront` plot
on the Binh and Korn multi-objective problem.

.. GENERATED FROM PYTHON SOURCE LINES 28-35

.. code-block:: default

    from __future__ import annotations

    from gemseo.algos.doe.doe_factory import DOEFactory
    from gemseo.api import configure_logger
    from gemseo.post.post_factory import PostFactory
    from gemseo.problems.analytical.binh_korn import BinhKorn








.. GENERATED FROM PYTHON SOURCE LINES 36-40

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

.. code-block:: default


    configure_logger()






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

 .. code-block:: none


    <RootLogger root (INFO)>



.. GENERATED FROM PYTHON SOURCE LINES 45-48

Import the optimization problem
-------------------------------
Then, we instantiate the BinkKorn optimization problem.

.. GENERATED FROM PYTHON SOURCE LINES 48-51

.. code-block:: default


    problem = BinhKorn()








.. GENERATED FROM PYTHON SOURCE LINES 52-56

Create and execute scenario
---------------------------
Then, we create a Design of Experiment factory,
and we request the execution a a full-factorial DOE using 100 samples

.. GENERATED FROM PYTHON SOURCE LINES 56-59

.. code-block:: default

    doe_factory = DOEFactory()
    doe_factory.execute(problem, algo_name="OT_OPT_LHS", n_samples=100)





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

 .. code-block:: none

        INFO - 16:55:30: Optimization problem:
        INFO - 16:55:30:    minimize compute_binhkorn(x, y) = (4*x**2+ 4*y**2, (x-5.)**2 + (y-5.)**2)
        INFO - 16:55:30:    with respect to x, y
        INFO - 16:55:30:    subject to constraints:
        INFO - 16:55:30:       ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
        INFO - 16:55:30:       ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
        INFO - 16:55:30:    over the design space:
        INFO - 16:55:30:    +------+-------------+-------+-------------+-------+
        INFO - 16:55:30:    | name | lower_bound | value | upper_bound | type  |
        INFO - 16:55:30:    +------+-------------+-------+-------------+-------+
        INFO - 16:55:30:    | x    |      0      |   1   |      5      | float |
        INFO - 16:55:30:    | y    |      0      |   1   |      3      | float |
        INFO - 16:55:30:    +------+-------------+-------+-------------+-------+
        INFO - 16:55:30: Solving optimization problem with algorithm OT_OPT_LHS:
        INFO - 16:55:30: ...   0%|          | 0/100 [00:00<?, ?it]
        INFO - 16:55:30: ...   1%|          | 1/100 [00:00<00:00, 3127.74 it/sec, obj=[98.50728236 17.96432773]]
        INFO - 16:55:30: ...   2%|▏         | 2/100 [00:00<00:00, 2247.75 it/sec, obj=[ 5.10086961 36.30607322]]
        INFO - 16:55:30: ...   3%|▎         | 3/100 [00:00<00:00, 2155.72 it/sec, obj=[ 5.73344945 34.53104096]]
        INFO - 16:55:30: ...   4%|▍         | 4/100 [00:00<00:00, 2125.85 it/sec, obj=[ 2.99274814 38.89733089]]
        INFO - 16:55:30: ...   5%|▌         | 5/100 [00:00<00:00, 2110.66 it/sec, obj=[ 1.45355254 43.72252794]]
        INFO - 16:55:30: ...   6%|▌         | 6/100 [00:00<00:00, 2038.87 it/sec, obj=[26.6936315  22.26390516]]
        INFO - 16:55:30: ...   7%|▋         | 7/100 [00:00<00:00, 1995.66 it/sec, obj=[95.06326723 19.21865528]]
        INFO - 16:55:30: ...   8%|▊         | 8/100 [00:00<00:00, 1963.51 it/sec, obj=[11.37121324 33.49445422]]
        INFO - 16:55:30: ...   9%|▉         | 9/100 [00:00<00:00, 1969.47 it/sec, obj=[33.52756721 17.72645069]]
        INFO - 16:55:30: ...  10%|█         | 10/100 [00:00<00:00, 1977.33 it/sec, obj=[ 4.39737722 39.0292696 ]]
        INFO - 16:55:30: ...  11%|█         | 11/100 [00:00<00:00, 1983.98 it/sec, obj=[83.81080473  6.44933463]]
        INFO - 16:55:30: ...  12%|█▏        | 12/100 [00:00<00:00, 1990.26 it/sec, obj=[64.12902701 21.25106868]]
        INFO - 16:55:30: ...  13%|█▎        | 13/100 [00:00<00:00, 1994.66 it/sec, obj=[71.23914138  9.90898899]]
        INFO - 16:55:30: ...  14%|█▍        | 14/100 [00:00<00:00, 1970.48 it/sec, obj=[58.61362754 10.51836124]]
        INFO - 16:55:30: ...  15%|█▌        | 15/100 [00:00<00:00, 1955.99 it/sec, obj=[17.1313131  28.27257144]]
        INFO - 16:55:30: ...  16%|█▌        | 16/100 [00:00<00:00, 1946.14 it/sec, obj=[38.07815452 16.26141839]]
        INFO - 16:55:30: ...  17%|█▋        | 17/100 [00:00<00:00, 1951.69 it/sec, obj=[45.29614269 15.87407947]]
        INFO - 16:55:30: ...  18%|█▊        | 18/100 [00:00<00:00, 1957.62 it/sec, obj=[78.26474617 13.97242459]]
        INFO - 16:55:30: ...  19%|█▉        | 19/100 [00:00<00:00, 1963.53 it/sec, obj=[58.47852866 16.87161502]]
        INFO - 16:55:30: ...  20%|██        | 20/100 [00:00<00:00, 1969.30 it/sec, obj=[ 7.27923171 33.31157275]]
        INFO - 16:55:30: ...  21%|██        | 21/100 [00:00<00:00, 1962.93 it/sec, obj=[15.973784   25.86458657]]
        INFO - 16:55:30: ...  22%|██▏       | 22/100 [00:00<00:00, 1952.32 it/sec, obj=[47.66727367 15.96717933]]
        INFO - 16:55:30: ...  23%|██▎       | 23/100 [00:00<00:00, 1940.28 it/sec, obj=[88.93580164  6.7712213 ]]
        INFO - 16:55:30: ...  24%|██▍       | 24/100 [00:00<00:00, 1943.76 it/sec, obj=[22.10070032 24.30208587]]
        INFO - 16:55:30: ...  25%|██▌       | 25/100 [00:00<00:00, 1948.08 it/sec, obj=[16.24176378 31.32352362]]
        INFO - 16:55:30: ...  26%|██▌       | 26/100 [00:00<00:00, 1952.10 it/sec, obj=[29.38105263 25.45614728]]
        INFO - 16:55:30: ...  27%|██▋       | 27/100 [00:00<00:00, 1955.52 it/sec, obj=[27.91051427 19.91664614]]
        INFO - 16:55:30: ...  28%|██▊       | 28/100 [00:00<00:00, 1959.27 it/sec, obj=[37.28915469 18.0103703 ]]
        INFO - 16:55:30: ...  29%|██▉       | 29/100 [00:00<00:00, 1953.53 it/sec, obj=[ 1.01177424 43.45819442]]
        INFO - 16:55:30: ...  30%|███       | 30/100 [00:00<00:00, 1946.40 it/sec, obj=[85.96675181 10.67473104]]
        INFO - 16:55:30: ...  31%|███       | 31/100 [00:00<00:00, 1937.64 it/sec, obj=[87.34794862  6.40751572]]
        INFO - 16:55:30: ...  32%|███▏      | 32/100 [00:00<00:00, 1940.91 it/sec, obj=[25.57943091 20.74700556]]
        INFO - 16:55:30: ...  33%|███▎      | 33/100 [00:00<00:00, 1944.54 it/sec, obj=[76.39229092 11.95231974]]
        INFO - 16:55:30: ...  34%|███▍      | 34/100 [00:00<00:00, 1947.91 it/sec, obj=[51.64657249 21.13171808]]
        INFO - 16:55:30: ...  35%|███▌      | 35/100 [00:00<00:00, 1951.69 it/sec, obj=[95.8196038   9.31540942]]
        INFO - 16:55:30: ...  36%|███▌      | 36/100 [00:00<00:00, 1955.18 it/sec, obj=[117.01496061   6.96792478]]
        INFO - 16:55:30: ...  37%|███▋      | 37/100 [00:00<00:00, 1949.22 it/sec, obj=[78.37371654 21.95951365]]
        INFO - 16:55:30: ...  38%|███▊      | 38/100 [00:00<00:00, 1944.75 it/sec, obj=[38.21112297 20.46923673]]
        INFO - 16:55:31: ...  39%|███▉      | 39/100 [00:00<00:00, 1941.51 it/sec, obj=[77.87964266 17.24338363]]
        INFO - 16:55:31: ...  40%|████      | 40/100 [00:00<00:00, 1944.33 it/sec, obj=[59.79323468 22.48208172]]
        INFO - 16:55:31: ...  41%|████      | 41/100 [00:00<00:00, 1946.89 it/sec, obj=[34.56674439 21.30358144]]
        INFO - 16:55:31: ...  42%|████▏     | 42/100 [00:00<00:00, 1949.72 it/sec, obj=[ 0.46290462 45.68910522]]
        INFO - 16:55:31: ...  43%|████▎     | 43/100 [00:00<00:00, 1952.19 it/sec, obj=[100.73843409  11.99829359]]
        INFO - 16:55:31: ...  44%|████▍     | 44/100 [00:00<00:00, 1954.83 it/sec, obj=[21.18253098 23.10192252]]
        INFO - 16:55:31: ...  45%|████▌     | 45/100 [00:00<00:00, 1948.38 it/sec, obj=[51.0570421  12.66142788]]
        INFO - 16:55:31: ...  46%|████▌     | 46/100 [00:00<00:00, 1944.12 it/sec, obj=[90.5172837  14.59877358]]
        INFO - 16:55:31: ...  47%|████▋     | 47/100 [00:00<00:00, 1942.50 it/sec, obj=[19.5728568  23.86163594]]
        INFO - 16:55:31: ...  48%|████▊     | 48/100 [00:00<00:00, 1944.96 it/sec, obj=[47.52592344 19.80950276]]
        INFO - 16:55:31: ...  49%|████▉     | 49/100 [00:00<00:00, 1947.31 it/sec, obj=[113.47621648   6.0816867 ]]
        INFO - 16:55:31: ...  50%|█████     | 50/100 [00:00<00:00, 1949.70 it/sec, obj=[37.3534796  16.14867587]]
        INFO - 16:55:31: ...  51%|█████     | 51/100 [00:00<00:00, 1951.92 it/sec, obj=[41.42991507 25.64561402]]
        INFO - 16:55:31: ...  52%|█████▏    | 52/100 [00:00<00:00, 1951.48 it/sec, obj=[11.19222105 30.06624599]]
        INFO - 16:55:31: ...  53%|█████▎    | 53/100 [00:00<00:00, 1947.46 it/sec, obj=[86.67968625 15.70216699]]
        INFO - 16:55:31: ...  54%|█████▍    | 54/100 [00:00<00:00, 1943.22 it/sec, obj=[90.96123503  7.41344494]]
        INFO - 16:55:31: ...  55%|█████▌    | 55/100 [00:00<00:00, 1943.53 it/sec, obj=[34.69157131 20.23247304]]
        INFO - 16:55:31: ...  56%|█████▌    | 56/100 [00:00<00:00, 1945.35 it/sec, obj=[56.73716821 11.35918483]]
        INFO - 16:55:31: ...  57%|█████▋    | 57/100 [00:00<00:00, 1947.53 it/sec, obj=[58.33058905 10.70125973]]
        INFO - 16:55:31: ...  58%|█████▊    | 58/100 [00:00<00:00, 1945.44 it/sec, obj=[19.48697889 32.41188233]]
        INFO - 16:55:31: ...  59%|█████▉    | 59/100 [00:00<00:00, 1946.54 it/sec, obj=[ 9.51315338 32.51410675]]
        INFO - 16:55:31: ...  60%|██████    | 60/100 [00:00<00:00, 1943.65 it/sec, obj=[14.29503909 27.34400969]]
        INFO - 16:55:31: ...  61%|██████    | 61/100 [00:00<00:00, 1940.23 it/sec, obj=[21.01471588 31.54351738]]
        INFO - 16:55:31: ...  62%|██████▏   | 62/100 [00:00<00:00, 1936.80 it/sec, obj=[32.62517737 26.61449355]]
        INFO - 16:55:31: ...  63%|██████▎   | 63/100 [00:00<00:00, 1938.32 it/sec, obj=[54.44405197 15.56114393]]
        INFO - 16:55:31: ...  64%|██████▍   | 64/100 [00:00<00:00, 1940.18 it/sec, obj=[31.67860311 19.0346702 ]]
        INFO - 16:55:31: ...  65%|██████▌   | 65/100 [00:00<00:00, 1941.96 it/sec, obj=[12.95434975 27.79409607]]
        INFO - 16:55:31: ...  66%|██████▌   | 66/100 [00:00<00:00, 1943.65 it/sec, obj=[33.53586239 22.39126054]]
        INFO - 16:55:31: ...  67%|██████▋   | 67/100 [00:00<00:00, 1945.71 it/sec, obj=[13.77290554 32.89122016]]
        INFO - 16:55:31: ...  68%|██████▊   | 68/100 [00:00<00:00, 1942.38 it/sec, obj=[24.55394276 22.04978627]]
        INFO - 16:55:31: ...  69%|██████▉   | 69/100 [00:00<00:00, 1939.76 it/sec, obj=[38.26118694 15.93322028]]
        INFO - 16:55:31: ...  70%|███████   | 70/100 [00:00<00:00, 1937.96 it/sec, obj=[69.87762034 24.69083953]]
        INFO - 16:55:31: ...  71%|███████   | 71/100 [00:00<00:00, 1939.37 it/sec, obj=[117.708931     4.51806144]]
        INFO - 16:55:31: ...  72%|███████▏  | 72/100 [00:00<00:00, 1941.12 it/sec, obj=[24.18294227 29.58378487]]
        INFO - 16:55:31: ...  73%|███████▎  | 73/100 [00:00<00:00, 1942.79 it/sec, obj=[65.41388231 10.29184587]]
        INFO - 16:55:31: ...  74%|███████▍  | 74/100 [00:00<00:00, 1944.50 it/sec, obj=[12.24819086 35.25245724]]
        INFO - 16:55:31: ...  75%|███████▌  | 75/100 [00:00<00:00, 1945.75 it/sec, obj=[25.31003032 24.72346394]]
        INFO - 16:55:31: ...  76%|███████▌  | 76/100 [00:00<00:00, 1942.29 it/sec, obj=[49.90530196 13.18421478]]
        INFO - 16:55:31: ...  77%|███████▋  | 77/100 [00:00<00:00, 1940.09 it/sec, obj=[13.58752584 29.15968355]]
        INFO - 16:55:31: ...  78%|███████▊  | 78/100 [00:00<00:00, 1939.25 it/sec, obj=[ 7.8547085  33.00337593]]
        INFO - 16:55:31: ...  79%|███████▉  | 79/100 [00:00<00:00, 1940.52 it/sec, obj=[34.10323875 25.61801713]]
        INFO - 16:55:31: ...  80%|████████  | 80/100 [00:00<00:00, 1942.10 it/sec, obj=[50.66487602 26.86160736]]
        INFO - 16:55:31: ...  81%|████████  | 81/100 [00:00<00:00, 1943.74 it/sec, obj=[19.81349957 25.42966204]]
        INFO - 16:55:31: ...  82%|████████▏ | 82/100 [00:00<00:00, 1945.43 it/sec, obj=[ 7.4546377  34.97394601]]
        INFO - 16:55:31: ...  83%|████████▎ | 83/100 [00:00<00:00, 1944.79 it/sec, obj=[106.96361831   5.44850984]]
        INFO - 16:55:31: ...  84%|████████▍ | 84/100 [00:00<00:00, 1942.50 it/sec, obj=[ 5.36344422 38.7517519 ]]
        INFO - 16:55:31: ...  85%|████████▌ | 85/100 [00:00<00:00, 1940.15 it/sec, obj=[49.116045   14.30366519]]
        INFO - 16:55:31: ...  86%|████████▌ | 86/100 [00:00<00:00, 1940.40 it/sec, obj=[65.46370949  9.53265669]]
        INFO - 16:55:31: ...  87%|████████▋ | 87/100 [00:00<00:00, 1941.95 it/sec, obj=[34.43204876 18.46542147]]
        INFO - 16:55:31: ...  88%|████████▊ | 88/100 [00:00<00:00, 1943.38 it/sec, obj=[89.05896719 23.04880076]]
        INFO - 16:55:31: ...  89%|████████▉ | 89/100 [00:00<00:00, 1944.71 it/sec, obj=[28.30965373 24.49805668]]
        INFO - 16:55:31: ...  90%|█████████ | 90/100 [00:00<00:00, 1945.94 it/sec, obj=[39.10277528 23.64680578]]
        INFO - 16:55:31: ...  91%|█████████ | 91/100 [00:00<00:00, 1944.42 it/sec, obj=[117.48545685   8.18526365]]
        INFO - 16:55:31: ...  92%|█████████▏| 92/100 [00:00<00:00, 1942.29 it/sec, obj=[66.02068516 15.50045449]]
        INFO - 16:55:31: ...  93%|█████████▎| 93/100 [00:00<00:00, 1940.03 it/sec, obj=[67.39952092 12.92818287]]
        INFO - 16:55:31: ...  94%|█████████▍| 94/100 [00:00<00:00, 1941.05 it/sec, obj=[55.28167693 13.5771757 ]]
        INFO - 16:55:31: ...  95%|█████████▌| 95/100 [00:00<00:00, 1942.30 it/sec, obj=[17.335399   26.84155459]]
        INFO - 16:55:31: ...  96%|█████████▌| 96/100 [00:00<00:00, 1943.79 it/sec, obj=[ 4.93473731 35.93049995]]
        INFO - 16:55:31: ...  97%|█████████▋| 97/100 [00:00<00:00, 1945.21 it/sec, obj=[25.16386975 26.89754332]]
        INFO - 16:55:31: ...  98%|█████████▊| 98/100 [00:00<00:00, 1946.61 it/sec, obj=[81.61206359 10.43246637]]
        INFO - 16:55:31: ...  99%|█████████▉| 99/100 [00:00<00:00, 1944.71 it/sec, obj=[40.83406128 15.59204339]]
        INFO - 16:55:31: ... 100%|██████████| 100/100 [00:00<00:00, 1942.77 it/sec, obj=[52.02690192 12.02601654]]
        INFO - 16:55:31: Optimization result:
        INFO - 16:55:31:    Optimizer info:
        INFO - 16:55:31:       Status: None
        INFO - 16:55:31:       Message: None
        INFO - 16:55:31:       Number of calls to the objective function by the optimizer: 100
        INFO - 16:55:31:    Solution:
        INFO - 16:55:31:       The solution is feasible.
        INFO - 16:55:31:       Objective: 30.39964825035985
        INFO - 16:55:31:       Standardized constraints:
        INFO - 16:55:31:          ineq1 = [-11.77907222]
        INFO - 16:55:31:          ineq2 = [-38.26307397]
        INFO - 16:55:31:       Design space: 
        INFO - 16:55:31:       +------+-------------+-------------------+-------------+-------+
        INFO - 16:55:31:       | name | lower_bound |       value       | upper_bound | type  |
        INFO - 16:55:31:       +------+-------------+-------------------+-------------+-------+
        INFO - 16:55:31:       | x    |      0      | 1.542975634014225 |      5      | float |
        INFO - 16:55:31:       | y    |      0      | 1.269910308585573 |      3      | float |
        INFO - 16:55:31:       +------+-------------+-------------------+-------------+-------+

    Optimization result:
       Design variables: [1.54297563 1.26991031]
       Objective function: 30.39964825035985
       Feasible solution: True



.. GENERATED FROM PYTHON SOURCE LINES 60-68

Post-process scenario
---------------------
Lastly, we post-process the scenario by means of the :class:`.ParetoFront`
plot which generates a plot or a matrix of plots if there are more than
2 objectives, plots in blue the locally non dominated points for the current
two objectives, plots in green the globally (all objectives) Pareto optimal
points. The plots in green denotes non-feasible points. Note that the user
can avoid the display of the non-feasible points.

.. GENERATED FROM PYTHON SOURCE LINES 68-88

.. code-block:: default


    PostFactory().execute(
        problem,
        "ParetoFront",
        show_non_feasible=False,
        objectives=["compute_binhkorn"],
        objectives_labels=["f1", "f2"],
        save=False,
        show=True,
    )

    PostFactory().execute(
        problem,
        "ParetoFront",
        show_non_feasible=True,
        objectives=["compute_binhkorn"],
        objectives_labels=["f1", "f2"],
        save=False,
        show=True,
    )



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


    *

      .. image-sg:: /examples/post_process/algorithms/images/sphx_glr_plot_pareto_front_binhkorn_001.png
         :alt: Pareto front
         :srcset: /examples/post_process/algorithms/images/sphx_glr_plot_pareto_front_binhkorn_001.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/post_process/algorithms/images/sphx_glr_plot_pareto_front_binhkorn_002.png
         :alt: Pareto front
         :srcset: /examples/post_process/algorithms/images/sphx_glr_plot_pareto_front_binhkorn_002.png
         :class: sphx-glr-multi-img


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

 .. code-block:: none


    <gemseo.post.pareto_front.ParetoFront object at 0x7fbc55732f40>




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

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


.. _sphx_glr_download_examples_post_process_algorithms_plot_pareto_front_binhkorn.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_pareto_front_binhkorn.py <plot_pareto_front_binhkorn.py>`

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

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


.. only:: html

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

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