CRA() [4/5]

vec ChineseRemainder::CRA ( const PolyRing * R, vec f, vec g, mpz_srcptr um, mpz_srcptr vn, mpz_srcptr mn )

static

Definition at line 181 of file cra.cpp.

{
  vecterm head;
  vec result = &head;
 
  while (1)
    {
      if (g == nullptr)
        {
          // mult each term of f by n:
          for (; f != nullptr; f = f->next)
            {
              result->next = R->new_vec();
              result = result->next;
              result->next = nullptr;
              result->coeff = CRA(R, f->coeff, nullptr, um, vn, mn);
              result->comp = f->comp;
            }
          break;
        }
      if (f == nullptr)
        {
          // mult each term of g by n:
          for (; g != nullptr; g = g->next)
            {
              result->next = R->new_vec();
              result = result->next;
              result->next = nullptr;
              result->coeff = CRA(R, nullptr, g->coeff, um, vn, mn);
              result->comp = g->comp;
            }
          break;
        }
      if (f->comp < g->comp)
        {
          result->next = R->new_vec();
          result = result->next;
          result->next = nullptr;
          result->coeff = CRA(R, nullptr, g->coeff, um, vn, mn);
          result->comp = g->comp;
          g = g->next;
        }
      else if (f->comp > g->comp)
        {
          result->next = R->new_vec();
          result = result->next;
          result->next = nullptr;
          result->coeff = CRA(R, f->coeff, nullptr, um, vn, mn);
          result->comp = f->comp;
          f = f->next;
        }
      else
        {
          result->next = R->new_vec();
          result = result->next;
          result->next = nullptr;
          result->coeff = CRA(R, f->coeff, g->coeff, um, vn, mn);
          result->comp = f->comp;
          f = f->next;
          g = g->next;
        }
    }
  result->next = nullptr;
  return head.next;
}

References CRA(), Ring::new_vec(), vec::next, and result().