After preconditioning the defining linear map, global minimization of the dual log-partition function certifies moment body membership, and L-BFGS handles dense n=m=1000 instances in seconds.
Shannon-and von neumann-entropy regularizations of linear and semidefinite programs
1 Pith paper cite this work. Polarity classification is still indexing.
abstract
We consider the LP in standard form min {c T x\,: Ax = b; x $\ge$ 0} and inspired by $\epsilon$-regularization in Optimal Transport, we introduce its $\epsilon$-regularization ''min {c T x + $\epsilon$ f (x)\,: Ax = b; x $\ge$ 0}'' via the (convex) Boltzmann-Shannon entropy f (x)\,:= i x i ln x i . We also provide a similar regularization for the semidefinite program ''min {Tr(C $\bullet$ X)\,: A(X) = b; X 0}'' but with now the so-called Von Neumann entropy, as in Quantum Optimal Transport. Importantly, both are not barriers of the LP and SDP cones respectively. We show that this problem admits an equivalent unconstrained convex problem max $\lambda$$\in$R m G$\epsilon$($\lambda$) for an explicit concave differentiable function G$\epsilon$ in dual variables $\lambda$ $\in$ R m . As $\epsilon$ goes to zero, its optimal value converges to the optimal value of the initial LP. While it resembles the log-barrier formulation of interior point algorithm for the initial LP, it has a distinguishing advantage. Namely for fixed $\lambda$, G$\epsilon$($\lambda$) is obtained as a minimization over the whole space x $\in$ R d (and not over x $\ge$ 0) to still obtain a nonnegative solution x($\lambda$) $\ge$ 0, whence an explicit form of G$\epsilon$ very useful for its unconstrained maximization over R m .
citation-role summary
citation-polarity summary
fields
math.OC 1years
2025 1verdicts
CONDITIONAL 1roles
other 1polarities
unclear 1representative citing papers
citing papers explorer
-
Maximal entropy in the moment body
After preconditioning the defining linear map, global minimization of the dual log-partition function certifies moment body membership, and L-BFGS handles dense n=m=1000 instances in seconds.