OSQP: An Operator Splitting Solver for Quadratic Programs

B. Stellato, G. Banjac, P. J. Goulart, A. Bemporad and S. Boyd

Mathematical Programming Computation, February 2020.
BibTeX  URL  Preprint  Code 

@article{SBGetal:2020,
  author = {B. Stellato and G. Banjac and P. J. Goulart and A. Bemporad and S. Boyd},
  title = {OSQP: An Operator Splitting Solver for Quadratic Programs},
  journal = {Mathematical Programming Computation},
  year = {2020},
  url = {https://doi.org/10.1007/s12532-020-00179-2},
  doi = {10.1007/s12532-020-00179-2}
}

We present a general purpose solver for quadratic programs based on the alternating direction method of multipliers, employing a novel operator splitting technique that requires the solution of a quasi-definite linear system with the same coefficient matrix in each iteration. Our algorithm is very robust, placing no requirements on the problem data such as positive definiteness of the objective function or linear independence of the constraint functions. It is division-free once an initial matrix factorization is carried out, making it suitable for real-time applications in embedded systems. In addition, our technique is the first operator splitting method for quadratic programs able to reliably detect primal and dual infeasible problems from the algorithm iterates. The method also supports factorization caching and warm starting, making it particularly efficient when solving parametrized problems arising in finance, control, and machine learning. Our open-source C implementation OSQP has a small footprint, is library-free, and has been extensively tested on many problem instances from a wide variety of application areas. It is typically ten times faster than competing interior point methods, and sometimes much more when factorization caching or warm start is used.