(DNN)2: Doubly Non-Negative Relaxations for Deep Neural Networks
Authors: Hanna Jiamei Zhang, Alan Papalia, Michael Everett, David M. Rosen
Organizations: Departments of Computer Science, Electrical and Computer Engineering, and Mathematics Northeastern University, Boston, MA. · Department of Naval Architecture and Marine Engineering, University of Michigan, Ann Arbor, MI.
Existing linear program (LP) and semidefinite program (SDP) relaxations for rectified linear unit (ReLU) neural network (NN) verification yield overly-conservative safety guarantees due to significant relaxation gaps. While the completely positive program (CPP) formulation closes this gap, it is NP-hard to solve. Its cheapest tractable relaxation, the doubly non-negative program (DNN), retains critical constraints as an SDP, but one whose size exceeds the reach of interior-point methods at practical scale. While Burer-Monteiro (BM) factorization has been applied to make SDP-based verification scalable, no such result exists for the strictly tighter DNN formulation. A key obstacle is that additional non-negativity constraints in the DNN cause dual multipliers for optimality certification to be non-unique, making standard certification methods inapplicable. We propose a novel eigenvalue maximization procedure that searches the non-unique multiplier space for a valid certificate, i.e. a global optimality guarantee. Experiments demonstrate that our approach (DNN)2 produces bounds consistently tighter than the standard SDP method, often matching the exact solution, and that our certification procedure confirms global optimality when a valid certificate exists. These results are a key step toward providing tight, certifiable, and computationally scalable verification guarantees needed to deploy neural network controllers and perception modules in safety-critical autonomous systems.
Figures & tables
Fig. 1 : (top) A comparison of the tightness of various verification methods against the ground-truth MILP solution across 100 ReLUNNs . While the standard c-SDP [ 9 , 2 ] leaves significant relaxation gaps, our proposed rank-2 BM-factorized DNN, (DNN)2 , consistently produces tighter bounds, frequently matching the exact MILP solution 95% of all verification instances. Our proposed optimality certification procedure (Eq. 10 ) reliably distinguishes between global and local optima, with the latter marked with \APLstar . (bottom) Wall time computational performance comparison between the direct DNN (solved via interior-point methods, MOSEK ) and (DNN)2 (solved as a non-linear program, KNITRO ) with increasing network size parameter k , using the Julia JuMP interface. As expected, while interior-point solve times grow rapidly (cubically) with problem size, our (DNN)2 exhibits significantly more favorable linear scaling.
Fig. 2 : Analysis of how verification accuracy degrades as NNs increase in height (top) and depth (bottom). The canonical c- SDP bounds (blue) become progressively looser as the network grows, making safety verification more difficult. In contrast, our (DNN)2 formulation (orange) remains highly tight to the MILP ground truth across all tested scales, demonstrating its robustness for larger, more complex network architectures. Our proposed certification procedure (Eq. 10 ) successfully identifies instances where (DNN)2 has reached a local rather than global optimum; these points are marked with \APLstar and are excluded from the box plots. The percentage of instances certified as globally optimal is annotated for each configuration.
Fig. 3 : Visualizing why the DNN relaxation is tighter than c- SDP . Let z1=ReLU(x1+x2) and z2=ReLU(x1−x2) . (left) The feasible set of each relaxation projected onto the post-activation variables (z1,z2) for fixed input x1=x2=0.5 . (right) The same projection across all feasible inputs, under the objective maxz1+z2 . Because the DNN formulation enforces additional element-wise non-negativity constraints not present in the SDP , it produces a strictly smaller feasible set (red) than the SDP (blue), leading to tighter outer NN bounds and more accurate safety guarantees.