Changeset 40df4d in git

Timestamp:

Dec 17, 2013, 1:44:17 PM (10 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

035192fb77c255d934b21cd0a32eb12359bda176

Parents:

a1e807a999ad81c88502c684ee93531e7840ea82

Message:

fix: new grobcov.lib (a.montes)

File:

• Singular/LIB/grobcov.lib (modified) (6 diffs)

Singular/LIB/grobcov.lib

@@ -113 +113 @@
                    Can be used without calling presiouly setglobalrings(),
 
-locus(G);     Special routine for determining the locus of points
+locus(G):     Special routine for determining the locus of points
                    of  objects. Given a parametric ideal J with
                    parameters (a_1,..a_m) and variables (x_1,..,xn),
@@ -131 +131 @@
                    Can be used without calling presiouly setglobalrings(),
 
-locusdg(G);  Is a special routine for computing the locus in dinamical geometry.
+locusdg(G):  Is a special routine for computing the locus in dinamical geometry.
                     It detects automatically a possible point that is to be avoided by the mover,
                     whose coordinates must be the last two coordinates in the definition of the ring.
@@ -139 +139 @@
                     Can be used without calling presiouly setglobalrings(),
 
-locusto(L);  Transforms the output of locus to a string that
+locusto(L):  Transforms the output of locus to a string that
                   can be reed from different computational systems.
                   Can be used without calling presiouly setglobalrings(),
 
-addcons(L); Given a disjoint set of locally closed subsets in P-representation,
+addcons(L): Given a disjoint set of locally closed subsets in P-representation,
                   it returns the canonical P-representation of the constructible set
                   formed by the union of them.
@@ -1222 +1222 @@
 proc PtoCrep(list L)
 {
+  int te=0;
   def RR=basering;
+  if(defined(@P)){te=1;}
+  else{te=0; setglobalrings();}
   setring(@P);
   def Lp=imap(RR,L);
@@ -1240 +1243 @@
   def La=list(ida,idb);
   setring(RR);
-  return(imap(@P,La));
+  def LL=imap(@P,La);
+  if(te==0){kill @P; kill @R; kill @RP;}
+  return(LL);
 }
 
@@ -3895 +3900 @@
   int i; int j; int k; int l;
   intvec v; intvec v1; intvec v0;
-  ideal J; list K;
+  ideal J; list K; list L1;
+  for(i=1;i<=size(L);i++)
+  {
+    if(equalideals(L[i][1][1],ideal(1))==0)
+    {L1[size(L1)+1]=L[i];}
+  }
+  L=L1;
   for(i=1;i<=size(L);i++)
   {