pp_Mult_mm__T.cc

Go to the documentation of this file.

/****************************************
*  Computer Algebra System SINGULAR     *
****************************************/
/***************************************************************
 *  File:    pp_Mult_mm__Template.cc
 *  Purpose: template for p_Mult_n
 *  Author:  obachman (Olaf Bachmann)
 *  Created: 8/00
 *******************************************************************/
 
/***************************************************************
 *
 *   Returns:  p*m
 *   Const:    p, m
 *
 ***************************************************************/
LINKAGE poly pp_Mult_mm__T(poly p, const poly m, const ring ri)
{
  p_Test(p, ri);
  p_LmTest(m, ri);
  if (p == NULL)
  {
    return NULL;
  }
  spolyrec rp;
#ifdef HAVE_ZERODIVISORS
  rp.next = NULL;
#endif
  poly q = &rp;
  number ln = pGetCoeff(m);
  omBin bin = ri->PolyBin;
  DECLARE_LENGTH(const unsigned long length = ri->ExpL_Size);
  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);
  number tmp;
 
  do
  {
    tmp = n_Mult__T(ln, pGetCoeff(p), ri->cf);
#ifdef HAVE_ZERODIVISORS
    if (! n_IsZero__T(tmp, ri->cf))
    {
#endif
      p_AllocBin( pNext(q), bin, ri);
      q = pNext(q);
      pSetCoeff0(q, tmp);
      p_MemSum__T(q->exp, p->exp, m_e, length);
      p_MemAddAdjust__T(q, ri);
#ifdef HAVE_ZERODIVISORS
    }
    else n_Delete__T(&tmp, ri->cf);
#endif
    p = pNext(p);
  }
  while (p != NULL);
  pNext(q) = NULL;
 
  p_Test(pNext(&rp), ri);
  return pNext(&rp);
}