{"id":68,"date":"2020-01-16T15:39:28","date_gmt":"2020-01-16T05:09:28","guid":{"rendered":"http:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/?page_id=68"},"modified":"2020-01-16T15:40:16","modified_gmt":"2020-01-16T05:10:16","slug":"determinant-interior-point-algorithm","status":"publish","type":"page","link":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/hcp-approaches\/determinant-interior-point-algorithm\/","title":{"rendered":"Determinant interior point algorithm"},"content":{"rendered":"<div class=\"wpb-content-wrapper\"><p>[vc_row][vc_column][vc_empty_space][vc_column_text]<\/p>\n<h1>Determinant interior point algorithm<\/h1>\n<p>In 2004, Borkar et al [1] presented the mathematical programming problems, (where <b>G<\/b> is a suitably perturbed\u00a0fundamental matrix)<\/p>\n<p align=\"center\">min <b>G<\/b><sub>11<\/sub><br \/>\nover the set of feasible doubly stochastic policies,<\/p>\n<p>and showed that it was equivalent to the Hamiltonian cycle problem. It was later proved that the gap between optimal values of\u00a0the above problem\u00a0for a Hamiltonian graph and for a non-Hamiltonian graph of the same order is strictly positive [2]. This strengthens the potential usefulness of the optimisation model. However, in implementing an efficient algorithm, numerical difficulties arise. The first major challenge is the evaluation of the <b>G<\/b><sub>11<\/sub> functional, as the fundamental matrix is an inverse. The second major challenge is that relevant Hessian matrices, an important ingredient in many optimisation algorithms, are dense due to perturbation and consequently expensive to compute.<\/p>\n<p>Looking for other functionals to determine Hamiltonicity, Ejov et al [3] proved that the determinant of the inverse of the fundamental matrix, too, is an appropriate optimisation functional. Furthermore, Ejov et al proved that the determinant objective function is more robust than the function <b>G<\/b><sub>11<\/sub>, as it is maximised at Hamiltonian cycles with or without the linear symmetric perturbation required by <b>G<\/b><sub>11<\/sub>.<\/p>\n<p>The following result was proved in [3]. Consider <b>P<\/b>, a probability transition matrix for a given graph, containing continuous variables corresponding to the arcs in the graph. Then, if we restrict <b>P<\/b> to being doubly-stochastic &#8211; that is, <b>P<\/b> has row and column sums of 1, and contains only nonnegative values &#8211; the following optimisation\u00a0problem (where <b>I<\/b> is the <i>N <\/i>\u00d7 <i>N<\/i>\u00a0identity matrix and <b>J<\/b> is an <i>N<\/i>\u00a0\u00d7 <i>N<\/i> matrix with unit entries)<\/p>\n<p align=\"center\">det(<b>I <\/b>&#8211; <b>P<\/b> + <sup>1<\/sup>\/<sub>N<\/sub><b>J<\/b>)<\/p>\n<p align=\"center\">s.t.<\/p>\n<p align=\"center\"><b>P<\/b> is doubly-stochastic,<\/p>\n<p align=\"left\">is bounded below by 0 and above by <i>k<\/i>, where <i>k<\/i> is the length of the longest cycle in the graph. For Hamiltonian graphs, <i>k = N<\/i>, and for non-Hamiltonian graphs, <i>k &lt;= N-1<\/i>. Therefore there is a positive gap of at least 1 for this formulation.<\/p>\n<p align=\"left\">Taking advantage of the inherent sparsity in difficult instances of HCP, Filar et al have developed efficient methods of computing derivatives of the objective function [4]. This method is used in an algorithm called\u00a0DIPA (Determinant interior point algorithm) that seeks to solve the above formulation, detailed in Haythorpe\u00a0[5,6].<\/p>\n<h3>PUBLICATIONS<\/h3>\n<p><span class=\"doi\"><span class=\"value\">[1]\u00a0Borkar, V.S., Ejov, V. and Filar, J.A. &#8220;Directed graphs, Hamiltonicity and doubly stochastic matrices.&#8221; <i>Random Structures and Algorithms<\/i>, 25(4):376-395, 2004. <a href=\"http:\/\/onlinelibrary.wiley.com\/doi\/10.1002\/rsa.20034\/abstract\" target=\"_blank\" rel=\"noopener noreferrer\">Wiley Online Library<\/a><\/span><\/span><\/p>\n<p>[2] Borkar, V.S., Ejov, V. and Filar, J.A. &#8220;On the Hamiltonicity gap and doubly stochastic matrices.&#8221; <i>Random Structures and Algorithms<\/i>, 34(4):502-519, 2009. <a href=\"http:\/\/onlinelibrary.wiley.com\/doi\/10.1002\/rsa.20237\/abstract\" target=\"_blank\" rel=\"noopener noreferrer\">Wiley Online Library<\/a><\/p>\n<p>[3] Ejov, V., Filar, J.A., Murray, W. and Nguyen, G.T. &#8220;Determinants and longest cycles of graphs&#8221;. <i>SIAM Journal on Discrete Mathematics<\/i>, 22(3):1215-1225, 2008. <a href=\"http:\/\/epubs.siam.org\/sidma\/resource\/1\/sjdmec\/v22\/i3\/p1215_s1\" target=\"_blank\" rel=\"noopener noreferrer\">SIAM<\/a><\/p>\n<p>[4] Filar, J.A., Haythorpe, M. and Murray, W. &#8220;On the Determinant and its Derivatives of the Rank-one corrected Generator of a Markov Chain on a Graph.&#8221; <i>Journal of Global Optimization<\/i>, to appear. <a href=\"http:\/\/www.springerlink.com\/content\/x351480m6862m533\/\" target=\"_blank\" rel=\"noopener noreferrer\">SpringerLink<\/a><\/p>\n<p>[5] Haythorpe, M. &#8220;Markov Chain based algorithms for the Hamiltonian cycle problem.&#8221; PhD thesis, University of South Australia, 2010. <a href=\"http:\/\/www.stanford.edu\/group\/SOL\/dissertations.html\" target=\"_blank\" rel=\"noopener noreferrer\">Systems Optimization Laboratory<\/a><\/p>\n<p>[6] Haythorpe, M. &#8220;Finding Hamiltonian cycles using an interior point method.&#8221; <i>Australian Mathematics Society Gazette<\/i>, 37:170-179, 2010. <a href=\"https:\/\/www.austms.org.au\/Publ\/Gazette\/2010\/Jul10\/TechPaperHaythorpe.pdf\" target=\"_blank\" rel=\"noopener noreferrer\">PDF available<\/a><\/p>\n<p>[\/vc_column_text][vc_empty_space][\/vc_column][\/vc_row]<\/p>\n<\/div>","protected":false},"excerpt":{"rendered":"<p>[vc_row][vc_column][vc_empty_space][vc_column_text] Determinant interior point algorithm In 2004, Borkar et al [1] presented the mathematical programming problems, (where G is a suitably perturbed\u00a0fundamental matrix) min G11 over the set of feasible doubly stochastic policies, and showed that it was equivalent to the Hamiltonian cycle problem. It was later proved that the gap between optimal values of\u00a0the [&hellip;]<\/p>\n","protected":false},"author":2,"featured_media":0,"parent":63,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-68","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/pages\/68","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/comments?post=68"}],"version-history":[{"count":0,"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/pages\/68\/revisions"}],"up":[{"embeddable":true,"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/pages\/63"}],"wp:attachment":[{"href":"https:\/\/sites.flinders.edu.au\/flinders-hamiltonian-cycle-project\/wp-json\/wp\/v2\/media?parent=68"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}