[BACK]Return to minimal.k CVS log [TXT][DIR] Up to [local] / OpenXM / src / k097 / lib / minimal

Diff for /OpenXM/src/k097/lib/minimal/minimal.k between version 1.6 and 1.9

version 1.6, 2000/05/06 07:58:37 version 1.9, 2000/05/06 13:41:12
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/k097/lib/minimal/minimal.k,v 1.5 2000/05/05 08:13:49 takayama Exp $ */  /* $OpenXM: OpenXM/src/k097/lib/minimal/minimal.k,v 1.8 2000/05/06 10:45:43 takayama Exp $ */
 #define DEBUG 1  #define DEBUG 1
 /* #define ORDINARY 1 */  /* #define ORDINARY 1 */
 /* If you run this program on openxm version 1.1.2 (FreeBSD),  /* If you run this program on openxm version 1.1.2 (FreeBSD),
Line 37  def load_tower() {
Line 37  def load_tower() {
     sm1(" [(parse) (k0-tower.sm1) pushfile ] extension ");      sm1(" [(parse) (k0-tower.sm1) pushfile ] extension ");
     sm1(" /k0-tower.sm1.loaded 1 def ");      sm1(" /k0-tower.sm1.loaded 1 def ");
   }    }
     sm1(" oxNoX ");
 }  }
 load_tower();  load_tower();
 SonAutoReduce = true;  SonAutoReduce = true;
Line 1152  def Sschreyer(g) {
Line 1153  def Sschreyer(g) {
                   /* i must be equal to f[2], I think. Double check. */                    /* i must be equal to f[2], I think. Double check. */
   
                   /* Correction Of Constant */                    /* Correction Of Constant */
                   c2 = -f[6];  /* or f[6]?  Double check. */                    /* Correction of syzygy */
                     c2 = f[6];  /* or -f[6]?  Double check. */
                     Print("c2="); Println(c2);
                   nn = Length(bases);                    nn = Length(bases);
                   for (ii=0; ii<nn;ii++) {                    for (ii=0; ii<nn;ii++) {
                      if ((ii != place) && (! IsNull(bases[ii]))) {                       if ((ii != i) && (! IsNull(bases[ii]))) {
                        bases[ii] = bases[ii]*c2;                         m = Length(bases[ii]);
                          for (jj=0; jj<m; jj++) {
                            if (jj != place) {
                              bases[ii,jj] = bases[ii,jj]*c2;
                            }
                          }
                      }                       }
                   }                    }
   
                     Print("Old freeRes[level] = "); sm1_pmat(freeRes[level]);
                   freeRes[level] = bases;                    freeRes[level] = bases;
                     Print("New freeRes[level] = "); sm1_pmat(freeRes[level]);
   
                  /* Update the freeRes[level-1] */                   /* Update the freeRes[level-1] */
                     Print("Old freeRes[level-1] = "); sm1_pmat(freeRes[level-1]);
                   bases = freeRes[level-1];                    bases = freeRes[level-1];
                   bases[place] = f[0];                    bases[place] = f[0];
                   freeRes[level-1] = bases;                    freeRes[level-1] = bases;
                     Print("New freeRes[level-1] = "); sm1_pmat(freeRes[level-1]);
   
                   reducer[level-1,place] = f[1];                    reducer[level-1,place] = f[1]-SunitOfFormat(place,f[1]);
                      /* This reducer is different from that of SlaScala(). */
   
                     reducerBasis = reducer[level-1];
                     nn = Length(reducerBasis);
                     for (ii=0; ii<nn;ii++) {
                        if ((ii != place) && (! IsNull(reducerBasis[ii]))) {
                          m = Length(reducerBasis[ii]);
                          for (jj=0; jj<m; jj++) {
                            if (jj != place) {
                              reducerBasis[ii,jj] = reducerBasis[ii,jj]*c2;
                            }
                          }
                        }
                     }
                     reducer[level-1] = reducerBasis;
   
                }else{                 }else{
                   /* redundantTable[level,i] = 0; */                    /* redundantTable[level,i] = 0; */
                   bases = freeRes[level];                    bases = freeRes[level];
Line 1182  def Sschreyer(g) {
Line 1210  def Sschreyer(g) {
     if (level >= 1) {      if (level >= 1) {
       Println(" ");        Println(" ");
       Print("Triangulating reducer at level "); Println(level-1);        Print("Triangulating reducer at level "); Println(level-1);
         Println("freeRes[level]="); sm1_pmat(freeRes[level]);
       reducerBase = reducer[level-1];        reducerBase = reducer[level-1];
       Print("reducerBase=");  Println(reducerBase);        Print("reducerBase=");  Println(reducerBase);
         Println("Compare freeRes[level] and reducerBase (put -1)");
       m = Length(reducerBase);        m = Length(reducerBase);
       for (ii=m-1; ii>=0; ii--) {        for (ii=m-1; ii>=0; ii--) {
         if (!IsNull(reducerBase[ii])) {          if (!IsNull(reducerBase[ii])) {
            for (jj=ii-1; jj>=0; jj--) {             for (jj=ii-1; jj>=0; jj--) {
              if (!IsNull(reducerBase[jj])) {               if (!IsNull(reducerBase[jj])) {
               if (!IsZero(reducerBase[jj,ii])) {                if (!IsZero(reducerBase[jj,ii])) {
                 reducerBase[jj] = reducerBase[jj]-reducerBase[jj,ii]*reducerBase[ii];                  /* reducerBase[ii,ii] should be always constant. */
                   reducerBase[jj] = reducerBase[ii,ii]*reducerBase[jj]-reducerBase[jj,ii]*reducerBase[ii];
               }                }
              }               }
            }             }
Line 1337  def Sminimal_v(g) {
Line 1368  def Sminimal_v(g) {
           if (level < maxLevel-1) {            if (level < maxLevel-1) {
             bases = freeRes[level+1];              bases = freeRes[level+1];
             dr = reducer[level,q];              dr = reducer[level,q];
             dr[q] = -1;              /* dr[q] = -1;  We do not need this in our reducer format. */
               /* dr[q] should be a non-zero constant. */
             newbases = SnewArrayOfFormat(bases);              newbases = SnewArrayOfFormat(bases);
             betti_levelplus = Length(bases);              betti_levelplus = Length(bases);
             /*              /*
                bases[i,j] ---> bases[i,j]+bases[i,q]*dr[j]                 bases[i,j] ---> bases[i,j]+bases[i,q]*dr[j]
             */              */
             for (i=0; i<betti_levelplus; i++) {              for (i=0; i<betti_levelplus; i++) {
               newbases[i] = bases[i] + bases[i,q]*dr;                newbases[i] = dr[q]*bases[i] - bases[i,q]*dr;
             }              }
             Println(["level, q =", level,q]);              Println(["level, q =", level,q]);
             Println("bases="); sm1_pmat(bases);              Println("bases="); sm1_pmat(bases);
Line 1353  def Sminimal_v(g) {
Line 1385  def Sminimal_v(g) {
             minRes[level+1] = newbases;              minRes[level+1] = newbases;
             freeRes = minRes;              freeRes = minRes;
 #ifdef DEBUG  #ifdef DEBUG
 /*  Do it later.              for (qq=q; qq<betti; qq++) {
             for (qq=0; qq<betti; qq++) {  
                 for (i=0; i<betti_levelplus; i++) {                  for (i=0; i<betti_levelplus; i++) {
                   if (!IsZero(newbases[i,qq])) {                    if ((!IsZero(newbases[i,qq])) && (redundantTable[level,qq] >0)) {
                     Println(["[i,qq]=",[i,qq]," is not zero in newbases."]);                      Println(["[i,qq]=",[i,qq]," is not zero in newbases."]);
                     Print("redundantTable ="); sm1_pmat(redundantTable[level]);                      Print("redundantTable ="); sm1_pmat(redundantTable[level]);
                     Error("Stop in Sminimal for debugging.");                      Error("Stop in Sminimal for debugging.");
                   }                    }
                 }                  }
             }              }
 */  
 #endif  #endif
           }            }
         }          }
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.9
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>