Changeset aef5db in git for libpolys/polys/matpol.cc

Timestamp:

Jun 29, 2018, 2:13:21 PM (6 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '6e5adcba05493683b94648c659a729c189812c77')

Children:

96cd73ba8e350c13e7964b70073c3fec90b296c4

Parents:

13673122266edb13d76ff215c7d27bf6d7aa7ce3

Message:

add: system("unflatten",<smatrix>,<cols>)

File:

• libpolys/polys/matpol.cc (modified) (1 diff)

libpolys/polys/matpol.cc

@@ -1843 +1843 @@
   return res;
 }
+
+ideal sm_UnFlatten(ideal a, int col, const ring R)
+{
+  if ((IDELEMS(a)!=1)
+  ||((a->rank % col)!=0))
+  {
+    Werror("wrong format: %d x %d for unflatten",(int)a->rank,IDELEMS(a));
+    return NULL;
+  }
+  int row=a->rank/col;
+  ideal res=idInit(col,row);
+  poly p=a->m[0];
+  while(p!=NULL)
+  {
+    poly h=p_Head(p,R);
+    int comp=p_GetComp(h,R);
+    int c=(comp-1)/row;
+    int r=comp%row; if (r==0) r=row;
+    p_SetComp(h,r,R); p_SetmComp(h,R);
+    res->m[c]=p_Add_q(res->m[c],h,R);
+    pIter(p);
+  }
+  return res;
+}
+
 /*2
 *returns the trace of matrix a