Changeset f11ea16 in git for libpolys/polys

Timestamp:

Aug 14, 2011, 6:11:52 PM (13 years ago)

Author:

Oleksandr Motsak <motsak@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

2cae9214dc99498e560e328e02d738ab9e013739

Parents:

d8c7a64f3bcb71205b5fc350abe3a065f1b1576b

git-author:

Oleksandr Motsak <motsak@mathematik.uni-kl.de>2011-08-14 18:11:52+02:00

git-committer:

Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 16:12:19+01:00

Message:

FIX: duplicate of id_ChineseRemainder

File:

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

libpolys/polys/simpleideals.cc

@@ -1717 +1717 @@
 #endif
 
-/*2
-* xx,q: arrays of length 0..rl-1
-* xx[i]: SB mod q[i]
-* assume: char=0
-* assume: q[i]!=0
-* destroys xx
-*/
-#ifdef HAVE_FACTORY
-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;
-  poly r,h,hh,res_p;
-  number *x=(number *)omAlloc(rl*sizeof(number));
-  for(i=cnt-1;i>=0;i--)
-  {
-    res_p=NULL;
-    loop
-    {
-      r=NULL;
-      for(j=rl-1;j>=0;j--)
-      {
-        h=xx[j]->m[i];
-        if ((h!=NULL)
-        &&((r==NULL)||(p_LmCmp(r,h,R)==-1)))
-          r=h;
-      }
-      if (r==NULL) break;
-      h=p_Head(r, R);
-      for(j=rl-1;j>=0;j--)
-      {
-        hh=xx[j]->m[i];
-        if ((hh!=NULL) && (p_LmCmp(r,hh, R)==0))
-        {
-          x[j]=p_GetCoeff(hh, R);
-          hh=p_LmFreeAndNext(hh, R);
-          xx[j]->m[i]=hh;
-        }
-        else
-          x[j]=n_Init(0, R->cf); // is R->cf really n_Q???, yes!
-      }
-       
-      number n=nChineseRemainder(x,q,rl, R->cf);
-
-      for(j=rl-1;j>=0;j--)
-      {
-        x[j]=NULL; // nlInit(0...) takes no memory
-      }
-      if (n_IsZero(n, R->cf)) p_Delete(&h, R);
-      else
-      {
-        p_SetCoeff(h,n, R);
-        //Print("new mon:");pWrite(h);
-        res_p=p_Add_q(res_p, h, R);
-      }
-    }
-    result->m[i]=res_p;
-  }
-  omFree(x);
-  for(i=rl-1;i>=0;i--) id_Delete(&(xx[i]), R);
-  omFree(xx);
-  return result;
-}
-#endif
 
 /*2