Changeset a5bbb5 in git for Singular/iparith.cc

Timestamp:

Jan 9, 2017, 2:02:30 PM (7 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

b8b76ef3092b3162aa763b4c5b7705ceab2203fa

Parents:

a418d9fffbccc54b2388850e64bfc084cd0e61dcd29e550951bb89f260d4afb0da084c6a1a83e101

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2017-01-09 14:02:30+01:00

git-committer:

GitHub <noreply@github.com>2017-01-09 14:02:30+01:00

Message:

Merge pull request #810 from YueRen/tropicallib

chg: kghostview -> xdg-open + line breaks after commata in tex file

File:

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

Singular/iparith.cc

@@ -24 +24 @@
 #include <misc/prime.h>
 
-#include <polys/prCopy.h>
 #include <polys/matpol.h>
 #include <polys/monomials/maps.h>
@@ -1871 +1870 @@
   if (rField_is_Ring(currRing))
   {
-    //ring origR = currRing;
-    //ring tempR = rCopy(origR);
-    //coeffs new_cf=nInitChar(n_Q,NULL);
-    //nKillChar(tempR->cf);
-    //tempR->cf=new_cf;
-    //rComplete(tempR);
     ideal vid = (ideal)v->Data();
     int i = idPosConstant(vid);
@@ -1884 +1877 @@
       return FALSE;
     }
-    //rChangeCurrRing(tempR);
-    //ideal vv = idrCopyR(vid, origR, currRing);
     ideal vv = id_Copy(vid, currRing);
-    //ideal ww = idrCopyR((ideal)w->Data(), origR, currRing);
     ideal ww = id_Copy((ideal)w->Data(), currRing);
     /* drop degree zero generator from vv (if any) */
@@ -1895 +1885 @@
     res->data = (char *)d;
     idDelete(&vv); idDelete(&ww);
-    //rChangeCurrRing(origR);
-    //rDelete(tempR);
     return FALSE;
   }
@@ -2309 +2297 @@
   if (rField_is_Ring_Z(currRing))
   {
-    ring origR = currRing;
-    ring tempR = rCopy(origR);
-    coeffs new_cf=nInitChar(n_Q,NULL);
-    nKillChar(tempR->cf);
-    tempR->cf=new_cf;
-    rComplete(tempR);
-    ideal uid = (ideal)u->Data();
-    rChangeCurrRing(tempR);
-    ideal uu = idrCopyR(uid, origR, currRing);
-    sleftv uuAsLeftv; memset(&uuAsLeftv, 0, sizeof(uuAsLeftv));
-    uuAsLeftv.rtyp = IDEAL_CMD;
-    uuAsLeftv.data = uu; uuAsLeftv.next = NULL;
-    if (hasFlag(u, FLAG_STD)) setFlag(&uuAsLeftv,FLAG_STD);
-    assumeStdFlag(&uuAsLeftv);
     PrintS("// NOTE: computation of Hilbert series etc. is being\n");
     PrintS("//       performed for generic fibre, that is, over Q\n");
-    intvec *module_w=(intvec*)atGet(&uuAsLeftv,"isHomog",INTVEC_CMD);
-    intvec *iv=hFirstSeries(uu,module_w,currRing->qideal);
-    int returnWithTrue = 1;
-    switch((int)(long)v->Data())
-    {
-      case 1:
-        res->data=(void *)iv;
-        returnWithTrue = 0;
-      case 2:
-        res->data=(void *)hSecondSeries(iv);
-        delete iv;
-        returnWithTrue = 0;
-    }
-    if (returnWithTrue)
-    {
-      WerrorS(feNotImplemented);
-      delete iv;
-    }
-    idDelete(&uu);
-    rChangeCurrRing(origR);
-    rDelete(tempR);
-    if (returnWithTrue) return TRUE; else return FALSE;
   }
 #endif
@@ -3768 +3720 @@
 {
   SPrintStart();
+#ifdef HAVE_RINGS
   if (rField_is_Ring_Z(currRing))
   {
-    ring origR = currRing;
-    ring tempR = rCopy(origR);
-    coeffs new_cf=nInitChar(n_Q,NULL);
-    nKillChar(tempR->cf);
-    tempR->cf=new_cf;
-    rComplete(tempR);
-    ideal vid = (ideal)v->Data();
-    rChangeCurrRing(tempR);
-    ideal vv = idrCopyR(vid, origR, currRing);
-    sleftv vvAsLeftv; memset(&vvAsLeftv, 0, sizeof(vvAsLeftv));
-    vvAsLeftv.rtyp = IDEAL_CMD;
-    vvAsLeftv.data = vv; vvAsLeftv.next = NULL;
-    if (hasFlag(v, FLAG_STD)) setFlag(&vvAsLeftv,FLAG_STD);
-    assumeStdFlag(&vvAsLeftv);
     PrintS("// NOTE: computation of degree is being performed for\n");
     PrintS("//       generic fibre, that is, over Q\n");
-    intvec *module_w=(intvec*)atGet(&vvAsLeftv,"isHomog",INTVEC_CMD);
-    scDegree(vv,module_w,currRing->qideal);
-    idDelete(&vv);
-    rChangeCurrRing(origR);
-    rDelete(tempR);
-  }
+  }
+#endif
   assumeStdFlag(v);
   intvec *module_w=(intvec*)atGet(v,"isHomog",INTVEC_CMD);
@@ -3918 +3853 @@
      Warn("dim(%s) may be wrong because the mixed monomial ordering",v->Name());
   }
+#ifdef HAVE_RINGS
   if (rField_is_Ring(currRing))
   {
@@ -3982 +3918 @@
     return FALSE;
   }
+#endif
   res->data = (char *)(long)scDimInt((ideal)(v->Data()),currRing->qideal);
   return FALSE;
@@ -4146 +4083 @@
 static BOOLEAN jjHILBERT(leftv, leftv v)
 {
+#ifdef HAVE_RINGS
   if (rField_is_Ring_Z(currRing))
   {
-    ring origR = currRing;
-    ring tempR = rCopy(origR);
-    coeffs new_cf=nInitChar(n_Q,NULL);
-    nKillChar(tempR->cf);
-    tempR->cf=new_cf;
-    rComplete(tempR);
-    ideal vid = (ideal)v->Data();
-    rChangeCurrRing(tempR);
-    ideal vv = idrCopyR(vid, origR, currRing);
-    sleftv vvAsLeftv; memset(&vvAsLeftv, 0, sizeof(vvAsLeftv));
-    vvAsLeftv.rtyp = IDEAL_CMD;
-    vvAsLeftv.data = vv; vvAsLeftv.next = NULL;
-    if (hasFlag(v, FLAG_STD)) setFlag(&vvAsLeftv,FLAG_STD);
-    assumeStdFlag(&vvAsLeftv);
     PrintS("// NOTE: computation of Hilbert series etc. is being\n");
     PrintS("//       performed for generic fibre, that is, over Q\n");
-    intvec *module_w=(intvec*)atGet(&vvAsLeftv,"isHomog",INTVEC_CMD);
-    //scHilbertPoly(vv,currRing->qideal);
-    hLookSeries(vv,module_w,currRing->qideal);
-    idDelete(&vv);
-    rChangeCurrRing(origR);
-    rDelete(tempR);
-    return FALSE;
-  }
+  }
+#endif
   assumeStdFlag(v);
   intvec *module_w=(intvec*)atGet(v,"isHomog",INTVEC_CMD);
@@ -4180 +4098 @@
 static BOOLEAN jjHILBERT_IV(leftv res, leftv v)
 {
+#ifdef HAVE_RINGS
   if (rField_is_Ring_Z(currRing))
   {
@@ -4185 +4104 @@
     PrintS("//       performed for generic fibre, that is, over Q\n");
   }
+#endif
   res->data=(void *)hSecondSeries((intvec *)v->Data());
   return FALSE;
@@ -4293 +4213 @@
 {
   ideal result=kInterRed((ideal)(v->Data()), currRing->qideal);
+#ifdef HAVE_RINGS
   if(rField_is_Ring(currRing))
-    Warn("interred: this command is experimental over the integers");
+    WarnS("interred: this command is experimental over the integers");
+#endif
   if (TEST_OPT_PROT) { PrintLn(); mflush(); }
   res->data = result;
@@ -5782 +5704 @@
     return TRUE;
   }
+#ifdef HAVE_RINGS
   if (rField_is_Ring_Z(currRing))
   {
-    ring origR = currRing;
-    ring tempR = rCopy(origR);
-    coeffs new_cf=nInitChar(n_Q,NULL);
-    nKillChar(tempR->cf);
-    tempR->cf=new_cf;
-    rComplete(tempR);
-    ideal uid = (ideal)u->Data();
-    rChangeCurrRing(tempR);
-    ideal uu = idrCopyR(uid, origR, currRing);
-    sleftv uuAsLeftv; memset(&uuAsLeftv, 0, sizeof(uuAsLeftv));
-    uuAsLeftv.rtyp = IDEAL_CMD;
-    uuAsLeftv.data = uu; uuAsLeftv.next = NULL;
-    if (hasFlag(u, FLAG_STD)) setFlag(&uuAsLeftv,FLAG_STD);
-    assumeStdFlag(&uuAsLeftv);
     PrintS("// NOTE: computation of Hilbert series etc. is being\n");
     PrintS("//       performed for generic fibre, that is, over Q\n");
-    intvec *module_w=(intvec*)atGet(&uuAsLeftv,"isHomog",INTVEC_CMD);
-    intvec *iv=hFirstSeries(uu,module_w,currRing->qideal,wdegree);
-    int returnWithTrue = 1;
-    switch((int)(long)v->Data())
-    {
-      case 1:
-        res->data=(void *)iv;
-        returnWithTrue = 0;
-      case 2:
-        res->data=(void *)hSecondSeries(iv);
-        delete iv;
-        returnWithTrue = 0;
-    }
-    if (returnWithTrue)
-    {
-      WerrorS(feNotImplemented);
-      delete iv;
-    }
-    idDelete(&uu);
-    rChangeCurrRing(origR);
-    rDelete(tempR);
-    if (returnWithTrue) return TRUE; else return FALSE;
-  }
+  }
+#endif
   assumeStdFlag(u);
   intvec *module_w=(intvec *)atGet(u,"isHomog",INTVEC_CMD);
@@ -5836 +5724 @@
       return FALSE;
   }
+  delete iv;
   WerrorS(feNotImplemented);
-  delete iv;
   return TRUE;
 }
@@ -5923 +5811 @@
                                (matrix)v->CopyD());
   return FALSE;
-}
-static BOOLEAN currRingIsOverIntegralDomain ()
-{
-  /* true for fields and Z, false otherwise */
-  if (rField_is_Ring_PtoM(currRing)) return FALSE;
-  if (rField_is_Ring_2toM(currRing)) return FALSE;
-  if (rField_is_Ring_ModN(currRing)) return FALSE;
-  return TRUE;
 }
 static BOOLEAN jjMINOR_M(leftv res, leftv v)
@@ -6094 +5974 @@
   }
   if ((!noAlgorithm) && (strcmp(algorithm, "Bareiss") == 0)
-      && (!currRingIsOverIntegralDomain()))
+      && (!rField_is_Domain(currRing)))
   {
     Werror("Bareiss algorithm not defined over coefficient rings %s",
@@ -9161 +9041 @@
   }
   #endif
+#ifdef HAVE_RINGS
   if (rField_is_Ring(currRing))
   {
@@ -9181 +9062 @@
     }
   }
+#endif
   return FALSE;
 }