OpenXM/src/kan96xx/Kan/poly4.c - diff

Return to poly4.c CVS log

Up to [local] / OpenXM / src / kan96xx / Kan

Diff for /OpenXM/src/kan96xx/Kan/poly4.c between version 1.9 and 1.11

-version 1.9, 2003/08/21 12:28:58
+version 1.11, 2003/08/23 02:28:38
 Line 1
 Line 1
 Line 1
- /* $OpenXM: OpenXM/src/kan96xx/Kan/poly4.c,v 1.8 2003/08/19 08:02:10 takayama Exp $ */
+ /* $OpenXM: OpenXM/src/kan96xx/Kan/poly4.c,v 1.10 2003/08/22 11:47:03 takayama Exp $ */
  #include <stdio.h>
  #include "datatype.h"
  #include "stackm.h"
 Line 8
 Line 8
 Line 8
  static void shell(int v[],int n);
  static int degreeOfPrincipalPart(POLY f);
  static int degreeOfInitW(POLY f,int w[]);
+ static int degreeOfInitWS(POLY f,int w[],int s[]);
  static void shell(v,n)
-Line 437  POLY POLYToInitW(f,w)
+Line 438  POLY POLYToInitW(f,w)
 Line 437  POLY POLYToInitW(f,w)
 Line 438  POLY POLYToInitW(f,w)
       POLY f;
       int w[]; /* weight vector */
  {
-   POLY node;
-   struct listPoly nod;
    POLY h;
    POLY g;
    int maxd;
    int dd;
-   node = &nod; node->next = POLYNULL; h = node;
+   h = POLYNULL;
+   /*printf("1:%s\n",POLYToString(f,'*',1));*/
    if (f ISZERO) return(f);
    maxd = degreeOfInitW(f,w);
-   g = pCopy(f); /* shallow copy */
+   g = f;
    while (!(f ISZERO)) {
      dd = degreeOfInitW(f,w);
      if (dd > maxd) maxd = dd;
 Line 456  POLY POLYToInitW(f,w)
 Line 456  POLY POLYToInitW(f,w)
 Line 456  POLY POLYToInitW(f,w)
    while (!(g ISZERO)) {
      dd = degreeOfInitW(g,w);
      if (dd == maxd) {
-       h->next = g;
+       h = ppAdd(h,newCell(g->coeffp,g->m)); /* it might be slow. */
-       h = h->next;
      }
      g = g->next;
    }
-   h->next = POLYNULL;
+   /*printf("2:%s\n",POLYToString(h,'*',1));*/
-   return(node->next);
+   return(h);
  }
+ static int degreeOfInitWS(f,w,s)
+      POLY f;
+      int w[];
+          int s[];
+ {
+   int n,i,dd;
+   if (f ISZERO) {
+     errorPoly("degreeOfInitWS(0,w) ");
+   }
+   if (s == (int *) NULL) return degreeOfInitW(f,w);
+   n = f->m->ringp->n; dd = 0;
+   for (i=0; i<n-1; i++) {
+     dd += (f->m->e[i].D)*w[n+i];
+     dd += (f->m->e[i].x)*w[i];
+   }
+   dd += s[(f->m->e[n-1].x)];
+   return(dd);
+ }
+ POLY POLYToInitWS(f,w,s)
+      POLY f;
+      int w[]; /* weight vector */
+          int s[]; /* shift vector */
+ {
+   POLY h;
+   POLY g;
+   int maxd;
+   int dd;
+   h = POLYNULL;
+   /*printf("1s:%s\n",POLYToString(f,'*',1));*/
+   if (f ISZERO) return(f);
+   maxd = degreeOfInitWS(f,w,s);
+   g = f;
+   while (!(f ISZERO)) {
+     dd = degreeOfInitWS(f,w,s);
+     if (dd > maxd) maxd = dd;
+     f = f->next;
+   }
+   while (!(g ISZERO)) {
+     dd = degreeOfInitWS(g,w,s);
+     if (dd == maxd) {
+       h = ppAdd(h,newCell(g->coeffp,g->m)); /* it might be slow. */
+     }
+     g = g->next;
+   }
+   /*printf("2s:%s\n",POLYToString(h,'*',1));*/
+   return(h);
+ }
+ int ordWsAll(f,w,s)
+      POLY f;
+      int w[]; /* weight vector */
+          int s[]; /* shift vector */
+ {
+   int maxd;
+   int dd;
+   if (f ISZERO)  errorPoly("ordWsAll(0,w,s) ");
+   maxd = degreeOfInitWS(f,w,s);
+   while (!(f ISZERO)) {
+     dd = degreeOfInitWS(f,w,s);
+     if (dd > maxd) maxd = dd;
+     f = f->next;
+   }
+   return maxd;
+ }
  /*
 .The substitution  "ringp->multiplication = ...." is allowed only in
    KsetUpRing(), so the check in KswitchFunction is not necessary.
-Line 559  POLY goDeHomogenizeS(POLY f) {
+Line 626  POLY goDeHomogenizeS(POLY f) {
 Line 559  POLY goDeHomogenizeS(POLY f) {
 Line 626  POLY goDeHomogenizeS(POLY f) {
            cp = f->coeffp;
            if (cp->tag == POLY_COEFF) {
                  t = goDeHomogenizeS((cp->val).f);
-                 nc = newCell(polyToCoeff(t,f->m->ringp),f->m);
+                 nc = newCell(polyToCoeff(t,f->m->ringp),monomialCopy(f->m));
                  ans = ppAddv(ans,nc);
                  f = f->next;
            }else{

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>