Hypatia.jl
Hypatia.jl copied to clipboard
interior point solver for general convex conic optimization problems
after #252
wlog we should be able to use use qr(P).Q instead of P (as opposed to only using the qr for point picking) leaving as a reminder although this is more...
Hypatia should support the `SOSPolynomialSet` set: https://github.com/JuliaOpt/SumOfSquares.jl/blob/master/src/sos_polynomial.jl#L16-L29 so that when using the Hypatia solver with SumOfSquares, the PSD reformulation is not done.
need to first rewrite loop code to not depend on index incrementing. need to be wary of threading loops with lapack calls inside due to https://github.com/JuliaLang/julia/issues/32786
If these will contain data from previous solve, it's possible they could cause failures during the next solve
shouldn't be too hard given the infrastructure already built. should involve some refactoring of nativeinterface.jl. based on CVXOPT's "cone QP" algorithm and the similar method in ConicIP.jl solver
through model or through solver?
in particular: - cone barrier procedures (one per primitive cone) - for WSOS cones specifically, the n+1 indepedent cholesky factorizations - ~~linear system solves (two independent sub-solves)~~ LAPACK does it...