The group has a strong experience in research and development of computational techniques. This is reflected in the production and implementation of computational software.

24 parallel sparse PCA codes based on alternating maximization by Peter Richtárik and Martin Takáč.

Accelerated coordinate descent methods for minimizing composite functions by Peter Richtárik and Martin Takáč.

*EMSOL* is an implementation of the simplex method written and
maintained by Julian Hall.

*HighOrderReverse* (High Order Reverse Automatic Differentiation)
is a C++ implementation of reverse Automatic Differentiation routines that
calculate the Hessian matrix, its sparsity pattern and the directional
derivative of the Hessian matrix. It is an additional driver of ADOL-C package
and was written by Robert Gower.

More information about HighOrderReverse

*HOPDM* (Higher Order Primal-Dual Method) is Jacek Gondzio's
implementation of an infeasible primal-dual path-following interior point
method for linear, convex quadratic and convex nonlinear programming
problems.

HOPDM allows solving large scale linear, convex quadratic and convex nonlinear programming problems. The code algorithm uses multiple centrality correctors; their number is chosen appropriately for a given problem in order to reduce the overall solution time. HOPDM automatically chooses the most efficient factorization method for a given problem (either normal equations or augmented system). The code compares favourably with commercial LP, QP and NLP packages.

*MFIPMCS* (Matrix-free Interior Point Method for Compressed Sensing)
is an interior point method implemented in MATLAB for the solution of real
valued compressed sensing problems.

The "matrix-free" implies that only matrix-vector products operations are allowed and the process is memoryless. The solver employs an efficient preconditioning technique along with a Krylov subspace method (i.e conjugate gradient) for the fast solution of linear systems at every iteration.

More information about MFIPMCS

*OOPS* (Object-Oriented Parallel Solver) is a parallel interior point
code that exploits any special structure in the Hessian and Jacobian matrices,
developed by Jacek Gondzio, Andreas Grothey and Robert Sarkissian.

The solver is an implementation of the primal-dual interior point method with multiple centrality correctors. The solver is implemented using object-oriented programming techniques. It solves linear (LP), quadratic (QP) and nonlinear (NLP) problems. The sequential code compares favourably with commercial packages and the parallel code shows perfect speed-ups.

*pdNCG* (primal-dual Newton Conjugate Gradients) is a MATLAB
implementation for the solution of L1-regularized strongly-convex problems.
The solver is memoryless, it requires only matrix-vector product operations,
hence it is appropriate for large-scale instances.

*quNac* (quasi-Newton action constrained) is a MATLAB implementation of
a Newton-CG method that uses quNac matrices as preconditioners, as described
in this article. The software was
written by Robert Gower.

*S2GD* is an efficient implementation of Semi Stochastic Gradient
Descent for logistic regression by Jakub Konečný.

*SML* (Structured Modelling Language) is an implementation of a
structure-conveying extension to the AMPL modelling language.

SML extends AMPL with object-oriented features that allow to construct models from sub-models. Unlike traditional modelling languages, the new approach does not scramble the block structure of the problem, and thus it enables the passing of this structure on to the solver. Its design allows the problem generation phase to be parallelisable.

*SVM-OOPS*, developed by Kristian Woodsend and Jacek Gondzio,
is for training standard linear support vector machines. It uses the
*OOPS* solver, and incorporates SVM reformulations designed for
interior point methods.

The algorithm is best suited for large-scale classification problems, where the number of samples greatly exceeds the number of features. It is highly efficient on multicore processors, solving problems involving tens or hundreds of thousands of support vectors. Larger problems can be solved on computing clusters using a combination of OpenMP and MPI. Our approach has the additional benefit of providing a good early approximation to the separating hyperplane.

More information about SVM-OOPS

Instance generators for l1-regularized over- and underdetermined least squares, written by Kimon Fountoulakis and Jacek Gondzio.