MPC 2024 - Issue 2
Refine
Year of publication
- 2024 (4)
Document Type
- Article (4)
Language
- English (4)
Has Fulltext
- yes (4)
This paper describes the computational challenge developed for a computational com-
petition held in 2023 for the 20th anniversary of the Mixed Integer Programming
Workshop. The topic of this competition was reoptimization, also known as warm
starting, of mixed integer linear optimization problems after slight changes to the
input data for a common formulation. The challenge was to accelerate the proof of
optimality of the modified instances by leveraging the information from the solving
processes of previously solved instances, all while creating high-quality primal solu-
tions. Specifically, we discuss the competition’s format, the creation of public and
hidden datasets, and the evaluation criteria. Our goal is to establish a methodology
for the generation of benchmark instances and an evaluation framework, along with
benchmark datasets, to foster future research on reoptimization of mixed integer linear
optimization problems.
The method of nonlinear conjugate gradients (NCG) is widely used in practice for
unconstrained optimization, but it satisfies weak complexity bounds at best when
applied to smooth convex functions. In contrast, Nesterov’s accelerated gradient
(AG) method is optimal up to constant factors for this class. However, when spe-
cialized to quadratic function, conjugate gradient is optimal in a strong sense among
function-gradient methods. Therefore, there is seemingly a gap in the menu of avail-
able algorithms: NCG, the optimal algorithm for quadratic functions that also exhibits
good practical performance for general functions, has poor complexity bounds com-
pared to AG. We propose an NCG method called C+AG (“conjugate plus accelerated
gradient”) to close this gap, that is, it is optimal for quadratic functions and still satisfies
the best possible complexity bound for more general smooth convex functions. It takes
conjugate gradient steps until insufficient progress is made, at which time it switches to
accelerated gradient steps, and later retries conjugate gradient. The proposed method
has the following theoretical properties: (i) It is identical to linear conjugate gradient
(and hence terminates finitely) if the objective function is quadratic; (ii) Its running-
time bound is O(e−1/2) gradient evaluations for an L-smooth convex function, where
is the desired residual reduction, (iii) Its running-time bound is O(√L/l ln(1/e)) if
the function is both L-smooth and -strongly convex. We also conjecture and outline
a proof that a variant of the method has the property: (iv) It is n-step quadratically
convergent for a function whose second derivative is smooth and invertible at the
optimizer. Note that the bounds in (ii) and (iii) match AG and are the best possible,
i.e., they match lower bounds up to constant factors for the classes of functions under
consideration. On the other hand, (i) and (iv) match NCG. In computational tests,
the function-gradient evaluation count for the C+AG method typically behaves as
whichever is better of AG or classical NCG. In some test cases it outperforms both.
AbstractIn statistics, log-concave density estimation is a central problem within the field of nonparametric inference under shape constraints. Despite great progress in recent years on the statistical theory of the canonical estimator, namely the log-concave maximum likelihood estimator, adoption of this method has been hampered by the complexities of the non-smooth convex optimization problem that underpins its computation. We provide enhanced understanding of the structural properties of this optimization problem, which motivates the proposal of new algorithms, based on both randomized and Nesterov smoothing, combined with an appropriate integral discretization of increasing accuracy. We prove that these methods enjoy, both with high probability and in expectation, a convergence rate of order 1/T up to logarithmic factors on the objective function scale, where T denotes the number of iterations. The benefits of our new computational framework are demonstrated on both synthetic and real data, and our implementation is available in a github repository (Log-Concave Computation).