gbvector_normal_form() [2/2]

void QRingInfo_field_QQ::gbvector_normal_form ( const FreeModule * F, gbvector *& f, bool use_denom, ring_elem & denom ) const

virtual

Reimplemented from QRingInfo.

Definition at line 275 of file qring.cpp.

{
  exponents_t EXP1 = ALLOCATE_EXPONENTS(exp_size);
 
  GBRing *GR = R->get_gb_ring();
  gbvector head;
  gbvector *result = &head;
  result->next = nullptr;
  gbvector *t = f;
  while (t != nullptr)
    {
      GR->gbvector_get_lead_exponents(F, t, EXP1);
      int x = ringtable->find_divisor(EXP1, 1);
      if (x >= 0)
        {
          const gbvector *r = quotient_gbvector(x);
          gbvector *zero = nullptr;
          GR->gbvector_reduce_lead_term(
              F, F, head.next, t, zero, r, zero, use_denom, denom);
        }
      else
        {
          result->next = t;
          t = t->next;
          result = result->next;
          result->next = nullptr;
        }
    }
  f = head.next;
}

References ALLOCATE_EXPONENTS, QRingInfo::exp_size, GBRing::gbvector_get_lead_exponents(), GBRing::gbvector_reduce_lead_term(), gbvector::next, QRingInfo::quotient_gbvector(), QRingInfo::R, result(), QRingInfo_field::ringtable, x, and zero.