Changeset b71400a in git for Singular/dyn_modules/gfanlib/initial.cc

Timestamp:

Sep 15, 2014, 11:08:21 AM (10 years ago)

Author:

Yue Ren <ren@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

dffd1541cc74564f509bd5d284948852484430b0

Parents:

9c03260ac0404405e4994f73b30ef46f079535d8

git-author:

Yue Ren <ren@mathematik.uni-kl.de>2014-09-15 11:08:21+02:00

git-committer:

Yue Ren <ren@mathematik.uni-kl.de>2015-02-06 13:47:05+01:00

Message:

chg: status update 15.09.

File:

• Singular/dyn_modules/gfanlib/initial.cc (modified) (3 diffs)

Singular/dyn_modules/gfanlib/initial.cc

@@ -36 +36 @@
   for (int i=0; i<W.getHeight(); i++)
     d[i] = wDeg(p,r,W[i]);
+  return d;
+}
+
+gfan::ZVector WDeg(const poly p, const ring r, const gfan::ZVector w, const gfan::ZMatrix W)
+{
+  gfan::ZVector d = gfan::ZVector(W.getHeight()+1);
+  d[0] = wDeg(p,r,w);
+  for (int i=0; i<W.getHeight(); i++)
+    d[i+1] = wDeg(p,r,W[i]);
   return d;
 }
@@ -166 +175 @@
 }
 
+poly initial(const poly p, const ring r, const gfan::ZVector w, const gfan::ZMatrix W)
+{
+  int n = rVar(r);
+  poly q0 = p_Head(p,r);
+  poly q1 = q0;
+  gfan::ZVector d = WDeg(p,r,w,W);
+  for (poly currentTerm = p->next; currentTerm; pIter(currentTerm))
+  {
+    gfan::ZVector e = WDeg(currentTerm,r,w,W);
+    if (d<e)
+    {
+      p_Delete(&q0,r);
+      q0 = p_Head(p,r);
+      q1 = q0;
+      d = e;
+    }
+    else
+      if (d==e)
+      {
+        pNext(q1) = p_Head(currentTerm,r);
+        pIter(q1);
+      }
+  }
+  return q0;
+}
+
+
 /***
  * Runs the above procedure over all generators of an ideal.
@@ -179 +215 @@
 }
 
+ideal initial(const ideal I, const ring r, const gfan::ZVector w, const gfan::ZMatrix W)
+{
+  int k = idSize(I); ideal inI = idInit(k);
+  for (int i=0; i<k; i++)
+    inI->m[i] = initial(I->m[i],r,w,W);
+  return inI;
+}
 
 #ifndef NDEBUG