PSD-Relaxed Multiplicative Gate
Implementation & benchmark of arXiv:2608.16836 — Nonnegative Quadratics over a Quadrant with a Bilinear Constraint
Source paper: Nonnegative Quadratics over a Quadrant with a Bilinear Constraint arXiv:2608.16836 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace an unconstrained multiplicative interaction between two nonnegative neural features by a lifted gate whose first and second moments satisfy the paper's semidefinite relaxation for the set F = {(x1,x2): x1,x2 >= 0, x1 x2 <= 1}. Insert the gate into an MLP, attention score, or MoE router to prevent explosive feature products while retaining a tractable convex feasible set.
Formulas
Mathematical statement
The paper studies F = {(x1,x2) in R^2: x1 >= 0, x2 >= 0, x1 x2 <= 1} and represents a lifted point by Y = [[1,x1,x2],[x1,X11,X12],[x2,X12,X22]], where x1 and x2 are first moments and X11, X12, X22 are second moments. Theorem 2.4 introduces shadow variables s1 and s2 for cubic moments x1^2 x2 and x1 x2^2. The required constraints are M1 = [[x2,X12],[X12,s1]] >= 0, M2 = [[x1,X12],[X12,s2]] >= 0, M3 = [[X11,s1,x1-s1],[s1,X12,0],[x1-s1,0,1-X12]] >= 0, and M4 = [[X22,s2,x2-s2],[s2,X12,0],[x2-s2,0,1-X12]] >= 0. Basic Shor-RLT consequences are 0 <= X12 <= 1, x1,x2 >= 0, X11 >= x1^2, and x1 x2 >= X12^2. The transferable property is that nonlinear product constraints are represented by four constant-size PSD constraints.
Implementation notes
Integrate this as a constrained multiplicative feature block rather than constraining every hidden neuron. Given nonnegative feature vectors a and b, produce a = softplus(u)/c and b = softplus(v)/c so typical coordinates are O(1). For each selected channel pair, introduce lifted variables X11, X12, X22, s1, and s2 and define the interaction output g = X12; use h = W1 a + W2 b + W3 g or h = g * value. During training, solve a projection problem min_z ||z-z_raw||^2 over z = (X11,X12,X22,s1,s2) subject to the four PSD constraints, with x1=a and x2=b. For an MVP, use projected gradient: eigendecompose each 2x2 or 3x3 matrix, clip negative eigenvalues to zero, map the corrected matrices back to shared variables, and repeat 3-10 times. A cvxpylayer implementation can provide a reference result. The paper's matrices and constraints are exact; projection iterations and feature scaling are engineering choices. First test a two-layer CIFAR-10 MLP with ordinary multiplication, clipped multiplication, and this gate at equal parameter count. Measure task loss, minimum eigenvalue, gradient norms, activation tails, and accuracy. Success means fewer exploding-gradient events and equal or higher accuracy at the same training budget.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.