Clash of MINLP Relaxations: Piecewise Linear vs. Global Parabolic

TL;DR

This paper compares piecewise linear and global parabolic relaxations for solving MINLPs, introducing an iterative PARA method that outperforms PWL relaxations at tighter tolerances, especially in problems with sine constraints.

Contribution

The paper presents a novel iterative PARA relaxation method and provides a comprehensive computational comparison with PWL relaxations, highlighting advantages at tighter tolerances.

Findings

PARA relaxations outperform PWL at high tolerances.

PARA maintains problem size without increasing variables.

PARA is particularly effective for sine-constrained problems.

Abstract

Solving mixed-integer nonlinear programs (MINLPs) typically relies on constructing relaxations that are easier to tackle than the original problem. Recently, global parabolic (PARA) relaxations were introduced, featuring separable quadratic functions -- paraboloids -- as global under- or overestimators of general nonlinear constraint functions. So far, the paraboloids are all computed at once by solving a mixed-integer linear program (MIP). For small tolerances or wide function domains, the corresponding MIP grows in size and is eventually intractable, reventing a meaningful comparison with established relaxation techniques. We therefore propose a novel iterative method to compute PARA approximations that succeeds on all tolerance-domain combinations where the original one has failed. The computational study is preceded by a thorough theoretical explanation and analysis. Finally, the improved method enables a computational comparison with piecewise linear (PWL) relaxations in terms of runtime on general MINLP instances. The results show that the modern solver SCIP can solve PWL relaxations faster when the tolerance is high, shifting strongly in favor of PARA for tighter tolerances. We attribute the effect to the difference in the corresponding problem size: PWL relaxations introduce binary variables to identify the active linear piece and their number grows with decreasing tolerance. PARA, on the other hand, does not require additional variables such that the dimension is maintained. For problems with at least one (co)sine constraint, the effect significantly amplifies. Thereby, for medium tolerances, PARA relaxations outperform SCIP stand-alone. Applied problems like alternating current optimal power flow (AC-OPF) feature such constraint types, leaving PARA a viable relaxation strategy.