geo_remainder()

void MarkedGB::geo_remainder	(	gbvector *&	f,
		bool	use_denom,
		ring_elem &	denom )

Definition at line 224 of file reducedgb-marked.cpp.

{
  gbvector head;
  gbvector *frem = &head;
  (void) use_denom;
  (void) denom;
  frem->next = nullptr;
 
  gbvectorHeap fb(R, F);
  gbvectorHeap zero(R, Fsyz);
  fb.add(f);
 
  const gbvector *lead;
  exponents_t EXP = ALLOCATE_EXPONENTS(R->exponent_byte_size());
  while ((lead = fb.get_lead_term()) != nullptr)
    {
      R->gbvector_get_lead_exponents(F, lead, EXP);
      int x = lead->comp;
      int w = T->find_divisor(EXP, x);
      if (w < 0)
        {
          frem->next = fb.remove_lead_term();
          frem = frem->next;
          frem->next = nullptr;
        }
      else
        {
          POLY g = polys[w];
          R->reduce_marked_lead_term_heap(
              F, Fsyz, lead, EXP, head.next, fb, zero, leadterms[w], g.f, nullptr);
        }
    }
  f = head.next;
  R->gbvector_sort(F, f);
}

References gbvectorHeap::add(), ALLOCATE_EXPONENTS, gbvector::comp, ReducedGB::F, POLY::f, ReducedGB::Fsyz, gbvectorHeap::get_lead_term(), leadterms, gbvector::next, POLY, ReducedGB::R, gbvectorHeap::remove_lead_term(), T, x, and zero.