libqsopt_ex.cpp

 
#include "libqsopt_ex.h"
 
namespace dlinear::qsopt_ex {
 
mpq_class *StringToMpqPtr(const std::string &str) { return CStringToMpqPtr(str.c_str()); }
mpq_class StringToMpq(const std::string &str) { return CStringToMpq(str.c_str()); }
mpq_class *CStringToMpqPtr(const char str[]) {
  mpq_t val;
  mpq_init(val);
  mpq_EGlpNumReadStr(val, str);
  auto result = new mpq_class(val);
  mpq_clear(val);
  return result;
}
mpq_class CStringToMpq(const char str[]) {
  mpq_t val;
  mpq_init(val);
  mpq_EGlpNumReadStr(val, str);
  mpq_class result(val);
  mpq_clear(val);
  return result;
}
 
void MpqArray::AllocateMpqArray(size_t n_elements) {
  auto const memSize = static_cast<size_t>(sizeof(mpq_t) * n_elements + sizeof(size_t));
  void *newArray = nullptr;
  if (memSize) {
    newArray = calloc(1, memSize);
    if (!newArray) {
      fprintf(stderr, "EXIT: Not enough memory while allocating %zd bytes", memSize);
      exit(1);
    }
  }
  size_t *sizeArray = n_elements ? static_cast<size_t *>(newArray) : nullptr;
  if (n_elements) sizeArray[0] = n_elements;
 
  array_ = reinterpret_cast<mpq_t *>(n_elements ? (sizeArray + 1) : nullptr);
  for (size_t i = 0; i < n_elements; ++i) mpq_init(array_[i]);
}
 
void MpqArray::FreeMpqArray() {
  auto *sizeArray = reinterpret_cast<size_t *>(array_);
  if (sizeArray) sizeArray--;
  size_t nElements = sizeArray ? sizeArray[0] : 0;
 
  for (size_t i = 0; i < nElements; ++i) mpq_clear(array_[i]);
  free(sizeArray);
  array_ = nullptr;
}
 
MpqArray::MpqArray(size_t n_elements) : array_{nullptr} { AllocateMpqArray(n_elements); }
 
MpqArray::~MpqArray() { FreeMpqArray(); }
 
void MpqArray::Resize(size_t nElements) {
  {
    FreeMpqArray();
    AllocateMpqArray(nElements);
  }
}
 
void QSXStart() { QSexactStart(); }
 
void QSXFinish() { QSexactClear(); }
 
}  // namespace dlinear::qsopt_ex