gb_geo_reduce()

void gb2_comp::gb_geo_reduce ( gbvector *& f, gbvector *& fsyz )

private

Definition at line 397 of file res-a2-gb.cpp.

{
  gbvector head;
  gbvector *result = &head;
  result->next = nullptr;
 
  exponents_t div_totalexp = newarray_atomic(int, M->n_vars());
  int count = 0;
 
  gbvectorHeap fb(GR, F);
  gbvectorHeap fsyzb(GR, Fsyz);
  fb.add(f);
  fsyzb.add(fsyz);
  const gbvector *lead;
  while ((lead = fb.get_lead_term()) != nullptr)
    {
      Bag *b;
      GR->gbvector_get_lead_exponents(F, lead, div_totalexp);
      if (originalR->is_quotient_ring() &&
          originalR->get_quotient_monomials()->search_expvector(div_totalexp,
                                                                b))
        {
          const gbvector *g = originalR->quotient_gbvector(b->basis_elem());
          GR->reduce_lead_term_heap(F,
                                    Fsyz,
                                    lead,
                                    div_totalexp,  // are these two needed
                                    head.next,
                                    fb,
                                    fsyzb,
                                    g,
                                    nullptr);
          count++;
        }
      else if (monideals[lead->comp]->mi_search->search_expvector(div_totalexp,
                                                                  b))
        {
          gb_elem *q = reinterpret_cast<gb_elem *>(b->basis_ptr());
          GR->reduce_lead_term_heap(
              F, Fsyz, lead, div_totalexp, head.next, fb, fsyzb, q->f, q->fsyz);
          count++;
        }
      else
        {
          result->next = fb.remove_lead_term();
          result = result->next;
          result->next = nullptr;
        }
    }
 
  if (M2_gbTrace >= 4)
    {
      buffer o;
      o << "." << count;
      emit_wrapped(o.str());
    }
  f = head.next;
 
  fsyz = fsyzb.value();
  freemem(div_totalexp);
}

References gbvectorHeap::add(), int_bag::basis_elem(), int_bag::basis_ptr(), gbvector::comp, emit_wrapped(), F, gb_elem::f, freemem(), Fsyz, gb_elem::fsyz, gbvectorHeap::get_lead_term(), GR, M, M2_gbTrace, newarray_atomic, gbvector::next, originalR, gbvectorHeap::remove_lead_term(), result(), buffer::str(), and gbvectorHeap::value().

Referenced by gb_reduce().