pp_Mult_mm_Noether__T.cc File Reference

Go to the source code of this file.

Functions
LINKAGE poly	pp_Mult_mm_Noether__T (poly p, const poly m, const poly spNoether, int &ll, const ring ri)

Function Documentation

pp_Mult_mm_Noether__T()

LINKAGE poly pp_Mult_mm_Noether__T	(	poly	p,
		const poly	m,
		const poly	spNoether,
		int &	ll,
		const ring	ri
	)

Definition at line 19 of file pp_Mult_mm_Noether__T.cc.

{
  p_Test(p, ri);
  p_LmTest(m, ri);
  assume(spNoether != NULL);
  if (p == NULL)
  {
    ll = 0;
    return NULL;
  }
  spolyrec rp;
  poly q = &rp, r;
  number n;
  const unsigned long *spNoether_exp = spNoether->exp;
  number ln = pGetCoeff(m);
  omBin bin = ri->PolyBin;
  DECLARE_LENGTH(const unsigned long length = ri->ExpL_Size);
  DECLARE_ORDSGN(const long* ordsgn = ri->ordsgn);
  const unsigned long* m_e = m->exp;
  pAssume(!n_IsZero__T(ln,ri->cf));
  pAssume1(p_GetComp(m, ri) == 0 || p_MaxComp(p, ri) == 0);
  int l = 0;
 
  do
  {
    p_AllocBin(r, bin, ri);
    p_MemSum__T(r->exp, p->exp, m_e, length);
    p_MemAddAdjust__T(r, ri);
 
    p_MemCmp__T(r->exp, spNoether_exp, length, ordsgn, goto Continue, goto Continue, goto Break);
 
    Break:
    p_FreeBinAddr(r, ri);
    break;
 
    Continue:
 
    n = n_Mult__T(ln, pGetCoeff(p), ri->cf);
 
    #ifdef HAVE_RINGS
    if(n_IsZero__T(n, ri->cf))
    {
      n_Delete__T(&n, ri->cf);
      p_FreeBinAddr(r, ri);
    } else
    #endif
    {
      l++;
      q = pNext(q) = r;
      pSetCoeff0(q, n);
    }
 
    pIter(p);
 
  } while (p != NULL);
 
  if (ll < 0)
    ll = l;
  else
    ll = pLength(p);
 
  pNext(q) = NULL;
 
  p_Test(pNext(&rp), ri);
  return pNext(&rp);
}