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

Diff for /OpenXM/src/kan96xx/Kan/resol.c between version 1.2 and 1.5

version 1.2, 2000/01/16 07:55:41 version 1.5, 2000/07/26 02:21:30
Line 1 
Line 1 
 /* $OpenXM$ */  /* $OpenXM: OpenXM/src/kan96xx/Kan/resol.c,v 1.4 2000/06/26 11:15:04 takayama Exp $ */
 /* resol.c */  /* resol.c */
 #include <stdio.h>  #include <stdio.h>
 #include "datatype.h"  #include "datatype.h"
Line 12  static struct arrayOfMonomialSyz schreyerSkelton0(stru
Line 12  static struct arrayOfMonomialSyz schreyerSkelton0(stru
 static struct arrayOfMonomialSyz putMonomialSyz(struct arrayOfMonomialSyz a,  static struct arrayOfMonomialSyz putMonomialSyz(struct arrayOfMonomialSyz a,
                                                 struct monomialSyz *s);                                                  struct monomialSyz *s);
   
   static int RemoveRedundantInSchreyerSkelton = 1;
   
 struct monomialSyz *newMonomialSyz(void)  struct monomialSyz *newMonomialSyz(void)
 {  {
   struct monomialSyz *s;    struct monomialSyz *s;
Line 44  struct arrayOfMonomialSyz enlargeArrayOfMonomialSyz(st
Line 46  struct arrayOfMonomialSyz enlargeArrayOfMonomialSyz(st
 static struct arrayOfMonomialSyz schreyerSkelton0(struct arrayOfPOLY g,int i)  static struct arrayOfMonomialSyz schreyerSkelton0(struct arrayOfPOLY g,int i)
 /* return value will be changed by the next call of this function. */  /* return value will be changed by the next call of this function. */
 {  {
   int m,j,k;    int m,j,k,flag;
   static int s_ij_size = 0;    static int s_ij_size = 0;
   static struct monomialSyz **s_ij = NULL;    static struct monomialSyz **s_ij = NULL;
   struct monomialSyz *s;    struct monomialSyz *s;
   struct spValue sv;    struct spValue sv;
   struct arrayOfMonomialSyz ans;    struct arrayOfMonomialSyz ans;
     extern int RemoveRedundantInSchreyerSkelton;
     /* It was 1. */
   
   m = g.n;    m = g.n;
   if (m > s_ij_size) {    if (m > s_ij_size) {
Line 70  static struct arrayOfMonomialSyz schreyerSkelton0(stru
Line 74  static struct arrayOfMonomialSyz schreyerSkelton0(stru
     }      }
   }    }
   shellForMonomialSyz(s_ij,m-i-1);    shellForMonomialSyz(s_ij,m-i-1);
   for (j=0; j<m-i-1;j++) {    if (RemoveRedundantInSchreyerSkelton) {
     s = s_ij[j];          do {
     if (s->deleted != 1) {            flag = 0;
       for (k=j+1; k<m-i-1; k++) {        for (j=0; j<m-i-1;j++) {
         if (s_ij[k]->deleted != 1) {                  if (s_ij[j]->deleted != 1) {
           if ((*isReducible)(s_ij[k]->a,s->a)) s_ij[k]->deleted = 1;                    for (k=0; k<m-i-1;k++) {
         }                          if ((j != k) && (s_ij[k]->deleted != 1)) {
       }                            if ((*isReducible)(s_ij[k]->a,s_ij[j]->a)) {
     }                                  s_ij[k]->deleted = 1;
                                   flag = 1;
                             }
                           }
                     }
                   }
             }
           }while (flag);
   }    }
   ans.size = m-i-1;    ans.size = m-i-1;
   ans.limit = s_ij_size;    ans.limit = s_ij_size;
Line 113  struct arrayOfMonomialSyz schreyerSkelton(struct array
Line 124  struct arrayOfMonomialSyz schreyerSkelton(struct array
     ipart = schreyerSkelton0(g,i);      ipart = schreyerSkelton0(g,i);
     for (k=0; k< ipart.size; k++) {      for (k=0; k< ipart.size; k++) {
       if ((ipart.p)[k]->deleted != 1) {        if ((ipart.p)[k]->deleted != 1) {
         ans = putMonomialSyz(ans,(ipart.p)[k]);          ans = putMonomialSyz(ans,(ipart.p)[k]);
       }        }
     }      }
   }    }
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>