source: git/kernel/pp_Mult_mm__T.cc @ 35aab3

spielwiese
Last change on this file since 35aab3 was 35aab3, checked in by Hans Schönemann <hannes@…>, 20 years ago
This commit was generated by cvs2svn to compensate for changes in r6879, which included commits to RCS files with non-trunk default branches. git-svn-id: file:///usr/local/Singular/svn/trunk@6880 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 1.4 KB
Line 
1/****************************************
2*  Computer Algebra System SINGULAR     *
3****************************************/
4/***************************************************************
5 *  File:    pp_Mult_mm__Template.cc
6 *  Purpose: template for p_Mult_n
7 *  Author:  obachman (Olaf Bachmann)
8 *  Created: 8/00
9 *  Version: $Id: pp_Mult_mm__T.cc,v 1.1.1.1 2003-10-06 12:16:02 Singular Exp $
10 *******************************************************************/
11
12/***************************************************************
13 *
14 *   Returns:  p*m
15 *   Const:    p, m
16 *
17 ***************************************************************/
18LINKAGE poly pp_Mult_mm(poly p, const poly m, const ring ri, poly &last)
19{
20  p_Test(p, ri);
21  p_LmTest(m, ri);
22  if (p == NULL)
23  {
24    last = NULL;
25    return NULL;
26  }
27  spolyrec rp;
28  poly q = &rp, r;
29  number ln = pGetCoeff(m);
30  omBin bin = ri->PolyBin;
31  DECLARE_LENGTH(const unsigned long length = ri->ExpL_Size);
32  const unsigned long* m_e = m->exp;
33  pAssume(!n_IsZero(ln,ri));
34  pAssume1(p_GetComp(m, ri) == 0 || p_MaxComp(p, ri) == 0);
35
36
37  do
38  {
39    p_AllocBin( pNext(q), bin, ri);
40    q = pNext(q);
41    pSetCoeff0(q, n_Mult(ln, pGetCoeff(p), ri));
42    p_MemSum(q->exp, p->exp, m_e, length);
43    p_MemAddAdjust(q, ri);
44    p = pNext(p);
45  }
46  while (p != NULL);
47  last = q;
48  pNext(q) = NULL;
49
50  p_Test(pNext(&rp), ri);
51  return pNext(&rp);
52}
53
54
Note: See TracBrowser for help on using the repository browser.