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

Return to ecart.c CVS log

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

Diff for /OpenXM/src/kan96xx/Kan/ecart.c between version 1.20 and 1.21

version 1.20, 2003/09/21 02:19:43

version 1.21, 2004/09/13 11:24:11

Line 1

/* $OpenXM: OpenXM/src/kan96xx/Kan/ecart.c,v 1.19 2003/09/20 09:57:29 takayama Exp $ */

/* $OpenXM: OpenXM/src/kan96xx/Kan/ecart.c,v 1.20 2003/09/21 02:19:43 takayama Exp $ */

#include <stdio.h>

#include "datatype.h"

#include "extern2.h"

Line 26 static int ecartCheckPoly(POLY f); /* check if it doe

static int ecartCheckEnv(); /* check if the environment is OK for ecart div*/

static struct ecartPolyArray *ecartPutPolyInG(POLY g,struct ecartPolyArray *eparray,POLY cf, POLY syz);

static int ecartGetEll(POLY r,POLY g);

static int ecartGetEllPartial(POLY r,POLY g);

static POLY ecartDivideSv(POLY r,int *d);

/* No automatic homogenization and s is used as a standart var. */

static POLY reduction_ecart0(POLY r,struct gradedPolySet *gset,

Line 117 static int ecartGetEll(POLY f,POLY g) {

Line 118 static int ecartGetEll(POLY f,POLY g) {

else return(-p);

}

static int ecartGetEllPartial(POLY f,POLY g) {

int n,i,p;

MONOMIAL tf;

MONOMIAL tg;

int nglob; int *glob; int k;

if (f ISZERO) return(-1);

if (g ISZERO) return(-1);

checkRingIsR(f,g);

if (!isSameComponent_x(f,g)) return(-1);

tf = f->m; tg = g->m; n = tf->ringp->n;

for (i=1; i<n; i++) {

if (tf->e[i].x < tg->e[i].x) return(-1);

if (tf->e[i].D < tg->e[i].D) return(-1);

}

if (tf->e[0].D < tg->e[0].D) return(-1); /* h */

/* Find u and ell s.t. s^\ell f = u g.

When u contains variables in glob (index x var), then returns -1.

x' = glob, x'' = complement of glob. Division in Q(x'')[x'].

nglob = tf->ringp->partialEcart;

glob = tf->ringp->partialEcartGlobalVarX;

for (i=0; i<nglob; i++) {

k = glob[i];

if (tf->e[k].x > tg->e[k].x) return(-1);

}

p = tf->e[0].x - tg->e[0].x; /* H, s */

if (p >=0 ) return 0;

else return(-p);

}

#define EP_SIZE 10

static struct ecartPolyArray *ecartPutPolyInG(POLY g,struct ecartPolyArray *eparray,POLY cf,POLY syz)

{

Line 203 static struct ecartReducer ecartFindReducer(POLY r,str

Line 239 static struct ecartReducer ecartFindReducer(POLY r,str

if (epa != NULL) {

/* Try to find in the second group. */

for (i=0; i< epa->size; i++) {

ell = ecartGetEll(r,(epa->pa)[i]);

ell = ecartGetEllPartial(r,(epa->pa)[i]);

if ((ell>=0) && (ell < ell2)) {

ell2 = ell;

minGgi = i;

Line 701 static struct ecartReducer ecartFindReducer_mod(POLY r

Line 737 static struct ecartReducer ecartFindReducer_mod(POLY r