Changeset e39e193 in git

Timestamp:

Nov 14, 2018, 5:40:40 PM (5 years ago)

Author:

Karim Abou Zeid <karim23697@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

93386435caf06320aebf4d6af894bff6839db8b4

Parents:

ecdeb553994eac8047af1b1065f7154fe0c5a0c5

Message:

Move divides and varAt to freegb

Location:

Singular

Files:

• LIB/fpaprops.lib (modified) (4 diffs)
• LIB/freegb.lib (modified) (2 diffs)
• extra.cc (modified) (1 diff)

Singular/LIB/fpaprops.lib

@@ -767 +767 @@
         w_overlap = w_overlap[1 .. (l-2)*lV];
       }
-      if (v_overlap == w_overlap && !lpdivides(G, v * lpVarAt(w, l - 1))) {
+      if (v_overlap == w_overlap && !lpLmDivides(G, v * lpVarAt(w, l - 1))) {
         UG[i,j] = 1;
       }
@@ -808 +808 @@
 }
 
-static proc lpVarAt(poly p, int pos) {
-  return (system("lpVarAt", p, pos));
-}
-
 static proc lpStandardWords(ideal G, int length)
 "ASSUME: G is simplified
@@ -822 +818 @@
   ideal words = maxideal(length);
   for (int i = ncols(words); i >= 1; i--) {
-    if (lpdivides(G, words[i])) {
+    if (lpLmDivides(G, words[i])) {
       words = delete(words, i);
     }
@@ -867 +863 @@
   } kill i;
   return (words);
-}
-
-proc lpdivides(ideal I, poly p) {
-  for (int i = 1; i <= ncols(I); i++) {
-    if (system("lpdivides", I[i], p)) {
-      return (1);
-    }
-  } kill i;
-  return (0);
 }
 

Singular/LIB/freegb.lib

@@ -32 +32 @@
 
 isVar(p);                        check whether p is a power of a single variable
+
+lpLmDivides(ideal I, poly p)     tests if there is a polynomial q in I with LM(q)|LM(p)
+lpVarAt(poly p, int pos)         returns the variable (as a poly) at position pos of the poly p
 
 SEE ALSO: fpadim_lib, fpaprops_lib, fpalgebras_lib, LETTERPLACE
@@ -354 +357 @@
   poly i = 1;
   isVar(i);
+}
+
+proc lpLmDivides(ideal I, poly p)
+"USAGE: lpLmDivides(I); I an ideal
+RETURN: boolean
+ASSUME: basering is a Letterplace ring
+PURPOSE: tests if there is a polynomial q in I with LM(q)|LM(p)
+EXAMPLE: example lpLmDivides; shows examples
+"
+{
+  for (int i = 1; i <= ncols(I); i++) {
+    if (system("lpLmDivides", I[i], p)) {
+      return (1);
+    }
+  } kill i;
+  return (0);
+}
+example
+{
+  "EXAMPLE:"; echo = 2;
+  ring r = 0,(x,y),dp;
+  def R = makeLetterplaceRing(5);
+  setring R;
+  poly p = x*y*y;
+  lpLmDivides(y*y, p);
+  lpLmDivides(y*x, p);
+  lpLmDivides(ideal(y*y, y*x), p);
+}
+
+proc lpVarAt(poly p, int pos)
+"USAGE: lpVarAt(p, pos); p a poly, pos an int
+RETURN: poly
+ASSUME: basering is a Letterplace ring
+PURPOSE: returns the variable (as a poly) at position pos of the poly p
+EXAMPLE: example lpVarAt; shows examples
+"
+{
+  return (system("lpVarAt", p, pos));
+}
+example
+{
+  "EXAMPLE:"; echo = 2;
+  ring r = 0,(x,y),dp;
+  def R = makeLetterplaceRing(5);
+  setring R;
+  poly p = y*y*x;
+  lpVarAt(p, 3);
 }
 

Singular/extra.cc

@@ -1211 +1211 @@
   /*==================== divide-test for freeGB  =================*/
   #ifdef HAVE_SHIFTBBA
-    if (strcmp(sys_cmd, "lpdivides") == 0)
+    if (strcmp(sys_cmd, "lpLmDivides") == 0)
     {
       const short t[]={2,POLY_CMD,POLY_CMD};