1
$\begingroup$

I'm implementing some algorithms in C++ and I need to solve some linear programs as an intermediate step within the algorithm.

I've managed to get the linear_solver from google or-ools working, as well as glpk.

My issue is that I need to work with and solve LPs with very high precision integer types - more than 256 bits - in order for the algorithm to be correct (it's complicated but the algorithm can give very incorrect answers if less precision is used). At the moment I'm using mpz_class and mpf_class from gmp to store my numbers.

Improve this question
$\endgroup$
1
  • $\begingroup$ I don't know an answer to your question, but just wanted to mention that although increased precision has its uses, we typically don't see as much benefit as one might expect - it will give you a handful more orders of magnitude of correct calculations and that's it. Mileage may very, but people are typically far better off reformulating the problem, trying to scale it well, and/or using a solver with good numerics. $\endgroup$
    – Nikos Kazazakis
    Commented Feb 12 at 10:13

2 Answers 2

Reset to default
1
$\begingroup$

The Soplex solver is a simplex-based LP solver, available as a C++ library. It supports regular (double precision) input, quad precision computation, and can provide exact solutions using iterative refinement and rational exact arithmetic.

While I have never used it myself directly, it is the default LP solver in SCIP. I would recommend checking the documentation to see how you can use the higher-precision functionality.

Improve this answer
$\endgroup$
1
  • $\begingroup$ This page details how to use Soplex as an exact solver: link. To use higher precision Soplex, I believe one needs GMP as well. $\endgroup$
    – J. Dionisio
    Commented Jul 13 at 9:21
0
$\begingroup$

A similar question has been asked here. The answers reference a number of exact solvers. I don't know if any solver supports higher but limited precision instead of exact.

GLPK itself supports exact arithmetic. It may be more difficult to use from the C++ interface than the command line, but is documented in the API.

Improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.