Changeset e645db in git for Singular/LIB/dmodapp.lib

Timestamp:

Apr 18, 2020, 1:52:40 AM (4 years ago)

Author:

Viktor Levandovskyy <levandov@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'c5facdfddea2addfd91babd8b9019161dea4b695')

Children:

933ed85f7a4b64fb6b04964c04a137915354e36b

Parents:

57fbc7a35840961d9ffe380d87f3cc30a4574b66161411bdc5cc30f7965d2c36a5250edcb33d293e

Message:

Merge remote-tracking branch 'upstream/spielwiese' into spielwiese

File:

• Singular/LIB/dmodapp.lib (modified) (4 diffs)

Singular/LIB/dmodapp.lib

@@ -140 +140 @@
   - changed charVariety: no homogenization is needed
   - changed inForm: code is much simpler using jet
+  
+  06.04.20 by VL:
+  - engine updated with 2 more modular choices
 
 */
@@ -153 +156 @@
 LIB "qhmoduli.lib"; // for Max
 LIB "sing.lib";     // for slocus
+LIB "ncModslimgb.lib"; // for modular noncomm GB
 
 
@@ -263 +267 @@
 
 proc engine(def I, int i)
+"USAGE:  engine(I,i);  I ideal/module/matrix, i an int
+RETURN:  the same type as I
+PURPOSE: compute the Groebner basis of I with the algorithm, chosen via i
+NOTE:    By default and if i=0, slimgb is used; 
+if i=1, std does the job;
+if i=2, modular slimgb is run without final exactness check;
+if i=3, modular slimgb is run with final exactness check.
+EXAMPLE: example engine; shows examples
+"
+{
+  /* std - slimgb mix */
+  def J;
+  //  ideal J;
+  if ((i<=0) || (i>=4))
+  {
+      i = 0;
+      J = slimgb(I);
+  }
+  if (i==1)
+  {
+    // without options -> strange! (ringlist?)
+    intvec v = option(get);
+    option(redSB);
+    option(redTail);
+    J = std(I);
+    option(set, v);
+  }
+  if (i==2)
+  {
+      // modular slimgb without final exactness check
+      J = ncmodslimgb(I,0);
+  }
+  if (i==3)
+  {
+      // modular slimgb with final exactness check
+      J = ncmodslimgb(I,1);
+  }  
+  return(J);
+}
+example
+{
+  "EXAMPLE:"; echo = 2;
+  ring r = 0,(x,y),Dp;
+  ideal I  = y*(x3-y2),x*(x3-y2);
+  engine(I,0); // uses slimgb
+  engine(I,1); // uses std
+  engine(I,2); // uses modular slimgb without exactness check
+  engine(I,3); // uses modular slimgb with exactness check
+}
+
+/* 
+proc engine(def I, int i)
 "USAGE:  engine(I,i);  I  ideal/module/matrix, i an int
 RETURN:  the same type as I
@@ -295 +351 @@
   engine(I,0); // uses slimgb
   engine(I,1); // uses std
-}
+} 
+*/
 
 proc poly2list (poly f)