find_smallest_coeff_divisor()

int MonomialTableZZ::find_smallest_coeff_divisor	(	exponents_t	exp,
		int	comp ) const

Definition at line 247 of file montableZZ.cpp.

{
  assert(comp >= 1);
  if (comp >= static_cast<int>(_head.size())) return -1;
  mon_term *head = _head[comp];
  mon_term *t;
 
  int smallest_val = -1;
  mpz_srcptr smallest = nullptr;
 
  unsigned long expmask = ~(monomial_mask(_nvars, exp));
 
  for (t = head->_next; t != head; t = t->_next)
    if ((expmask & t->_mask) == 0)
      {
        bool is_div = 1;
        for (int i = 0; i < _nvars; i++)
          if (exp[i] < t->_lead[i])
            {
              is_div = 0;
              break;
            }
        if (is_div)
          {
            if (smallest_val < 0 || (mpz_cmpabs(smallest, t->_coeff) > 0))
              {
                smallest_val = t->_val;
                smallest = t->_coeff;
              }
          }
      }
  return smallest_val;
}

References MonomialTableZZ::mon_term::_coeff, MonomialTableZZ::mon_term::_lead, MonomialTableZZ::mon_term::_mask, MonomialTableZZ::mon_term::_next, _nvars, MonomialTableZZ::mon_term::_val, and monomial_mask().