Changeset 1e44f4e in git for Singular/LIB/gitfan.lib

Timestamp:

Jun 8, 2017, 5:11:39 PM (7 years ago)

Author:

Janko Boehm <boehm@…>

Branches:

(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')

Children:

ca40dc8a6a76a26cf2c696db58f6554dc49bb33d

Parents:

aafc931bc91661f27f2b47edfb2d77d93d3df9d0

git-author:

Janko Boehm <boehm@mathematik.uni-kl.de>2017-06-08 17:11:39+02:00

git-committer:

Yue Ren <ren@mathematik.uni-kl.de>2017-06-08 18:10:27+02:00

Message:

Update gitfan.lib

File:

• Singular/LIB/gitfan.lib (modified) (7 diffs)

Singular/LIB/gitfan.lib

@@ -33 +33 @@
   permutation; Permutation in map representation.
 
-PROCEDURES (not using group action):
+PROCEDURES:
   isAface(ideal, intvec); Checks whether the given face is an a-face.
   afaces(ideal); Returns a list of intvecs that correspond to the set of all a-faces, optionally for given list of simplex faces.
@@ -45 +45 @@
   GKZfan(intmat); Returns the GKZ-fan of the matrix Q.
   movingCone(intmat); Compute the moving cone.
-
-PROCEDURES (using group action):
   computeAfaceOrbits(list, list); Compute orbits of a-faces under a permutation group action.
   minimalAfaceOrbits(list); Find the minimal a-face orbits.
@@ -362 +360 @@
 
   isAface(I,w); // should be 0
-  "-----------";
   isAface(I,v); // should be 1
 }
@@ -1051 +1048 @@
   intmat Qt = transpose(Q);
   for(int i = 1; i <= size(G); i++){
-    i;
     intvec sig = permutationToIntvec(G[i]);
     intmat Bsig = perm2mat(sig, indices,Qt);
@@ -1186 +1182 @@
                      permutationFromIntvec(intvec( 5, 7, 1, 6, 9, 2, 8, 4, 10, 3 ));
   list Asigmagens = groupActionOnQImage(generatorsG,Q);
-groupActionOnHashes(Asigmagens,OC);
-}
+groupActionOnHashes(Asigmagens,list(OC));
+
+
+  list orb = findOrbits(simplexSymmetryGroup,nrows(Q));
+  list simplexOrbitRepresentatives;
+  for (int i=1;i<=size(orb);i++){simplexOrbitRepresentatives[i]=orb[i][1];}
+  list afaceOrbitRepresentatives=afaces(J,simplexOrbitRepresentatives);
+  list fulldimAfaceOrbitRepresentatives=fullDimImages(afaceOrbitRepresentatives,Q);
+  list afaceOrbits=computeAfaceOrbits(fulldimAfaceOrbitRepresentatives,simplexSymmetryGroup);
+  list minAfaceOrbits = minimalAfaceOrbits(afaceOrbits);
+  list listOfOrbitConeOrbits = orbitConeOrbits(minAfaceOrbits,Q);
+  list listOfMinimalOrbitConeOrbits = minimalOrbitConeOrbits(listOfOrbitConeOrbits);
+  list Asigma = groupActionOnQImage(simplexSymmetryGroup,Q);
+  groupActionOnHashes(Asigma,listOfOrbitConeOrbits);
+
+}
+
+
 
 
@@ -1977 +1989 @@
 
 
-proc subdividelist(list OC,int ncores){
+static proc subdividelist(list OC,int ncores){
   if (ncores>size(OC)){ncores=size(OC);}
   int percore = size(OC) div (ncores);
@@ -2472 +2484 @@
   AF=fullDimImages(AF,Q);
   AF = minimalAfaces(AF);
-  list OC = orbitCones(AF,Q);
+  OC = orbitCones(AF,Q);
   cone Qgamma = coneViaPoints(transpose(Q));
-  list GIT = GITfanParallel(OC,Q,Qgamma);
+  GIT = GITfanParallel(OC,Q,Qgamma);
   }
   fan Sigma = hashesToFan(GIT,OC);