{"title":"小树宽稀疏半inite程序的和弦转换复杂性","authors":"Richard Y. Zhang","doi":"10.1007/s10107-024-02137-5","DOIUrl":null,"url":null,"abstract":"<p>If a sparse semidefinite program (SDP), specified over <span>\\(n\\times n\\)</span> matrices and subject to <i>m</i> linear constraints, has an aggregate sparsity graph <i>G</i> with small treewidth, then chordal conversion will sometimes allow an interior-point method to solve the SDP in just <span>\\(O(m+n)\\)</span> time per-iteration, which is a significant speedup over the <span>\\(\\varOmega (n^{3})\\)</span> time per-iteration for a direct application of the interior-point method. Unfortunately, the speedup is not guaranteed by an <i>O</i>(1) treewidth in <i>G</i> that is independent of <i>m</i> and <i>n</i>, as a diagonal SDP would have treewidth zero but can still necessitate up to <span>\\(\\varOmega (n^{3})\\)</span> time per-iteration. Instead, we construct an extended aggregate sparsity graph <span>\\(\\overline{G}\\supseteq G\\)</span> by forcing each constraint matrix <span>\\(A_{i}\\)</span> to be its own clique in <i>G</i>. We prove that a small treewidth in <span>\\(\\overline{G}\\)</span> does indeed guarantee that chordal conversion will solve the SDP in <span>\\(O(m+n)\\)</span> time per-iteration, to <span>\\(\\epsilon \\)</span>-accuracy in at most <span>\\(O(\\sqrt{m+n}\\log (1/\\epsilon ))\\)</span> iterations. This sufficient condition covers many successful applications of chordal conversion, including the MAX-<i>k</i>-CUT relaxation, the Lovász theta problem, sensor network localization, polynomial optimization, and the AC optimal power flow relaxation, thus allowing theory to match practical experience.</p>","PeriodicalId":2,"journal":{"name":"ACS Applied Bio Materials","volume":null,"pages":null},"PeriodicalIF":4.6000,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Complexity of chordal conversion for sparse semidefinite programs with small treewidth\",\"authors\":\"Richard Y. Zhang\",\"doi\":\"10.1007/s10107-024-02137-5\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<p>If a sparse semidefinite program (SDP), specified over <span>\\\\(n\\\\times n\\\\)</span> matrices and subject to <i>m</i> linear constraints, has an aggregate sparsity graph <i>G</i> with small treewidth, then chordal conversion will sometimes allow an interior-point method to solve the SDP in just <span>\\\\(O(m+n)\\\\)</span> time per-iteration, which is a significant speedup over the <span>\\\\(\\\\varOmega (n^{3})\\\\)</span> time per-iteration for a direct application of the interior-point method. Unfortunately, the speedup is not guaranteed by an <i>O</i>(1) treewidth in <i>G</i> that is independent of <i>m</i> and <i>n</i>, as a diagonal SDP would have treewidth zero but can still necessitate up to <span>\\\\(\\\\varOmega (n^{3})\\\\)</span> time per-iteration. Instead, we construct an extended aggregate sparsity graph <span>\\\\(\\\\overline{G}\\\\supseteq G\\\\)</span> by forcing each constraint matrix <span>\\\\(A_{i}\\\\)</span> to be its own clique in <i>G</i>. We prove that a small treewidth in <span>\\\\(\\\\overline{G}\\\\)</span> does indeed guarantee that chordal conversion will solve the SDP in <span>\\\\(O(m+n)\\\\)</span> time per-iteration, to <span>\\\\(\\\\epsilon \\\\)</span>-accuracy in at most <span>\\\\(O(\\\\sqrt{m+n}\\\\log (1/\\\\epsilon ))\\\\)</span> iterations. This sufficient condition covers many successful applications of chordal conversion, including the MAX-<i>k</i>-CUT relaxation, the Lovász theta problem, sensor network localization, polynomial optimization, and the AC optimal power flow relaxation, thus allowing theory to match practical experience.</p>\",\"PeriodicalId\":2,\"journal\":{\"name\":\"ACS Applied Bio Materials\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":4.6000,\"publicationDate\":\"2024-09-17\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"ACS Applied Bio Materials\",\"FirstCategoryId\":\"100\",\"ListUrlMain\":\"https://doi.org/10.1007/s10107-024-02137-5\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"MATERIALS SCIENCE, BIOMATERIALS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"ACS Applied Bio Materials","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1007/s10107-024-02137-5","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"MATERIALS SCIENCE, BIOMATERIALS","Score":null,"Total":0}
引用次数: 0
摘要
如果一个在 \(n\times n\) 矩阵上指定并受制于 m 个线性约束的稀疏半定式程序(SDP)有一个具有小树宽的总稀疏性图 G、那么和弦转换有时会让内点法在每次迭代时只需要 \(O(m+n)\) 时间就能求解 SDP,这比直接应用内点法每次迭代所需的\(\varOmega (n^{3})\) 时间大大加快了速度。不幸的是,这种加速并不能通过 G 中与 m 和 n 无关的 O(1) 树状宽度来保证,因为对角 SDP 的树状宽度为零,但每次迭代仍然需要多达 \(\varOmega (n^{3})\) 的时间。相反,我们通过强迫每个约束矩阵 \(A_{i}\) 成为 G 中自己的小块,来构建一个扩展的集合稀疏性图 \(\overline{G}\supseteq G\) 。我们证明,(\overline{G}\)中的小树宽确实可以保证弦变换在每次迭代中以(O(m+n)\)时间求解 SDP,最多以(O(\sqrt{m+n}\log (1/\epsilon ))\) 次迭代达到(\epsilon\)精度。这个充分条件涵盖了和弦转换的许多成功应用,包括 MAX-k-CUT 松弛、Lovász theta 问题、传感器网络定位、多项式优化和交流最优功率流松弛,从而使理论与实践经验相匹配。
Complexity of chordal conversion for sparse semidefinite programs with small treewidth
If a sparse semidefinite program (SDP), specified over \(n\times n\) matrices and subject to m linear constraints, has an aggregate sparsity graph G with small treewidth, then chordal conversion will sometimes allow an interior-point method to solve the SDP in just \(O(m+n)\) time per-iteration, which is a significant speedup over the \(\varOmega (n^{3})\) time per-iteration for a direct application of the interior-point method. Unfortunately, the speedup is not guaranteed by an O(1) treewidth in G that is independent of m and n, as a diagonal SDP would have treewidth zero but can still necessitate up to \(\varOmega (n^{3})\) time per-iteration. Instead, we construct an extended aggregate sparsity graph \(\overline{G}\supseteq G\) by forcing each constraint matrix \(A_{i}\) to be its own clique in G. We prove that a small treewidth in \(\overline{G}\) does indeed guarantee that chordal conversion will solve the SDP in \(O(m+n)\) time per-iteration, to \(\epsilon \)-accuracy in at most \(O(\sqrt{m+n}\log (1/\epsilon ))\) iterations. This sufficient condition covers many successful applications of chordal conversion, including the MAX-k-CUT relaxation, the Lovász theta problem, sensor network localization, polynomial optimization, and the AC optimal power flow relaxation, thus allowing theory to match practical experience.