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

Return to poly2.c CVS log

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

Diff for /OpenXM/src/kan96xx/Kan/poly2.c between version 1.3 and 1.7

-version 1.3, 2001/05/04 01:06:25
+version 1.7, 2005/07/03 11:08:54
 Line 1
 Line 1
 Line 1
- /* $OpenXM: OpenXM/src/kan96xx/Kan/poly2.c,v 1.2 2000/01/16 07:55:40 takayama Exp $ */
+ /* $OpenXM: OpenXM/src/kan96xx/Kan/poly2.c,v 1.6 2005/06/16 05:07:23 takayama Exp $ */
  #include <stdio.h>
+ #include <stdlib.h>
  #include "datatype.h"
  #include "stackm.h"
  #include "extern.h"
-Line 694  POLY modulo0(f,ringp)
+Line 695  POLY modulo0(f,ringp)
 Line 694  POLY modulo0(f,ringp)
 Line 695  POLY modulo0(f,ringp)
  struct object test(ob)  /* test3 */
           struct object ob;
  {
-   struct object rob;
+   struct object rob = OINIT;
    int k;
    static POLY f0;
    static POLY f1;
-Line 947  POLY mapZmonom(f,ringp)
+Line 948  POLY mapZmonom(f,ringp)
 Line 947  POLY mapZmonom(f,ringp)
 Line 948  POLY mapZmonom(f,ringp)
    node->coeffp->val.f = ff;
    return(node);
  }
+ POLY reduceContentOfPoly(POLY f,struct coeff **contp) {
+   struct coeff *cont;
+   struct coeff *cOne = NULL;
+   extern struct ring *SmallRingp;
+   if (cOne == NULL) cOne = intToCoeff(1,SmallRingp);
+   *contp = cOne;
+   if (f == POLYNULL) return f;
+   if (f->m->ringp->p != 0) return f;
+   if (f->coeffp->tag != MP_INTEGER) return f;
+   cont = gcdOfCoeff(f);
+   *contp = cont;
+   if (coeffGreater(cont,cOne)) {
+         f = quotientByNumber(f,cont).first;
+   }
+   return f;
+ }
+ int coeffSizeMin(POLY f) {
+   int size;
+   int t;
+   if (f == POLYNULL) return 0;
+   if (f->m->ringp->p != 0) return 0;
+   if (f->coeffp->tag != MP_INTEGER) return 0;
+   size = mpz_size(f->coeffp->val.bigp);
+   while (f != POLYNULL) {
+         t = mpz_size(f->coeffp->val.bigp);
+         if (t < size)  size = t;
+         if (size == 1) return size;
+         f = f->next;
+   }
+ }
+ struct coeff *gcdOfCoeff(POLY f) {
+   extern struct ring *SmallRingp;
+   struct coeff *t;
+   MP_INT *tmp;
+   MP_INT *tmp2;
+   static MP_INT *cOne = NULL;
+   if (cOne == NULL) {
+         cOne = newMP_INT();
+         mpz_set_si(cOne,(long) 1);
+   }
+   if (f == POLYNULL) return intToCoeff(0,SmallRingp);
+   if (f->m->ringp->p != 0) return intToCoeff(0,SmallRingp);
+   if (f->coeffp->tag != MP_INTEGER) return intToCoeff(0,SmallRingp);
+   tmp = f->coeffp->val.bigp;
+   tmp2 = newMP_INT();
+   while (f != POLYNULL) {
+     mpz_gcd(tmp2,tmp,f->coeffp->val.bigp); /* tmp = tmp2 OK? */
+         tmp = tmp2;
+         if (mpz_cmp(tmp,cOne)==0) return intToCoeff(1,SmallRingp);
+         f = f->next;
+   }
+   return mpintToCoeff(tmp,SmallRingp);
+ }
+ int shouldReduceContent(POLY f,int ss) {
+   extern DoCancel;
+   static int prevSize = 1;
+   int size;
+   if (f == POLYNULL) return 0;
+   if (f->m->ringp->p != 0) return 0;
+   if (f->coeffp->tag != MP_INTEGER) return 0;
+   if (DoCancel & 2) return 1;
+   /* Apply the Noro strategy to reduce content. */
+   size = mpz_size(f->coeffp->val.bigp);
+   if (ss > 0) {
+         prevSize = size;
+         return 0;
+   }
+   if (size > 2*prevSize) {
+         return 1;
+   }else{
+         return 0;
+   }
+ }

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