[BACK]Return to gb.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / kan96xx / Kan

Diff for /OpenXM/src/kan96xx/Kan/gb.c between version 1.5 and 1.13

version 1.5, 2001/05/04 01:06:23 version 1.13, 2005/07/07 02:59:47
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/kan96xx/Kan/gb.c,v 1.4 2000/03/15 01:31:17 takayama Exp $ */  /* $OpenXM: OpenXM/src/kan96xx/Kan/gb.c,v 1.12 2005/06/30 08:35:07 takayama Exp $ */
 /*  untabify on May 4, 2001 */  /*  untabify on May 4, 2001 */
 #include <stdio.h>  #include <stdio.h>
 #include "datatype.h"  #include "datatype.h"
Line 8 
Line 8 
   
 #define INITGRADE 4  #define INITGRADE 4
 #define INITSIZE 2  #define INITSIZE 2
   #define DMAX 100
   
 int KanGBmessage = 1;  int KanGBmessage = 1;
   
Line 22  extern int UseCriterion2B;
Line 23  extern int UseCriterion2B;
 extern int Spairs;  extern int Spairs;
 extern int Criterion2B, Criterion2F, Criterion2M;  extern int Criterion2B, Criterion2F, Criterion2M;
 extern int AutoReduce;  extern int AutoReduce;
   extern int TraceLift;
   extern struct ring *TraceLift_ringmod;
   static int MaxLength[DMAX];
   static int SpNumber[DMAX];
   
   
 struct gradedPairs *updatePairs(grD,gt,gtGrade,t,grG)  struct gradedPairs *updatePairs(grD,gt,gtGrade,t,grG)
      struct gradedPairs *grD;  /* set of pairs */       struct gradedPairs *grD;  /* set of pairs */
      POLY gt;                 /* new polynomial */       POLY gt;                 /* new polynomial */
Line 158  struct gradedPairs *updatePairs(grD,gt,gtGrade,t,grG)
Line 162  struct gradedPairs *updatePairs(grD,gt,gtGrade,t,grG)
   return(grD);    return(grD);
 }  }
   
 struct gradedPolySet *groebner_gen(f,needBack,needSyz,grP,countDown,forceReduction)  struct gradedPolySet *groebner_gen(f,needBack,needSyz,grP,countDown,forceReduction,reduceOnly,gbCheck)
      struct arrayOfPOLY *f;       struct arrayOfPOLY *f;
      int needBack;       int needBack;
      int needSyz;       int needSyz;
      struct pair **grP;  /* if (needSyz), it is set. */       struct pair **grP;  /* if (needSyz), it is set. */
      int countDown;       int countDown;
      int forceReduction;       int forceReduction;
        int reduceOnly;
        int gbCheck;
 {  {
   int r;    int r;
   struct gradedPolySet *g;    struct gradedPolySet *g;
Line 189  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 195  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
   extern int StopDegree;    extern int StopDegree;
   extern struct ring *CurrentRingp;    extern struct ring *CurrentRingp;
   extern char *F_mpMult;    extern char *F_mpMult;
     extern int Homogenize;
   struct ring *rp;    struct ring *rp;
     int first;
     int statisticsPL, statisticsCount;
   
     if (KanGBmessage) {
           printf("Homogenize=%d, countDown=%d, forceReduction=%d, reduceOnly=%d, gbCheck=%d\n",
                      Homogenize,countDown, forceReduction, reduceOnly, gbCheck);
     }
     if (Statistics) {
       for (i=0; i<DMAX; i++) MaxLength[i] = SpNumber[i] = 0;
     }
   r = f->n;    r = f->n;
   if (r<=0) return((struct gradedPolySet *)NULL);    if (r<=0) return((struct gradedPolySet *)NULL);
   if (UseCriterion1) {    if (UseCriterion1) {
Line 204  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 220  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
   
   Spairs = Criterion1 = Criterion2B = Criterion2F = Criterion2M = 0;    Spairs = Criterion1 = Criterion2B = Criterion2F = Criterion2M = 0;
   
   g = newGradedPolySet(INITGRADE);    g = newGradedPolySet(INITGRADE); g->gb = 1;
     if (reduceOnly) g->gb = 2; /* unknown */
   d = newGradedPairs(INITGRADE*2);    d = newGradedPairs(INITGRADE*2);
   for (i=0; i<g->lim; i++) {    for (i=0; i<g->lim; i++) {
     g->polys[i] = newPolySet(INITSIZE);      g->polys[i] = newPolySet(INITSIZE);
   }    }
   
     first = 1;
   for (i=0; i<r; i++) {    for (i=0; i<r; i++) {
     gt = getArrayOfPOLY(f,i);      gt = getArrayOfPOLY(f,i);
     if (gt ISZERO) { rp = CurrentRingp; } else { rp = gt->m->ringp; }      if (gt ISZERO) { rp = CurrentRingp; } else { rp = gt->m->ringp; }
           if (TraceLift && (!(gt ISZERO)) && first) {
             TraceLift_ringmod = newRingOverFp(rp,getPrime(TraceLift)); first = 0;
             if (KanGBmessage) printf("Prime number for the trace lift is %d.\n",
                                                              TraceLift_ringmod->p);
           }
     grade = -1; whereInG(g,gt,&grade,&indx,Sugar);      grade = -1; whereInG(g,gt,&grade,&indx,Sugar);
     if (KanGBmessage == 2) {      if (KanGBmessage == 2) {
       printf("init=%s, ",POLYToString(head(gt),'*',1));        printf("init=%s, ",POLYToString(head(gt),'*',1));
Line 230  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 253  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
       g = putPolyInG(g,gt,grade,indx,syzp,1,serial);        g = putPolyInG(g,gt,grade,indx,syzp,1,serial);
     }      }
   
     markRedundant0(g,grade,indx);      /* markRedundant0(g,grade,indx); ?*/
           markGeneratorInG(g,grade,indx);  /*?*/
     if (Debug) {      if (Debug) {
       outputGradedPairs(d); outputGradedPolySet(g,needSyz);        outputGradedPairs(d); outputGradedPolySet(g,needSyz);
     }      }
Line 243  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 267  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
   while ((top = getPair(d)) != (struct pair *)NULL) {    while ((top = getPair(d)) != (struct pair *)NULL) {
     ig = top->ig; ii = top->ii; /* [ig,ii] */      ig = top->ig; ii = top->ii; /* [ig,ii] */
     jg = top->jg; ji = top->ji; /* [jg,ji] */      jg = top->jg; ji = top->ji; /* [jg,ji] */
           /*
       if (g->polys[ig]->del[ii] || g->polys[jg]->del[ji]) {
         if (KanGBmessage) printf("p");
         continue;
       }  Don't do this.
           */
     gi = g->polys[ig]->g[ii];      gi = g->polys[ig]->g[ii];
     gj = g->polys[jg]->g[ji];      gj = g->polys[jg]->g[ji];
   
     Spairs++;      Spairs++;
     h = (*sp)(gi,gj);          if (reduceOnly && (!needSyz) && (!needBack)) h = spZero(); /* rd = 0 */
           else h = (*sp)(gi,gj);
     rd = ppAddv(ppMult(h.a,gi),ppMult(h.b,gj));      rd = ppAddv(ppMult(h.a,gi),ppMult(h.b,gj));
   
       if (Statistics) {
         if (top->grade >=0 && top->grade < DMAX) {
           statisticsPL = pLength(rd);
                   SpNumber[top->grade]++;
           if (MaxLength[top->grade] < statisticsPL) {
             MaxLength[top->grade] = statisticsPL;
           }
         }
       }
   
     if (!Sugar || forceReduction) {      if (!Sugar || forceReduction) {
       rd = (*reduction)(rd,g,needBack,&syz);        rd = (*reduction)(rd,g,needBack,&syz);
     }else{      }else{
Line 281  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 323  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
     }      }
   
     if (!(rd ISZERO)) {      if (!(rd ISZERO)) {
         if (gbCheck) {
                   /* Abort the calculation. */
                   g->gb = 0;
                   if (KanGBmessage) {
                     printf("gbCheck failed. \n");
                     printf("Note that the result is NOT groebner basis.\n");
                   }
                   break;
         }
       if (needBack || needSyz) {        if (needBack || needSyz) {
         syzp = newSyz0();          syzp = newSyz0();
         syzp->cf = syzCf; /* no meaning */          syzp->cf = syzCf; /* no meaning */
Line 320  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 371  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
           if (countDown == 0) {            if (countDown == 0) {
             printf("\nThe computation of the Groebner basis is suspended because of countDown==0.\n");              printf("\nThe computation of the Groebner basis is suspended because of countDown==0.\n");
             printf("Note that the result is NOT groebner basis.\n");              printf("Note that the result is NOT groebner basis.\n");
               g->gb = 0;
             break;              break;
           }            }
         }          }
Line 341  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 393  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
         printf("Computation of the Groebner basis is suspended bacause of StopDegree < computing grade.\n");          printf("Computation of the Groebner basis is suspended bacause of StopDegree < computing grade.\n");
         printf("Note that the result is NOT groebner basis.\n");          printf("Note that the result is NOT groebner basis.\n");
       }        }
         g->gb = 0;
       break;        break;
     }      }
   }    }
Line 361  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
Line 414  struct gradedPolySet *groebner_gen(f,needBack,needSyz,
     printf("Criterion1 is applied %d times.\n",Criterion1);      printf("Criterion1 is applied %d times.\n",Criterion1);
     printf("Criterions M,F and B are applied M=%d, F=%d, B=%d times.\n",Criterion2M,Criterion2F,Criterion2B);      printf("Criterions M,F and B are applied M=%d, F=%d, B=%d times.\n",Criterion2M,Criterion2F,Criterion2B);
     Spairs = Criterion1 = Criterion2M = Criterion2F = Criterion2B = 0;      Spairs = Criterion1 = Criterion2M = Criterion2F = Criterion2B = 0;
   
       printf("degree(number of spolys): maximal polynomial size\n");
       statisticsCount = 0;
       for (i=0; i<DMAX; i++) {
         if (SpNumber[i] > 0) {
           printf("%3d(%3d): %5d, ",i,SpNumber[i],MaxLength[i]);
                   if (statisticsCount % 4 == 3) {
                     printf("\n");
                     statisticsCount = 0;
                   }else{ statisticsCount++; }
         }
       }
       printf("\n");
   }    }
   
   if (AutoReduce) {    if (AutoReduce || reduceOnly) {
     toReducedBasis(g,needBack,needSyz);      toReducedBasis(g,needBack,needSyz);
   }    }
   
Line 460  void toReducedBasis(struct gradedPolySet *grP,int need
Line 526  void toReducedBasis(struct gradedPolySet *grP,int need
   if (KanGBmessage) {    if (KanGBmessage) {
     printf("Done(reduced basis)\n");      printf("Done(reduced basis)\n");
   }    }
     grP->reduced = 1;
 }  }
   
   
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.13
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>