Certifiable Optimization at the ICRA Workshop

GTSAM Posts

Author: Frank Dellaert

This post focuses on our chordal-sparsity paper, while last week’s post already covered CMC-Opt and our new legged-robot estimators, including Yetong Zhang’s workshop paper on constraint manifolds with corners. Both papers will be presented at the Frontiers of Optimization for Robotics workshop at ICRA shows how much momentum there is right now around certifiable optimization, convex relaxation, and structure-exploiting solvers for robotics.

Chordal sparsity

Factor graph lifted to a QCQP, converted to a Bayes tree, and decomposed into clique-wise semidefinite variables.
From a factor graph to a lifted QCQP, then through the Bayes tree to a chordally decomposed SDP relaxation.


The chordal-sparsity paper is here:

Exploiting Chordal Sparsity for Globally Optimal Estimation with Factor Graphs, by Avinash Subramanian, Connor Holmes, Timothy D. Barfoot, Frank Dellaert, and Frederike Dümbgen.

The paper combines factor-graphs and structured optimization via the Bayes-tree with the certifiable-estimation viewpoint that Connor, Tim, and Frederike have helped push forward at the University of Toronto. Local solvers such as Gauss-Newton and Levenberg-Marquardt are fast and exploit sparsity beautifully, but they do not promise that the answer is globally optimal. Convex SDP relaxations can give global solutions or certificates, but the naive monolithic SDP is usually too expensive.

Our contribution is to make the relaxation respect the graph structure. Starting from a GTSAM factor graph, we lift the problem to a QCQP, construct the Bayes tree through variable elimination, and use the resulting cliques to build a chordally decomposed SDP. In plain terms: instead of solving one enormous positive semidefinite matrix problem, we solve many smaller clique-wise matrix problems that follow the same sparse structure GTSAM already understands.

Solver time scaling for a 3D pose-graph SLAM ring factor graph comparing monolithic SDP, chordal SDP, and local solvers.
For the 3D ring pose-graph example, the chordal estimator scales much better than the monolithic SDP while retaining global-optimality guarantees when the relaxation is tight.


The graph above shows the benefit of using the cliques of the Bayes tree: the chordal approach is not as fast a a local solver, but it is globally optimal, at a cost that is vastly less than a monolithic SDP solver, which OOMs on problems of moderate size.

The certifiable wave

The workshop has several papers that combine certifiable optimization, convex relaxations, and structure-exploiting solvers, which is fast becoming a theme in both estimation and control. Our chordal-sparsity paper is part of the “certifiable wave”, but here are many other papers touching on the theme:

GTSAM preview

I am exited to announce that soon GTSAM will support certifiable estimation in a big way, based on the recent QP and QCQP support in GTSAM. QCQPs are a natural bridge between factor-graph models and many convex-relaxation pipelines: once a nonconvex estimation problem is lifted into a quadratically constrained quadratic form, semidefinite relaxations and certificates can be handled ina systematic way. Keep watching this space !

Further browsing

Disclosure: AI was used to help draft this post and prepare the figures.