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

Diff for /OpenXM/src/kan96xx/Kan/sugar.c between version 1.6 and 1.7

version 1.6, 2003/09/12 02:52:50 version 1.7, 2003/09/20 09:57:29
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/kan96xx/Kan/sugar.c,v 1.5 2003/08/21 02:30:23 takayama Exp $ */  /* $OpenXM: OpenXM/src/kan96xx/Kan/sugar.c,v 1.6 2003/09/12 02:52:50 takayama Exp $ */
 #include <stdio.h>  #include <stdio.h>
 #include "datatype.h"  #include "datatype.h"
 #include "extern2.h"  #include "extern2.h"
Line 11 
Line 11 
 extern DebugReductionRed;  extern DebugReductionRed;
 extern DebugContentReduction;  extern DebugContentReduction;
 extern Ecart;  extern Ecart;
   extern int Ecart_sugarGrade;
   static int ecart_grade_sugar(POLY f);
   
 POLY reduction_sugar(POLY f,struct gradedPolySet *gset,int needSyz,  POLY reduction_sugar(POLY f,struct gradedPolySet *gset,int needSyz,
                      struct syz0 *syzp,int sugarGrade)                       struct syz0 *syzp,int sugarGrade)
Line 28  POLY reduction_sugar(POLY f,struct gradedPolySet *gset
Line 30  POLY reduction_sugar(POLY f,struct gradedPolySet *gset
   struct ring *rp;    struct ring *rp;
   extern DoCancel;    extern DoCancel;
   
   if (Ecart) return reduction_ecart(f,gset,needSyz,syzp);    if (Ecart) {
           Ecart_sugarGrade = sugarGrade;
           return reduction_ecart(f,gset,needSyz,syzp);
     }
   if (needSyz) {    if (needSyz) {
     if (f ISZERO) { rp = CurrentRingp; } else { rp = f->m->ringp; }      if (f ISZERO) { rp = CurrentRingp; } else { rp = f->m->ringp; }
     cf = cxx(1,0,0,rp);      cf = cxx(1,0,0,rp);
Line 167  int grade_sugar(f)
Line 172  int grade_sugar(f)
   static int nn,mm,ll,cc,n,m,l,c;    static int nn,mm,ll,cc,n,m,l,c;
   static struct ring *cr = (struct ring *)NULL;    static struct ring *cr = (struct ring *)NULL;
   
     /* if (Ecart) return ecart_grade_sugar(f); */
   if (f ISZERO) return(-1);    if (f ISZERO) return(-1);
   tf = f->m;    tf = f->m;
   if (tf->ringp != cr) {    if (tf->ringp != cr) {
Line 201  int grade_sugar(f)
Line 207  int grade_sugar(f)
       r += tf->e[i].x;        r += tf->e[i].x;
       r += tf->e[i].D;        r += tf->e[i].D;
     }      }
       f = f->next;
       ans = (ans>r?ans:r);
     }
     return(ans);
   }
   
   static int ecart_grade_sugar(POLY f)  /* experimental */
   {
     int r;
     int i,ans;
     MONOMIAL tf;
     static int nn,mm,ll,cc,n,m,l,c;
     static struct ring *cr = (struct ring *)NULL;
     int soffset;
     int sdegree;
   
     if (f ISZERO) return(-1);
     tf = f->m;
     if (tf->ringp != cr) {
       n = tf->ringp->n;
       m = tf->ringp->m;
       l = tf->ringp->l;
       c = tf->ringp->c;
       nn = tf->ringp->nn;
       mm = tf->ringp->mm;
       ll = tf->ringp->ll;
       cc = tf->ringp->cc;
       cr = tf->ringp;
     }
   
     soffset = 0x7fffffff; /* large number */
     while (f != NULL) {
       r = 0;
       tf = f->m;
       for (i=m; i<nn; i++) {
         r -= tf->e[i].x;
         r += tf->e[i].D;
       }
       f = f->next;
       soffset = (r < soffset? r: soffset);
     }
   
     ans = 0;
     while (f != NULL) {
       r = 0;  sdegree = 0;
       tf = f->m;
       for (i=0; i<cc; i++) {
         r += tf->e[i].x;
         r += tf->e[i].D;
       }
       for (i=c; i<ll; i++) {
         r += tf->e[i].x;
         r += tf->e[i].D;
       }
       for (i=l; i<mm; i++) {
         r += tf->e[i].x;
         r += tf->e[i].D;
       }
       for (i=m; i<nn; i++) {
         r += tf->e[i].x;
         r += tf->e[i].D;
         sdegree -= tf->e[i].x;
         sdegree += tf->e[i].D;
   
       }
           r += sdegree-soffset;
     f = f->next;      f = f->next;
     ans = (ans>r?ans:r);      ans = (ans>r?ans:r);
   }    }
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>