res-poly-ring.cpp

Go to the documentation of this file.

/* Copyright 2015-2021, Michael E. Stillman */
 
#include "schreyer-resolution/res-poly-ring.hpp"
#include "schreyer-resolution/res-monomial-sorter.hpp"  // for ResMonomialSorter
struct ResSchreyerOrder;
 
long ResPolynomial::npoly_destructor = 0;
long ResPolynomialConstructor::ncalls = 0;
long ResPolynomialConstructor::ncalls_fromarray = 0;
 
void ResPolyRing::memUsage(const ResPolynomial& f,
                           long& nterms,
                           long& bytes_used,
                           long& bytes_alloc) const
{
  long sz = 0;
  // sz = f.len * sizeof(FieldElement); // TODO MES: what to add here?
  sz += f.monoms.size() * sizeof(res_monomial_word);
  nterms += f.len;
  bytes_used += sz;
  bytes_alloc += sz;
}
 
#if 0
bool check_poly(const ResPolyRing& R,
                const ResPolynomial& f,
                const ResSchreyerOrder& ord)
{
  // We loop through each monomial, checking it against the one before
  // The order used is the Schreyer order given.
  auto& M = R.monoid();
  ResPolynomialIterator i(R, f);
  ResPolynomialIterator end(R, f, 1);
  res_const_packed_monomial prev = nullptr;
  for (; i != end; ++i)
    {
      if (prev == nullptr)
        prev = i.monomial();
      else
        {
          // Now compare to previous monomial
          long comp1 = M.get_component(prev);
          long comp2 = M.get_component(i.monomial());
          int result = M.compare_schreyer(prev,
                                          i.monomial(),
                                          ord.mTotalMonom[comp1],
                                          ord.mTotalMonom[comp2],
                                          ord.mTieBreaker[comp1],
                                          ord.mTieBreaker[comp2]);
          if (result == EQ or result == GT)
            {
              return false;
            }
          prev = i.monomial();
        }
    }
  return true;
}
#endif
 
bool check_poly(const ResPolyRing& R,
                const ResPolynomial& f,
                const ResSchreyerOrder& ord)
{
  std::vector<int*> monoms;
  auto end = ResPolynomialIterator(R,f,1);
  for (auto i = ResPolynomialIterator(R,f); i != end; ++i)
    {
      monoms.push_back(const_cast<res_packed_monomial>(i.monomial()));
    }
  ResMonomialSorter S(R.originalMonoid(), R.monoid(), ord, monoms);
  return S.ordered();
}
 
// Local Variables:
// compile-command: "make -C $M2BUILDDIR/Macaulay2/e "
// indent-tabs-mode: nil
// End: