Changeset 4cb74e in git

Timestamp:

May 5, 2023, 1:59:17 PM (12 months ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')

Children:

098b1e2c086aa89e01d5f7984f9c67c3708320e4

Parents:

d891c5fb7848eaccdb188900d38e88c00cfaae1a

Message:

2 algorithm for Hilbert series, check degree

File:

• kernel/combinatorics/hilb.cc (modified) (1 diff)

kernel/combinatorics/hilb.cc

@@ -2050 +2050 @@
 intvec* hFirstSeries(ideal A,intvec *module_w,ideal Q, intvec *wdegree)
 {
-  void (*WerrorS_save)(const char *s) = WerrorS_callback;
-  WerrorS_callback=WerrorS_dummy;
-  intvec* res=hFirstSeries1(A,module_w,Q,wdegree);
-  WerrorS_callback=WerrorS_save;
-  if (errorreported==0)
-  {
-    return res;
-  }
-  else errorreported=0;// retry with other alg.:
+  // find degree bound
+  int a,b,prod;
+  a=rVar(currRing);
+  b=1;
+  prod=a;
+  while(prod<(1<<15) && (a>1))
+  {
+    a--;b++;
+    prod*=a;
+    prod/=b;
+  }
+  if (a==1) b=(1<<15);
+  // check degree bound
+  BOOLEAN large_deg=FALSE;
+  int max=0;
+  for(int i=IDELEMS(A)-1;i>=0;i--)
+  {
+    if (A->m[i]!=NULL)
+    {
+      int mm=p_Totaldegree(A->m[i],currRing);
+      if (mm>max)
+      {
+        max=mm;
+        if (max>=b)
+        {
+          large_deg=TRUE;
+          break;
+        }
+      }
+    }
+  }
+  intvec* res;
+  if (!large_deg)
+  {
+    void (*WerrorS_save)(const char *s) = WerrorS_callback;
+    WerrorS_callback=WerrorS_dummy;
+    res=hFirstSeries1(A,module_w,Q,wdegree);
+    WerrorS_callback=WerrorS_save;
+    if (errorreported==0)
+    {
+      return res;
+    }
+    else errorreported=0;// retry with other alg.:
+  }
 
   static ring hilb_Qt=NULL;