Hypatia is a highly customizable open source interior point solver for generic conic optimization problems, written in Julia.
For more information on Hypatia, please see:
- documentation for Hypatia's conic form, predefined cones, and interfaces
- cones reference for cone definitions and oracles
- examples folder for applied examples and instances
- benchmarks folder for scripts used to run and analyze various computational benchmarks
and preprints of our papers:
- Solving natural conic formulations with Hypatia.jl for computational arguments for expanding the class of cones recognized by conic solvers
- Performance enhancements for a generic conic interior point algorithm for a description of Hypatia's algorithm and our enhanced stepping procedures
- Sum of squares generalizations for conic sets for barriers and computational techniques for our generalized polynomial sum of squares cones
- Conic optimization with spectral functions on Euclidean Jordan algebras for barriers and computational techniques for many of our epigraph/hypograph cones
and corresponding raw results CSV files generated by our run scripts in the benchmarks folder.
License
Hypatia is licensed under the MIT License (see LICENSE).
Installation
To use Hypatia, install Julia, then at the Julia REPL, type:
using Hypatia
using Pkg
Pkg.add("Hypatia")Hypatia is an experimental solver and a work in progress, and may not run with older releases of Julia. Default options/parameters are not well-tuned, so we encourage you to experiment with these.
Usage
Hypatia can be accessed through a low-level native Julia interface or through open-source modeling tools such as JuMP and Convex.jl. The native interface is more expressive, allowing Hypatia to solve conic models expressed with generic real floating point types and structured matrices or linear operators, for example. However, it is typically sufficient and more convenient to use JuMP.
Using JuMP, we can model a simple D-optimal experiment design problem and call Hypatia:
using LinearAlgebra
using JuMP
using Hypatia
model = Model(() -> Hypatia.Optimizer(verbose = false))
@variable(model, x[1:3] >= 0)
@constraint(model, sum(x) == 5)
@variable(model, hypo)
@objective(model, Max, hypo)
V = rand(2, 3)
Q = V * diagm(x) * V'
aff = vcat(hypo, [Q[i, j] for i in 1:2 for j in 1:i]...)
@constraint(model, aff in MOI.RootDetConeTriangle(2))
# solve and query solution
optimize!(model)
termination_status(model)
objective_value(model)
value.(x)See our D-optimal design example for more information and references.
Many more examples using the native interface or JuMP can be found in the examples folder.
Contributing
Comments, questions, suggestions, and improvements/extensions to the code or documentation are welcomed. Please reach out on Discourse, or submit an issue or contribute a PR on our GitHub. If contributing code, try to maintain consistent style and add docstrings or comments for clarity. New examples are welcomed and should be implemented similarly to the existing examples.
Acknowledgements
This work has been partially funded by the National Science Foundation under grant OAC-1835443 and the Office of Naval Research under grant N00014-18-1-2079.
Citing Hypatia
If you find Hypatia solver useful, please cite our solver paper:
@article{coey2022solving,
title={Solving natural conic formulations with {H}ypatia.jl},
author={Chris Coey and Lea Kapelevich and Juan Pablo Vielma},
year={2022},
journal={INFORMS Journal on Computing},
publisher={INFORMS},
volume={34},
number={5},
pages={2686--2699},
doi={https://doi.org/10.1287/ijoc.2022.1202}
}If you find aspects of Hypatia's IPM implementation useful, please cite our algorithm paper:
@article{coey2022performance,
title={Performance enhancements for a generic conic interior point algorithm},
author={Chris Coey and Lea Kapelevich and Juan Pablo Vielma},
year={2023},
journal={Mathematical Programming Computation},
publisher={Springer},
volume={15},
pages={53--101},
doi={https://doi.org/10.1007/s12532-022-00226-0}
}