{"id":2663,"date":"2024-05-29T02:10:27","date_gmt":"2024-05-29T02:10:27","guid":{"rendered":"https:\/\/elo-x.eu\/?p=2663"},"modified":"2024-05-29T02:11:36","modified_gmt":"2024-05-29T02:11:36","slug":"algorithms-for-solving-hard-optimization-problems-2","status":"publish","type":"post","link":"https:\/\/elo-x.eu\/?p=2663","title":{"rendered":"Algorithms for solving hard optimization problems"},"content":{"rendered":"\t\t<div data-elementor-type=\"wp-post\" data-elementor-id=\"2663\" class=\"elementor elementor-2663\">\n\t\t\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-90bc3ec elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"90bc3ec\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-f2b556f\" data-id=\"f2b556f\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-86c9a5c elementor-widget elementor-widget-heading\" data-id=\"86c9a5c\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"heading.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t<h2 class=\"elementor-heading-title elementor-size-xl\">Algorithms for solving nonlinear programs<\/h2>\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t<div class=\"elementor-element elementor-element-1085760 my-divider elementor-widget-divider--view-line elementor-widget elementor-widget-divider\" data-id=\"1085760\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"divider.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t<div class=\"elementor-divider\">\n\t\t\t<span class=\"elementor-divider-separator\">\n\t\t\t\t\t\t<\/span>\n\t\t<\/div>\n\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-beeb054 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"beeb054\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-38efd90\" data-id=\"38efd90\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-a52d48f elementor-widget elementor-widget-text-editor\" data-id=\"a52d48f\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p><span style=\"color: #352a87;\"><span style=\"font-size: 24px;\">Lahcen El Bourkhissi,\u00a0<\/span><\/span><span style=\"color: #352a87;\"><span style=\"font-size: 24px;\">Polytechnic University of Bucharest<\/span><\/span><span style=\"color: #352a87;\"><span style=\"font-size: 24px;\"><br \/><\/span><\/span><\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-31bbb35 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"31bbb35\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-d843487\" data-id=\"d843487\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-391b193 elementor-widget elementor-widget-text-editor\" data-id=\"391b193\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p>In my research, I focus on addressing nonlinear programs, which commonly arise in Nonlinear Model Predictive Control and other areas like machine learning. Solving these problems is challenging due to their nonconvex and sometimes nonsmooth nature. My primary goal is to develop efficient methods based on penalty and augmented Lagrangian approaches capable of effectively handling such challenges. Specifically, I aim to devise methods that involve solving only simple subproblems, such as convex subproblems.<\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-6c37975 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"6c37975\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-2c92c57\" data-id=\"2c92c57\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-b7a88c2 my-divider elementor-widget-divider--view-line elementor-widget elementor-widget-divider\" data-id=\"b7a88c2\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"divider.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t<div class=\"elementor-divider\">\n\t\t\t<span class=\"elementor-divider-separator\">\n\t\t\t\t\t\t<\/span>\n\t\t<\/div>\n\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t<div class=\"elementor-element elementor-element-6f16671 elementor-widget elementor-widget-heading\" data-id=\"6f16671\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"heading.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t<h2 class=\"elementor-heading-title elementor-size-default\">Penalty-based method for solving nonlinear programs<\/h2>\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-80ad4c4 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"80ad4c4\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-66cbae4\" data-id=\"66cbae4\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-82d499d elementor-widget elementor-widget-text-editor\" data-id=\"82d499d\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p>[1] El Bourkhissi, Necoara &#8220;Complexity of linearized quadratic penalty for optimization with nonlinear equality constraints&#8221; <em>&#8211; submitted, 2023<\/em><\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-8722eb8 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"8722eb8\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-50 elementor-top-column elementor-element elementor-element-72a2899\" data-id=\"72a2899\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-957f7c1 elementor-widget elementor-widget-text-editor\" data-id=\"957f7c1\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p>\u00a0<\/p><p>IN [1], we considered a nonconvex optimization problem with nonlinear equality constraints. We assumed that both the objective function and the functional constraints were locally smooth. To solve this problem, we applied a linearized quadratic penalty method. That is, we linearized the objective function and the functional constraints in the penalty formulation at the current iterate and added a quadratic regularization, thus yielding a subproblem that is easy to solve (in fact, the solution of our subproblem is equiavlent to the solving a linear system), and whose solution was the next iterate. Under a dynamic regularization parameter choice, we derived convergence guarantees for the iterates of our method to an \u03b5 first-order optimal solution in O(1\/\u03b5^3) evaluations of the first order information of the objective function and functional constraints. Finally, we showed that when the problem data satisfied the Kurdyka-Lojasiewicz property, e.g., were semialgebraic, the whole sequence generated by our algorithm converged and we derived convergence rates. We validated the theory and the performance of the proposed algorithm by numerically comparing it with some existing methods from the literature.<\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t<div class=\"elementor-column elementor-col-50 elementor-top-column elementor-element elementor-element-73c1594\" data-id=\"73c1594\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-705973c elementor-widget__width-initial elementor-widget elementor-widget-image\" data-id=\"705973c\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"image.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t\t\t\t<figure class=\"wp-caption\">\n\t\t\t\t\t\t\t\t\t\t<img fetchpriority=\"high\" decoding=\"async\" width=\"525\" height=\"309\" src=\"https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/perfprof_latest-1-1024x603.png\" class=\"attachment-large size-large wp-image-2664\" alt=\"\" srcset=\"https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/perfprof_latest-1-1024x603.png 1024w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/perfprof_latest-1-300x177.png 300w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/perfprof_latest-1-768x452.png 768w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/perfprof_latest-1.png 1412w\" sizes=\"100vw\" \/>\t\t\t\t\t\t\t\t\t\t\t<figcaption class=\"widget-image-caption wp-caption-text\">Performance profiles for the computation time (left) and the number of iterations (right).<\/figcaption>\n\t\t\t\t\t\t\t\t\t\t<\/figure>\n\t\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-3d736fa elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"3d736fa\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-3aee274\" data-id=\"3aee274\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-0393100 my-divider elementor-widget-divider--view-line elementor-widget elementor-widget-divider\" data-id=\"0393100\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"divider.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t<div class=\"elementor-divider\">\n\t\t\t<span class=\"elementor-divider-separator\">\n\t\t\t\t\t\t<\/span>\n\t\t<\/div>\n\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t<div class=\"elementor-element elementor-element-2d893c3 elementor-widget elementor-widget-heading\" data-id=\"2d893c3\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"heading.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t<h2 class=\"elementor-heading-title elementor-size-default\">ADMM scheme for nonsmooth nonlinear programs<\/h2>\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-3d870f3 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"3d870f3\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-100 elementor-top-column elementor-element elementor-element-2ecb947\" data-id=\"2ecb947\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-d38de93 elementor-widget elementor-widget-text-editor\" data-id=\"d38de93\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p>[2] El Bourkhissi, Necoara, Patrinos, &#8220;Linearized ADMM for nonsmooth nonconvex optimization with nonlinear equality constraints&#8221;, CDC 2023<\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<section class=\"elementor-section elementor-top-section elementor-element elementor-element-fff9796 elementor-section-boxed elementor-section-height-default elementor-section-height-default\" data-id=\"fff9796\" data-element_type=\"section\" data-e-type=\"section\">\n\t\t\t\t\t\t<div class=\"elementor-container elementor-column-gap-default\">\n\t\t\t\t\t<div class=\"elementor-column elementor-col-50 elementor-top-column elementor-element elementor-element-24920a1\" data-id=\"24920a1\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-5a37d73 elementor-widget elementor-widget-text-editor\" data-id=\"5a37d73\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"text-editor.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t<p><span style=\"color: var( --e-global-color-text ); font-weight: var( --e-global-typography-text-font-weight ); font-size: 1rem;\">In [2], we proposed a new approach for solving a structured nonsmooth nonconvex optimization problem with nonlinear equality constraints, where both the objective function and constraints were 2-block separable. Our method was based on a 2-block linearized ADMM, where we linearized the smooth part of the cost function and the nonlinear term of the functional constraints in the augmented Lagrangian at each outer iteration. This resulted in simple subproblems, whose solutions were used to update the iterates of the 2 blocks variables. We proved global convergence for the sequence generated by our method to a stationary point of the original problem. We applied our proposed algorithm as a solver for the nonlinear model predictive control problem of an inverted pendulum on a cart.<\/span><\/p>\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t<div class=\"elementor-column elementor-col-50 elementor-top-column elementor-element elementor-element-643dfc9\" data-id=\"643dfc9\" data-element_type=\"column\" data-e-type=\"column\">\n\t\t\t<div class=\"elementor-widget-wrap elementor-element-populated\">\n\t\t\t\t\t\t<div class=\"elementor-element elementor-element-ac5f642 elementor-widget__width-initial elementor-widget elementor-widget-image\" data-id=\"ac5f642\" data-element_type=\"widget\" data-e-type=\"widget\" data-widget_type=\"image.default\">\n\t\t\t\t<div class=\"elementor-widget-container\">\n\t\t\t\t\t\t\t\t\t\t\t\t<figure class=\"wp-caption\">\n\t\t\t\t\t\t\t\t\t\t<img decoding=\"async\" width=\"525\" height=\"257\" src=\"https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/result_final-1-1024x501.png\" class=\"attachment-large size-large wp-image-2665\" alt=\"\" srcset=\"https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/result_final-1-1024x501.png 1024w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/result_final-1-300x147.png 300w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/result_final-1-768x376.png 768w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/result_final-1-1536x752.png 1536w, https:\/\/elo-x.eu\/wp-content\/uploads\/2024\/05\/result_final-1.png 1657w\" sizes=\"100vw\" \/>\t\t\t\t\t\t\t\t\t\t\t<figcaption class=\"widget-image-caption wp-caption-text\">Nonlinear MPC inputs $u$ (top), cart lateral position $z_1$ (middle) and pendulum angular position $z_3$ (bottom) computed using Linearized ADMM<\/figcaption>\n\t\t\t\t\t\t\t\t\t\t<\/figure>\n\t\t\t\t\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/div>\n\t\t\t\t\t<\/div>\n\t\t<\/section>\n\t\t\t\t<\/div>\n\t\t","protected":false},"excerpt":{"rendered":"<p>Algorithms for solving nonlinear programs Lahcen El Bourkhissi,\u00a0Polytechnic University of Bucharest In my research, I focus on addressing nonlinear programs, which commonly arise in Nonlinear Model Predictive Control and other areas like machine learning. Solving these problems is challenging due to their nonconvex and sometimes nonsmooth nature. My primary goal is to develop efficient methods &hellip; <\/p>\n<p class=\"link-more\"><a href=\"https:\/\/elo-x.eu\/?p=2663\" class=\"more-link\">Continue reading<span class=\"screen-reader-text\"> &#8220;Algorithms for solving hard optimization problems&#8221;<\/span><\/a><\/p>\n","protected":false},"author":19,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"elementor_header_footer","format":"standard","meta":{"footnotes":""},"categories":[22],"tags":[],"class_list":["post-2663","post","type-post","status-publish","format-standard","hentry","category-projects"],"_links":{"self":[{"href":"https:\/\/elo-x.eu\/index.php?rest_route=\/wp\/v2\/posts\/2663","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/elo-x.eu\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/elo-x.eu\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/elo-x.eu\/index.php?rest_route=\/wp\/v2\/users\/19"}],"replies":[{"embeddable":true,"href":"https:\/\/elo-x.eu\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=2663"}],"version-history":[{"count":3,"href":"https:\/\/elo-x.eu\/index.php?rest_route=\/wp\/v2\/posts\/2663\/revisions"}],"predecessor-version":[{"id":2668,"href":"https:\/\/elo-x.eu\/index.php?rest_route=\/wp\/v2\/posts\/2663\/revisions\/2668"}],"wp:attachment":[{"href":"https:\/\/elo-x.eu\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=2663"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/elo-x.eu\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=2663"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/elo-x.eu\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=2663"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}