Changeset cb4700 in git

Timestamp:

Jan 31, 2024, 4:59:43 PM (3 months ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

a410042175ed269bc5f88a7a1239e49c4c28b706

Parents:

687d0d05ffa27de4df1418fdc38be90ed25c8aac

Message:

opt: idSkipZeroes, add: id_Jet0

Location:

libpolys/polys

Files:

• simpleideals.cc (modified) (4 diffs)
• simpleideals.h (modified) (1 diff)

libpolys/polys/simpleideals.cc

@@ -216 +216 @@
       {
         ide->m[j] = ide->m[k];
+        ide->m[k] = NULL;
       }
     }
@@ -227 +228 @@
     if (j == -1)
       j = 0;
-    else
-    {
-      for (k=j+1; k<idelems; k++)
-        ide->m[k] = NULL;
-    }
     j++;
     pEnlargeSet(&(ide->m),idelems,j-idelems);
@@ -238 +234 @@
 }
 
-int idSkipZeroes0 (ideal ide)
+int idSkipZeroes0 (ideal ide) /*idSkipZeroes without realloc*/
 {
   assume (ide != NULL);
@@ -1786 +1782 @@
 }
 
+ideal id_Jet0(const ideal i, const ring R)
+{
+  ideal r=idInit((i->nrows)*(i->ncols),i->rank);
+  r->nrows = i-> nrows;
+  r->ncols = i-> ncols;
+  //r->rank = i-> rank;
+
+  for(long k=((long)(i->nrows))*((long)(i->ncols))-1;k>=0; k--)
+    r->m[k]=pp_Jet0(i->m[k],R);
+
+  return r;
+}
+
 ideal id_JetW(const ideal i,int d, intvec * iv, const ring R)
 {

libpolys/polys/simpleideals.h

@@ -139 +139 @@
 BOOLEAN id_IsZeroDim(ideal I, const ring r);
 ideal id_Jet(const ideal i,int d, const ring R);
+ideal id_Jet0(const ideal i, const ring R);/*id_Jet(i,0,R)*/
 ideal id_JetW(const ideal i,int d, intvec * iv, const ring R);
 ideal  id_Subst(ideal id, int n, poly e, const ring r);