Changeset 8fe22c in git

Timestamp:

Feb 25, 2002, 7:13:45 PM (21 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')

Children:

962b1f0a4db237d3ca3ca0a559c3f5d44649ebaf

Parents:

2ca72f72ba5811479c336c3a36c9d120db4edfc7

Message:

*mschulze: added jjDIVISION4


git-svn-id: file:///usr/local/Singular/svn/trunk@5924 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• Singular/iparith.cc (modified) (2 diffs)

Singular/iparith.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: iparith.cc,v 1.274 2002-02-06 08:52:21 Singular Exp $ */
+/* $Id: iparith.cc,v 1.275 2002-02-25 18:13:45 mschulze Exp $ */
 
 /*
@@ -4814 +4814 @@
   return FALSE;
 }
+
+static BOOLEAN jjDIVISION4(leftv res, leftv v)
+{
+  int typ=v->Typ();
+  if(v->Typ()!=MODUL_CMD&&v->Typ()!=IDEAL_CMD)
+  {
+    Warn("1st argument must be an ideal or module!");
+    return TRUE;
+  }
+  ideal P=(ideal)v->Data();
+  v=v->next;
+  if(v->Typ()!=typ)
+  {
+    Warn("2nd argument must be an ideal or module!");
+    return TRUE;
+  }
+  assumeStdFlag(v);
+  ideal Q=(ideal)v->Data();
+  v=v->next;
+  if(v->Typ()!=INT_CMD)
+  {
+    Warn("3rd argument must be an int!");
+    return TRUE;
+  }
+  int n=(int)v->Data();
+  v=v->next;
+  if(v->Typ()!=INTVEC_CMD)
+  {
+    Warn("4th argument must be an intvec!");
+    return TRUE;
+  }
+  short *w=iv2array((intvec *)v->Data());
+
+  int N=0;
+  for(int i=IDELEMS(Q)-1;i>=0;i--)
+    N=max(N,pDegW(Q->m[i],w));
+  N+=n;
+
+  matrix T=mpNew(IDELEMS(Q),IDELEMS(P));
+  ideal R=idInit(IDELEMS(P),P->rank);
+  matrix U=mpNew(P->rank,P->rank);
+  for(int i=MATROWS(U);i>=1;i--)
+    MATELEM(U,i,i)=pOne();
+
+  for(int i=IDELEMS(P)-1;i>=0;i--)
+  {
+    poly p=ppJetW(P->m[i],N,w);
+
+    int j=IDELEMS(Q)-1;
+    while(p!=NULL)
+    {
+      if(pDivisibleBy(Q->m[j],p))
+      {
+        poly p0=pDivideM(pHead(p),pHead(Q->m[j]));
+        if(pDegW(p0,w)>n)
+          pDelete(&p0);
+        else
+          MATELEM(T,j+1,i+1)=pAdd(MATELEM(T,j+1,i+1),p0);
+        p0=ksOldSpolyRed(Q->m[j],p,0);
+        p=ppJetW(p0,N,w);
+        pDelete(&p0);
+        j=IDELEMS(Q)-1;
+      }
+      else
+      {
+        if(j==0)
+        {
+          poly p0=p;
+          pIter(p);
+          pNext(p0)=NULL;
+          if(pDegW(p0,w)>n)
+            pDelete(&p0);
+          else
+            R->m[i]=pAdd(R->m[i],p0);
+          j=IDELEMS(Q)-1;
+        }
+        else
+          j--;
+      }
+    }
+  }
+
+  omFree(w);
+
+  lists L=(lists)omAllocBin(slists_bin);
+  L->Init(3);
+  L->m[0].rtyp=MATRIX_CMD;
+  L->m[0].data=(void *)T;
+  L->m[1].rtyp=MATRIX_CMD;
+  L->m[1].data=(void *)U;
+  L->m[2].rtyp=typ;
+  L->m[2].data=(void *)R;
+  res->data=L;
+  res->rtyp=LIST_CMD;
+
+  return FALSE;
+}
+
 static BOOLEAN jjIDEAL_PL(leftv res, leftv v)
 {