Changeset ad25ace in git

Timestamp:

Apr 13, 2000, 3:01:49 PM (23 years ago)

Author:

Hans Schönemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '0604212ebb110535022efecad887940825b97c3f')

Children:

8e153e5516cee595e5037ab243fc9a0d54dad39e

Parents:

4fc736efe01670e3ed92d3911fa633d760f79b7c

Message:

*hannes: bug 19 (vector/poly)


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

File:

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

Singular/iparith.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: iparith.cc,v 1.208 2000-04-11 15:40:22 Singular Exp $ */
+/* $Id: iparith.cc,v 1.209 2000-04-13 13:01:49 Singular Exp $ */
 
 /*
@@ -117 +117 @@
 extern BOOLEAN expected_parms;
 
+#define ii_div_by_0 "div. by 0"
 int iiOp; /* the current operation*/
 
@@ -1018 +1019 @@
   if (b==0)
   {
-    WerrorS("div. by 0");
+    WerrorS(ii_div_by_0);
     return TRUE;
   }
@@ -1044 +1045 @@
   if (nIsZero(q))
   {
-    WerrorS("div. by 0");
+    WerrorS(ii_div_by_0);
     return TRUE;
   }
@@ -1057 +1058 @@
   if (q==NULL)
   {
-    WerrorS("div. by 0");
+    WerrorS(ii_div_by_0);
     return TRUE;
   }
@@ -1063 +1064 @@
   {
 #ifdef HAVE_FACTORY
-    res->data=(void*)(singclap_pdivide((poly)(u->Data()),
-                                       (poly)(v->Data())));
+    poly p=(poly)(u->Data());
+    if(pGetComp(p)==0)
+    {
+      res->data=(void*)(singclap_pdivide(p /*(poly)(u->Data())*/ ,
+                                         q /*(poly)(v->Data())*/ ));
+    }
+    else
+    {
+      int comps=pMaxComp(p);
+      ideal I=idInit(comps,1);
+      p=(poly)(u->CopyD());
+      poly h;
+      int i;
+      // conversio to a list of polys:
+      while (p!=NULL)
+      {
+        i=pGetComp(p)-1;
+        h=pNext(p);
+        pNext(p)=NULL;
+        pSetComp(p,0);
+        I->m[i]=pAdd(I->m[i],p);
+        p=h;
+      }
+      // division and conversion to vector:
+      h=NULL;
+      p=NULL;
+      for(i=comps-1;i>=0;i--)
+      {
+        if (I->m[i]!=NULL)
+        {
+          h=singclap_pdivide(I->m[i],q);
+          pSetCompP(h,i+1);
+          p=pAdd(p,h);
+        }
+      }
+      idDelete(&I);
+      res->data=(void *)p;
+    }
 #else
     WerrorS("division only by a monomial");