Changeset af7145 in git for Singular

Timestamp:

Dec 1, 2010, 1:56:35 PM (13 years ago)

Author:

Frank Seelisch <seelisch@…>

Branches:

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

Children:

35715cbb715ef8f8977fca4b923054147968c0d1

Parents:

4d22d3bbe140ec25c1d800b92606f9067563c581

Message:

command in extra.cc for hensel factors of a bivariate poly

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

File:

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

Singular/extra.cc

@@ -137 +137 @@
 
 #include <kernel/shiftgb.h>
+#include <kernel/linearAlgebra.h>
 
 #ifdef HAVE_EIGENVAL
@@ -629 +630 @@
         }
       }
+  /*==================== Hensel's lemma ======================*/
+      if(strcmp(sys_cmd, "henselfactors")==0)
+      {
+        if ((h != NULL) && (h->Typ() == POLY_CMD) &&
+            (h->next != NULL) && (h->next->Typ() == POLY_CMD) &&
+            (h->next->next != NULL) && (h->next->next->Typ() == POLY_CMD) &&
+            (h->next->next->next != NULL) &&
+                                     (h->next->next->next->Typ() == INT_CMD))
+        {
+          poly hh = (poly)h->Data();
+          poly f0 = (poly)h->next->Data();
+          poly g0 = (poly)h->next->next->Data();
+          int d   = (int)(long)h->next->next->next->Data();
+          poly f; poly g;
+          henselFactors(hh, f0, g0, d, f, g);
+          lists L = (lists)omAllocBin(slists_bin);
+          L->Init(2);
+          L->m[0].rtyp = POLY_CMD; L->m[0].data=(void*)f;
+          L->m[1].rtyp = POLY_CMD; L->m[1].data=(void*)g;
+          res->rtyp = LIST_CMD;
+          res->data = (char *)L;
+          return FALSE;
+        }
+        else
+        {
+          Werror( "expected argument list (poly, poly, poly, int)");
+          return TRUE;
+        }
+      }
   /*==================== forking experiments ======================*/
       if(strcmp(sys_cmd, "waitforssilinks")==0)