Esempio n. 1
0
 public final MVPerm compose(MVPerm perm) {
   MVPerm res = new MVPerm();
   for (int i = 0; i < this.elems.length; i++) {
     ModelValue mv = this.elems[i];
     if (mv == null) {
       res.put(i, perm.elems[i]);
     } else {
       ModelValue mv1 = perm.elems[mv.index];
       if (mv1 == null) {
         res.put(i, mv);
       } else if (!ModelValue.mvs[i].equals(mv1)) {
         res.put(i, mv1);
       }
     }
   }
   return res;
 }
Esempio n. 2
0
 public static final MVPerm[] permutationSubgroup(ValueEnumeration varEnum) {
   Set perms = new Set(20);
   Vect permVec = new Vect(20);
   // Compute the group generators:
   Value elem;
   while ((elem = varEnum.nextElement()) != null) {
     FcnRcdValue fcn = FcnRcdValue.convert(elem);
     if (fcn == null) {
       Assert.fail("The symmetry operator must specify a set of functions.");
     }
     MVPerm perm = new MVPerm();
     for (int i = 0; i < fcn.domain.length; i++) {
       Value dval = fcn.domain[i];
       Value rval = fcn.values[i];
       if ((dval instanceof ModelValue) && (rval instanceof ModelValue)) {
         perm.put((ModelValue) dval, (ModelValue) rval);
       } else {
         Assert.fail("Symmetry function must have model values as domain and range.");
       }
     }
     if (perm.size() > 0 && perms.put(perm) == null) {
       permVec.addElement(perm);
     }
   }
   // Compute the group generated by the generators:
   int gsz = permVec.size();
   int sz0 = 0;
   while (true) {
     int sz1 = permVec.size();
     for (int i = 0; i < gsz; i++) {
       MVPerm perm1 = (MVPerm) permVec.elementAt(i);
       for (int j = sz0; j < sz1; j++) {
         MVPerm perm = perm1.compose((MVPerm) permVec.elementAt(j));
         if (perm.size() > 0 && perms.put(perm) == null) {
           permVec.addElement(perm);
         }
       }
     }
     if (sz1 == permVec.size()) break;
     sz0 = sz1;
   }
   // Finally, put all the elements in an array ready for use:
   MVPerm[] res = new MVPerm[permVec.size()];
   Enumeration permEnum = permVec.elements();
   for (int i = 0; i < res.length; i++) {
     res[i] = (MVPerm) permEnum.nextElement();
   }
   return res;
 }