Changeset a29995 in git for Singular/pp_Mult_mm_Noether__Template.cc

Timestamp:

Nov 28, 2000, 12:50:56 PM (23 years ago)

Author:

Olaf Bachmann <obachman@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')

Children:

718add5b6d1cfeccd8e65b50d42db0fd515f126a

Parents:

5e16f9ae5b8928fda665888d965340387f46c85c

Message:

* towards tailRings for local case


git-svn-id: file:///usr/local/Singular/svn/trunk@4777 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• Singular/pp_Mult_mm_Noether__Template.cc (modified) (2 diffs)

Singular/pp_Mult_mm_Noether__Template.cc

@@ -7 +7 @@
  *  Author:  obachman (Olaf Bachmann)
  *  Created: 8/00
- *  Version: $Id: pp_Mult_mm_Noether__Template.cc,v 1.1 2000-11-25 11:59:57 Singular Exp $
+ *  Version: $Id: pp_Mult_mm_Noether__Template.cc,v 1.2 2000-11-28 11:50:56 obachman Exp $
  *******************************************************************/
 
 /***************************************************************
  *
- *   Returns:  p*m
+ *   Returns:  p*m, ll
+ *   ll == pLength(p*m) , if on input ll < 0
+ *   pLength(p) - pLength(p*m), if on input ll >= 0
  *   Const:    p, m
- *   monoms which are smaller than spNoether are cut
  *
  ***************************************************************/
-poly pp_Mult_mm_Noether(poly p, const poly m, const poly spNoether, int &ll, const ring ri, poly &last)
+poly pp_Mult_mm_Noether(poly p, const poly m, const poly spNoether, int &ll, 
+                        const ring ri, poly &last)
 {
   p_Test(p, ri);
   p_LmTest(m, ri);
+  assume(spNoether != NULL);
   if (p == NULL) 
   {
+    ll = 0;
     last = NULL;
-    ll = 0;
     return NULL;
   }
   spolyrec rp;
   poly q = &rp, r;
+  const unsigned long *spNoether_exp = spNoether->exp;
   number ln = pGetCoeff(m);
   omBin bin = ri->PolyBin;
-  assume(spNoether != NULL);
-  unsigned long *spNoether_exp = spNoether->exp;
   DECLARE_LENGTH(const unsigned long length = ri->ExpL_Size);
   DECLARE_ORDSGN(const long* ordsgn = ri->ordsgn);
@@ -40 +42 @@
   int l = 0;
   
-  while (p != NULL)
+  do
   {
-    p_AllocBin( r, bin, ri);
+    p_AllocBin(r, bin, ri);
     p_MemSum(r->exp, p->exp, m_e, length);
     p_MemAddAdjust(r, ri);
-    p_MemCmp(r->exp, spNoether->exp, length, ordsgn, goto Continue, goto Continue, goto Break);
+
+    p_MemCmp(r->exp, spNoether_exp, length, ordsgn, goto Continue, goto Continue, goto Break);
     
     Break:
     p_FreeBinAddr(r, ri);
     break;
-    
+
     Continue:
+    l++;
     q = pNext(q) = r;
     pSetCoeff0(q, n_Mult(ln, pGetCoeff(p), ri));
-    l++;
     pIter(p);
-  }
-  if (q  != &rp) last = q;
+  } while (p != NULL);
+
+  if (ll < 0)
+    ll = l;
+  else
+    ll = pLength(p);
+  
+  if (q != &rp) 
+    last = q;
   pNext(q) = NULL;
-  ll = l;
 
   p_Test(pNext(&rp), ri);