p_Merge_q__T.cc File Reference

Go to the source code of this file.

Functions
LINKAGE poly	p_Merge_q__T (poly p, poly q, const ring r)

Function Documentation

p_Merge_q__T()

LINKAGE poly p_Merge_q__T	(	poly	p,
		poly	q,
		const ring	r
	)

Definition at line 17 of file p_Merge_q__T.cc.

{
  assume(p != NULL && q != NULL);
  p_Test(p, r);
  p_Test(q, r);
#if PDEBUG > 0
  int l = pLength(p) + pLength(q);
#endif
 
  spolyrec rp;
  poly a = &rp;
  DECLARE_LENGTH(const unsigned long length = r->CmpL_Size);
  DECLARE_ORDSGN(const long* ordsgn = r->ordsgn);
 
  Top:     // compare p and q w.r.t. monomial ordering
  p_MemCmp__T(p->exp, q->exp, length, ordsgn, goto Equal, goto Greater , goto Smaller);
 
  Equal:
  // should never get here
  dReportError("Equal monomials in p_Merge_q");
  return NULL;
 
  Greater:
  a = pNext(a) = p;
  pIter(p);
  if (p==NULL) { pNext(a) = q; goto Finish;}
  goto Top;
 
  Smaller:
  a = pNext(a) = q;
  pIter(q);
  if (q==NULL) { pNext(a) = p; goto Finish;}
  goto Top;
 
  Finish:
 
  p_Test(pNext(&rp), r);
#if PDEBUG > 0
  pAssume1(l - pLength(pNext(&rp)) == 0);
#endif
  return pNext(&rp);
}