Changeset de27d8 in git for libpolys/polys

Timestamp:

Feb 18, 2013, 7:27:17 PM (11 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')

Children:

7d1ce3ef53e17b72b36010a52fcfbdbde745a3e7afc974f80049f26666102bf7d90f3aa6c1dfd957

Parents:

77ed18101c07cdbb767d8703e5ea1d0f7c8b05f6

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2013-02-18 19:27:17+01:00

git-committer:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2013-02-18 19:27:55+01:00

Message:

fix: intdiv/intmod/chinrem copied from master

Location:

libpolys/polys

Files:

• monomials/p_polys.cc (modified) (1 diff)
• simpleideals.cc (modified) (1 diff)
• simpleideals.h (modified) (1 diff)

libpolys/polys/monomials/p_polys.cc

@@ -128 +128 @@
         x[j]=n_Init(0, R);
     }
-    number n=n_ChineseRemainder(x,q,rl,R->cf);
+    number n=n_ChineseRemainderSym(x,q,rl,TRUE,R->cf);
     for(j=rl-1;j>=0;j--)
     {

libpolys/polys/simpleideals.cc

@@ -1768 +1768 @@
   return(idResult);
 }
+
+ideal id_ChineseRemainder(ideal *xx, number *q, int rl, const ring r)
+{
+  int cnt=IDELEMS(xx[0])*xx[0]->nrows;
+  ideal result=idInit(cnt,xx[0]->rank);
+  result->nrows=xx[0]->nrows; // for lifting matrices
+  result->ncols=xx[0]->ncols; // for lifting matrices
+  int i,j;
+  number *x=(number *)omAlloc(rl*sizeof(number));
+  poly *p=(poly *)omAlloc(rl*sizeof(poly));
+  for(i=cnt-1;i>=0;i--)
+  {
+    for(j=rl-1;j>=0;j--)
+    {
+      p[j]=xx[j]->m[i];
+    }
+    result->m[i]=p_ChineseRemainder(p,x,q,rl,r);
+    for(j=rl-1;j>=0;j--)
+    {
+      xx[j]->m[i]=p[j];
+    }
+  }
+  omFreeSize(p,rl*sizeof(poly));
+  omFreeSize(x,rl*sizeof(number));
+  for(i=rl-1;i>=0;i--) id_Delete(&(xx[i]),r);
+  omFreeSize(xx,rl*sizeof(ideal));
+  return result;
+}

libpolys/polys/simpleideals.h

@@ -147 +147 @@
 
 
+ideal id_ChineseRemainder(ideal *xx, number *q, int rl, const ring r);
 #endif