Sparse Semidefinite Programs with Guaranteed Near-Linear Time Complexity via Dualized Clique Tree Conversion

TL;DR

This paper introduces a dualized clique tree conversion method for large-scale semidefinite programs, achieving near-linear time complexity by exploiting sparsity patterns, enabling efficient solutions for problems like MAXCUT and power flow relaxations.

Contribution

The paper presents a novel dualization approach that guarantees sparse coupling in clique tree conversion, leading to near-linear time algorithms for certain large-scale semidefinite programs.

Findings

Per-iteration cost is linear in problem size $O(n)$

Achieves $oxed{ ext{near-linear}}$ total time complexity $O(n^{1.5} ext{log}(1/ extepsilon))$

Validated on large problems with up to $n=13659$ variables.

Abstract

Clique tree conversion solves large-scale semidefinite programs by splitting an $n\times n$ matrix variable into up to $n$ smaller matrix variables, each representing a principal submatrix of up to $\omega\times\omega$. Its fundamental weakness is the need to introduce overlap constraints that enforce agreement between different matrix variables, because these can result in dense coupling. In this paper, we show that by dualizing the clique tree conversion, the coupling due to the overlap constraints is guaranteed to be sparse over dense blocks, with a block sparsity pattern that coincides with the adjacency matrix of a tree. We consider two classes of semidefinite programs with favorable sparsity patterns that encompass the MAXCUT and MAX $k$-CUT relaxations, the Lovasz Theta problem, and the AC optimal power flow relaxation. Assuming that $\omega\ll n$, we prove that the per-iteration cost of an interior-point method is linear $O(n)$ time and memory, so an $\epsilon$-accurate and $\epsilon$-feasible iterate is obtained after $O(\sqrt{n}\log(1/\epsilon))$ iterations in near-linear $O(n^{1.5}\log(1/\epsilon))$ time. We confirm our theoretical insights with numerical results on semidefinite programs as large as $n=13659$. (Supporting code at https://github.com/ryz-codes/dual_ctc )