Changeset 69bc59 in git for kernel/ringgb.cc

Timestamp:

Jun 8, 2006, 11:56:54 PM (18 years ago)

Author:

Oliver Wienand <wienand@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

8e4858985c8a175c0ae6f5422cd1c577e587236d

Parents:

728a2c96465e8ac52df53f624fc9319cc328b4e0

Message:

[oliver]

ringgb.*:
* routines for testing of GB


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

File:

• kernel/ringgb.cc (modified) (5 diffs)

kernel/ringgb.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: ringgb.cc,v 1.7 2006-06-08 18:00:35 wienand Exp $ */
+/* $Id: ringgb.cc,v 1.8 2006-06-08 21:56:54 wienand Exp $ */
 /*
 * ABSTRACT: ringgb interface
@@ -120 +120 @@
 }
 
+poly plain_spoly(poly f, poly g) {
+  number cf = pGetCoeff(f), cg = pGetCoeff(g);
+  int ct = ksCheckCoeff(&cf, &cg); // gcd and zero divisors
+  poly fm, gm;
+  k_GetLeadTerms(f, g, currRing, fm, gm, currRing);
+  pSetCoeff0(fm, cg);
+  pSetCoeff0(gm, cf);  // and now, m1 * LT(p1) == m2 * LT(p2)
+  return(pSub(pMult_mm(f, fm), pMult_mm(g, gm)));
+}
+
+
 poly spolyRing2toM(poly f, poly g, ring r) {
   poly m1 = NULL;
@@ -127 +138 @@
   // printPolyMsg("m2=", m2, "");
   // PrintLn();
-  return pSub(pp_Mult_mm(f, m1, r), pp_Mult_mm(g, m2, r));
+  return pSub(p_Mult_mm(f, m1, r), p_Mult_mm(g, m2, r));
 }
 
@@ -143 +154 @@
     // wrp(G->m[i]);
     // PrintLn();
-    h = spolyRing2toM(h, G->m[i], r);
+    h = spolyRing2toM(h, pCopy(G->m[i]), r);
     // PrintS("=> h=");
     // wrp(h);
@@ -174 +185 @@
   return h;
 }
+
+int testGB(ideal GI) {
+  poly f, g, h;
+  int i = 0;
+  int j = 0;
+  for (i = 0; i < IDELEMS(GI) - 1; i++) {
+    for (j = i + 1; j < IDELEMS(GI); j++) {
+      f = pCopy(GI->m[i]);
+      g = pCopy(GI->m[j]);
+      h = plain_spoly(f, g);
+      if (ringNF(h, GI, currRing) != NULL) {
+        wrp(GI->m[i]);
+        PrintLn();
+        wrp(GI->m[j]);
+        PrintLn();
+        wrp(h);
+        PrintLn();
+        wrp(ringNF(h, GI, currRing));
+        PrintLn();
+        return(0);
+      }
+      pDelete(&h);
+    }
+  }
+  return(1);
+}
+
 #endif