Note

Go to the end to download the full example code

Pareto front on Binh and Korn problem

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

from __future__ import annotations

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

Import

The first step is to import a high-level function for logging.

configure_logger()

<RootLogger root (INFO)>

Import the optimization problem

Then, we instantiate the Binh and Korn optimization problem (see BinhKorn).

problem = BinhKorn()

Create and execute scenario

Then, we instantiate the design of experiment factory, and we request the execution of a 100-length LHS optimized by simulated annealing.

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

INFO - 08:27:52: Optimization problem:
INFO - 08:27:52:    minimize compute_binhkorn(x, y) = (4*x**2+ 4*y**2, (x-5.)**2 + (y-5.)**2)
INFO - 08:27:52:    with respect to x, y
INFO - 08:27:52:    subject to constraints:
INFO - 08:27:52:       ineq1(x, y): (x-5.)**2 + y**2 <= 25. <= 0.0
INFO - 08:27:52:       ineq2(x, y): (x-8.)**2 + (y+3)**2 >= 7.7 <= 0.0
INFO - 08:27:52:    over the design space:
INFO - 08:27:52:    +------+-------------+-------+-------------+-------+
INFO - 08:27:52:    | name | lower_bound | value | upper_bound | type  |
INFO - 08:27:52:    +------+-------------+-------+-------------+-------+
INFO - 08:27:52:    | x    |      0      |   1   |      5      | float |
INFO - 08:27:52:    | y    |      0      |   1   |      3      | float |
INFO - 08:27:52:    +------+-------------+-------+-------------+-------+
INFO - 08:27:52: Solving optimization problem with algorithm OT_OPT_LHS:
INFO - 08:27:52: ...   0%|          | 0/100 [00:00<?, ?it]
INFO - 08:27:52: ...   1%|          | 1/100 [00:00<00:00, 2799.94 it/sec, obj=[98.50728236 17.96432773]]
INFO - 08:27:52: ...   2%|▏         | 2/100 [00:00<00:00, 2055.53 it/sec, obj=[ 5.10086961 36.30607322]]
INFO - 08:27:52: ...   3%|▎         | 3/100 [00:00<00:00, 1980.00 it/sec, obj=[ 5.73344945 34.53104096]]
INFO - 08:27:52: ...   4%|▍         | 4/100 [00:00<00:00, 1958.58 it/sec, obj=[ 2.99274814 38.89733089]]
INFO - 08:27:52: ...   5%|▌         | 5/100 [00:00<00:00, 1954.48 it/sec, obj=[ 1.45355254 43.72252794]]
INFO - 08:27:52: ...   6%|▌         | 6/100 [00:00<00:00, 1954.63 it/sec, obj=[26.6936315  22.26390516]]
INFO - 08:27:52: ...   7%|▋         | 7/100 [00:00<00:00, 1956.17 it/sec, obj=[95.06326723 19.21865528]]
INFO - 08:27:52: ...   8%|▊         | 8/100 [00:00<00:00, 1956.41 it/sec, obj=[11.37121324 33.49445422]]
INFO - 08:27:52: ...   9%|▉         | 9/100 [00:00<00:00, 1927.23 it/sec, obj=[33.52756721 17.72645069]]
INFO - 08:27:52: ...  10%|█         | 10/100 [00:00<00:00, 1927.35 it/sec, obj=[ 4.39737722 39.0292696 ]]
INFO - 08:27:52: ...  11%|█         | 11/100 [00:00<00:00, 1922.87 it/sec, obj=[83.81080473  6.44933463]]
INFO - 08:27:52: ...  12%|█▏        | 12/100 [00:00<00:00, 1908.53 it/sec, obj=[64.12902701 21.25106868]]
INFO - 08:27:52: ...  13%|█▎        | 13/100 [00:00<00:00, 1906.17 it/sec, obj=[71.23914138  9.90898899]]
INFO - 08:27:52: ...  14%|█▍        | 14/100 [00:00<00:00, 1908.05 it/sec, obj=[58.61362754 10.51836124]]
INFO - 08:27:52: ...  15%|█▌        | 15/100 [00:00<00:00, 1909.74 it/sec, obj=[17.1313131  28.27257144]]
INFO - 08:27:52: ...  16%|█▌        | 16/100 [00:00<00:00, 1838.30 it/sec, obj=[38.07815452 16.26141839]]
INFO - 08:27:52: ...  17%|█▋        | 17/100 [00:00<00:00, 1836.62 it/sec, obj=[45.29614269 15.87407947]]
INFO - 08:27:52: ...  18%|█▊        | 18/100 [00:00<00:00, 1837.37 it/sec, obj=[78.26474617 13.97242459]]
INFO - 08:27:52: ...  19%|█▉        | 19/100 [00:00<00:00, 1842.20 it/sec, obj=[58.47852866 16.87161502]]
INFO - 08:27:52: ...  20%|██        | 20/100 [00:00<00:00, 1847.39 it/sec, obj=[ 7.27923171 33.31157275]]
INFO - 08:27:52: ...  21%|██        | 21/100 [00:00<00:00, 1852.76 it/sec, obj=[15.973784   25.86458657]]
INFO - 08:27:52: ...  22%|██▏       | 22/100 [00:00<00:00, 1857.57 it/sec, obj=[47.66727367 15.96717933]]
INFO - 08:27:52: ...  23%|██▎       | 23/100 [00:00<00:00, 1861.80 it/sec, obj=[88.93580164  6.7712213 ]]
INFO - 08:27:52: ...  24%|██▍       | 24/100 [00:00<00:00, 1836.05 it/sec, obj=[22.10070032 24.30208587]]
INFO - 08:27:52: ...  25%|██▌       | 25/100 [00:00<00:00, 1836.06 it/sec, obj=[16.24176378 31.32352362]]
INFO - 08:27:52: ...  26%|██▌       | 26/100 [00:00<00:00, 1838.96 it/sec, obj=[29.38105263 25.45614728]]
INFO - 08:27:52: ...  27%|██▋       | 27/100 [00:00<00:00, 1838.89 it/sec, obj=[27.91051427 19.91664614]]
INFO - 08:27:52: ...  28%|██▊       | 28/100 [00:00<00:00, 1839.26 it/sec, obj=[37.28915469 18.0103703 ]]
INFO - 08:27:52: ...  29%|██▉       | 29/100 [00:00<00:00, 1836.77 it/sec, obj=[ 1.01177424 43.45819442]]
INFO - 08:27:52: ...  30%|███       | 30/100 [00:00<00:00, 1839.04 it/sec, obj=[85.96675181 10.67473104]]
INFO - 08:27:52: ...  31%|███       | 31/100 [00:00<00:00, 1838.72 it/sec, obj=[87.34794862  6.40751572]]
INFO - 08:27:52: ...  32%|███▏      | 32/100 [00:00<00:00, 1841.60 it/sec, obj=[25.57943091 20.74700556]]
INFO - 08:27:52: ...  33%|███▎      | 33/100 [00:00<00:00, 1842.40 it/sec, obj=[76.39229092 11.95231974]]
INFO - 08:27:52: ...  34%|███▍      | 34/100 [00:00<00:00, 1844.87 it/sec, obj=[51.64657249 21.13171808]]
INFO - 08:27:52: ...  35%|███▌      | 35/100 [00:00<00:00, 1847.71 it/sec, obj=[95.8196038   9.31540942]]
INFO - 08:27:52: ...  36%|███▌      | 36/100 [00:00<00:00, 1850.59 it/sec, obj=[117.01496061   6.96792478]]
INFO - 08:27:52: ...  37%|███▋      | 37/100 [00:00<00:00, 1853.58 it/sec, obj=[78.37371654 21.95951365]]
INFO - 08:27:52: ...  38%|███▊      | 38/100 [00:00<00:00, 1854.07 it/sec, obj=[38.21112297 20.46923673]]
INFO - 08:27:52: ...  39%|███▉      | 39/100 [00:00<00:00, 1855.97 it/sec, obj=[77.87964266 17.24338363]]
INFO - 08:27:52: ...  40%|████      | 40/100 [00:00<00:00, 1857.14 it/sec, obj=[59.79323468 22.48208172]]
INFO - 08:27:52: ...  41%|████      | 41/100 [00:00<00:00, 1858.86 it/sec, obj=[34.56674439 21.30358144]]
INFO - 08:27:52: ...  42%|████▏     | 42/100 [00:00<00:00, 1861.18 it/sec, obj=[ 0.46290462 45.68910522]]
INFO - 08:27:52: ...  43%|████▎     | 43/100 [00:00<00:00, 1862.90 it/sec, obj=[100.73843409  11.99829359]]
INFO - 08:27:52: ...  44%|████▍     | 44/100 [00:00<00:00, 1864.96 it/sec, obj=[21.18253098 23.10192252]]
INFO - 08:27:52: ...  45%|████▌     | 45/100 [00:00<00:00, 1867.21 it/sec, obj=[51.0570421  12.66142788]]
INFO - 08:27:52: ...  46%|████▌     | 46/100 [00:00<00:00, 1866.89 it/sec, obj=[90.5172837  14.59877358]]
INFO - 08:27:52: ...  47%|████▋     | 47/100 [00:00<00:00, 1868.43 it/sec, obj=[19.5728568  23.86163594]]
INFO - 08:27:52: ...  48%|████▊     | 48/100 [00:00<00:00, 1866.76 it/sec, obj=[47.52592344 19.80950276]]
INFO - 08:27:52: ...  49%|████▉     | 49/100 [00:00<00:00, 1867.95 it/sec, obj=[113.47621648   6.0816867 ]]
INFO - 08:27:52: ...  50%|█████     | 50/100 [00:00<00:00, 1869.40 it/sec, obj=[37.3534796  16.14867587]]
INFO - 08:27:52: ...  51%|█████     | 51/100 [00:00<00:00, 1870.79 it/sec, obj=[41.42991507 25.64561402]]
INFO - 08:27:52: ...  52%|█████▏    | 52/100 [00:00<00:00, 1872.54 it/sec, obj=[11.19222105 30.06624599]]
INFO - 08:27:52: ...  53%|█████▎    | 53/100 [00:00<00:00, 1874.24 it/sec, obj=[86.67968625 15.70216699]]
INFO - 08:27:52: ...  54%|█████▍    | 54/100 [00:00<00:00, 1873.80 it/sec, obj=[90.96123503  7.41344494]]
INFO - 08:27:52: ...  55%|█████▌    | 55/100 [00:00<00:00, 1875.33 it/sec, obj=[34.69157131 20.23247304]]
INFO - 08:27:52: ...  56%|█████▌    | 56/100 [00:00<00:00, 1875.61 it/sec, obj=[56.73716821 11.35918483]]
INFO - 08:27:52: ...  57%|█████▋    | 57/100 [00:00<00:00, 1876.94 it/sec, obj=[58.33058905 10.70125973]]
INFO - 08:27:52: ...  58%|█████▊    | 58/100 [00:00<00:00, 1878.38 it/sec, obj=[19.48697889 32.41188233]]
INFO - 08:27:52: ...  59%|█████▉    | 59/100 [00:00<00:00, 1879.57 it/sec, obj=[ 9.51315338 32.51410675]]
INFO - 08:27:52: ...  60%|██████    | 60/100 [00:00<00:00, 1881.11 it/sec, obj=[14.29503909 27.34400969]]
INFO - 08:27:52: ...  61%|██████    | 61/100 [00:00<00:00, 1882.36 it/sec, obj=[21.01471588 31.54351738]]
INFO - 08:27:52: ...  62%|██████▏   | 62/100 [00:00<00:00, 1881.98 it/sec, obj=[32.62517737 26.61449355]]
INFO - 08:27:52: ...  63%|██████▎   | 63/100 [00:00<00:00, 1883.36 it/sec, obj=[54.44405197 15.56114393]]
INFO - 08:27:52: ...  64%|██████▍   | 64/100 [00:00<00:00, 1883.57 it/sec, obj=[31.67860311 19.0346702 ]]
INFO - 08:27:52: ...  65%|██████▌   | 65/100 [00:00<00:00, 1884.92 it/sec, obj=[12.95434975 27.79409607]]
INFO - 08:27:52: ...  66%|██████▌   | 66/100 [00:00<00:00, 1886.07 it/sec, obj=[33.53586239 22.39126054]]
INFO - 08:27:52: ...  67%|██████▋   | 67/100 [00:00<00:00, 1887.40 it/sec, obj=[13.77290554 32.89122016]]
INFO - 08:27:52: ...  68%|██████▊   | 68/100 [00:00<00:00, 1888.61 it/sec, obj=[24.55394276 22.04978627]]
INFO - 08:27:52: ...  69%|██████▉   | 69/100 [00:00<00:00, 1888.44 it/sec, obj=[38.26118694 15.93322028]]
INFO - 08:27:52: ...  70%|███████   | 70/100 [00:00<00:00, 1889.22 it/sec, obj=[69.87762034 24.69083953]]
INFO - 08:27:52: ...  71%|███████   | 71/100 [00:00<00:00, 1889.24 it/sec, obj=[117.708931     4.51806144]]
INFO - 08:27:52: ...  72%|███████▏  | 72/100 [00:00<00:00, 1889.93 it/sec, obj=[24.18294227 29.58378487]]
INFO - 08:27:52: ...  73%|███████▎  | 73/100 [00:00<00:00, 1890.97 it/sec, obj=[65.41388231 10.29184587]]
INFO - 08:27:52: ...  74%|███████▍  | 74/100 [00:00<00:00, 1891.89 it/sec, obj=[12.24819086 35.25245724]]
INFO - 08:27:52: ...  75%|███████▌  | 75/100 [00:00<00:00, 1892.92 it/sec, obj=[25.31003032 24.72346394]]
INFO - 08:27:52: ...  76%|███████▌  | 76/100 [00:00<00:00, 1893.75 it/sec, obj=[49.90530196 13.18421478]]
INFO - 08:27:52: ...  77%|███████▋  | 77/100 [00:00<00:00, 1893.40 it/sec, obj=[13.58752584 29.15968355]]
INFO - 08:27:52: ...  78%|███████▊  | 78/100 [00:00<00:00, 1894.00 it/sec, obj=[ 7.8547085  33.00337593]]
INFO - 08:27:52: ...  79%|███████▉  | 79/100 [00:00<00:00, 1893.76 it/sec, obj=[34.10323875 25.61801713]]
INFO - 08:27:52: ...  80%|████████  | 80/100 [00:00<00:00, 1894.69 it/sec, obj=[50.66487602 26.86160736]]
INFO - 08:27:52: ...  81%|████████  | 81/100 [00:00<00:00, 1895.73 it/sec, obj=[19.81349957 25.42966204]]
INFO - 08:27:52: ...  82%|████████▏ | 82/100 [00:00<00:00, 1896.72 it/sec, obj=[ 7.4546377  34.97394601]]
INFO - 08:27:52: ...  83%|████████▎ | 83/100 [00:00<00:00, 1897.65 it/sec, obj=[106.96361831   5.44850984]]
INFO - 08:27:52: ...  84%|████████▍ | 84/100 [00:00<00:00, 1898.67 it/sec, obj=[ 5.36344422 38.7517519 ]]
INFO - 08:27:52: ...  85%|████████▌ | 85/100 [00:00<00:00, 1897.83 it/sec, obj=[49.116045   14.30366519]]
INFO - 08:27:52: ...  86%|████████▌ | 86/100 [00:00<00:00, 1898.38 it/sec, obj=[65.46370949  9.53265669]]
INFO - 08:27:52: ...  87%|████████▋ | 87/100 [00:00<00:00, 1898.31 it/sec, obj=[34.43204876 18.46542147]]
INFO - 08:27:52: ...  88%|████████▊ | 88/100 [00:00<00:00, 1896.27 it/sec, obj=[89.05896719 23.04880076]]
INFO - 08:27:52: ...  89%|████████▉ | 89/100 [00:00<00:00, 1896.28 it/sec, obj=[28.30965373 24.49805668]]
INFO - 08:27:52: ...  90%|█████████ | 90/100 [00:00<00:00, 1896.85 it/sec, obj=[39.10277528 23.64680578]]
INFO - 08:27:52: ...  91%|█████████ | 91/100 [00:00<00:00, 1897.35 it/sec, obj=[117.48545685   8.18526365]]
INFO - 08:27:52: ...  92%|█████████▏| 92/100 [00:00<00:00, 1896.88 it/sec, obj=[66.02068516 15.50045449]]
INFO - 08:27:52: ...  93%|█████████▎| 93/100 [00:00<00:00, 1897.02 it/sec, obj=[67.39952092 12.92818287]]
INFO - 08:27:52: ...  94%|█████████▍| 94/100 [00:00<00:00, 1896.92 it/sec, obj=[55.28167693 13.5771757 ]]
INFO - 08:27:52: ...  95%|█████████▌| 95/100 [00:00<00:00, 1897.22 it/sec, obj=[17.335399   26.84155459]]
INFO - 08:27:52: ...  96%|█████████▌| 96/100 [00:00<00:00, 1897.91 it/sec, obj=[ 4.93473731 35.93049995]]
INFO - 08:27:52: ...  97%|█████████▋| 97/100 [00:00<00:00, 1898.71 it/sec, obj=[25.16386975 26.89754332]]
INFO - 08:27:52: ...  98%|█████████▊| 98/100 [00:00<00:00, 1899.47 it/sec, obj=[81.61206359 10.43246637]]
INFO - 08:27:52: ...  99%|█████████▉| 99/100 [00:00<00:00, 1899.99 it/sec, obj=[40.83406128 15.59204339]]
INFO - 08:27:52: ... 100%|██████████| 100/100 [00:00<00:00, 1899.58 it/sec, obj=[52.02690192 12.02601654]]
INFO - 08:27:52: Optimization result:
INFO - 08:27:52:    Optimizer info:
INFO - 08:27:52:       Status: None
INFO - 08:27:52:       Message: None
INFO - 08:27:52:       Number of calls to the objective function by the optimizer: 100
INFO - 08:27:52:    Solution:
INFO - 08:27:52:       The solution is feasible.
INFO - 08:27:52:       Objective: 30.39964825035985
INFO - 08:27:52:       Standardized constraints:
INFO - 08:27:52:          ineq1 = [-11.77907222]
INFO - 08:27:52:          ineq2 = [-38.26307397]
INFO - 08:27:52:       Design space:
INFO - 08:27:52:       +------+-------------+-------------------+-------------+-------+
INFO - 08:27:52:       | name | lower_bound |       value       | upper_bound | type  |
INFO - 08:27:52:       +------+-------------+-------------------+-------------+-------+
INFO - 08:27:52:       | x    |      0      | 1.542975634014225 |      5      | float |
INFO - 08:27:52:       | y    |      0      | 1.269910308585573 |      3      | float |
INFO - 08:27:52:       +------+-------------+-------------------+-------------+-------+
Optimization result:
  • Design variables: [1.54297563 1.26991031]
  • Objective function: 30.39964825035985
  • Feasible solution: True


Post-process scenario

Lastly, we post-process the scenario by means of the 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 denote non-feasible points. Note that the user can avoid the display of the non-feasible points.

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

PostFactory().execute(
    problem,
    "ParetoFront",
    objectives=["compute_binhkorn"],
    objectives_labels=["f1", "f2"],
    save=False,
    show=True,
)
  • Pareto front
  • Pareto front
<gemseo.post.pareto_front.ParetoFront object at 0x7f0cb13cc340>

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

Download Python source code: plot_pareto_front_binhkorn.py

Download Jupyter notebook: plot_pareto_front_binhkorn.ipynb

Gallery generated by Sphinx-Gallery