OpenXM_contrib2/asir2000/engine/nd.c - diff

Return to nd.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine

Diff for /OpenXM_contrib2/asir2000/engine/nd.c between version 1.183 and 1.192

version 1.183, 2010/02/22 05:27:53

version 1.192, 2010/12/14 05:29:37

Line 1

/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.182 2010/02/22 02:20:43 noro Exp $ */

/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.191 2010/09/27 05:05:58 noro Exp $ */

#include "nd.h"

Line 49 static int nd_demand;

static int nd_module,nd_ispot,nd_mpos,nd_pot_nelim;

static NODE nd_tracelist;

static NODE nd_alltracelist;

static int nd_gentrace,nd_gensyz,nd_nora;

static int nd_gentrace,nd_gensyz,nd_nora,nd_newelim;

static int *nd_gbblock;

NumberField get_numberfield();

UINT *nd_det_compute_bound(NDV **dm,int n,int j);

Line 1917 again:

Line 1918 again:

goto again;

} else if ( nf ) {

if ( checkonly || gensyz ) return 0;

if ( nd_newelim ) {

if ( nd_module ) {

if ( MPOS(HDL(nf)) > 1 ) return 0;

} else if ( !(HDL(nf)[nd_exporigin] & nd_mask[0]) ) return 0;

}

if ( DP_Print ) { printf("+"); fflush(stdout); }

hc = HCU(nf);

nd_removecont(m,nf);

Line 2300 ND_pairs nd_newpairs( NODE g, int t )

Line 2306 ND_pairs nd_newpairs( NODE g, int t )

{

NODE h;

UINT *dl;

int ts,s;

int ts,s,i,t0,min,max;

ND_pairs r,r0;

dl = DL(nd_psh[t]);

Line 2308 ND_pairs nd_newpairs( NODE g, int t )

Line 2314 ND_pairs nd_newpairs( NODE g, int t )

for ( r0 = 0, h = g; h; h = NEXT(h) ) {

if ( nd_module && (MPOS(DL(nd_psh[(long)BDY(h)])) != MPOS(dl)) )

continue;

if ( nd_gbblock ) {

t0 = (long)BDY(h);

for ( i = 0; nd_gbblock[i] >= 0; i += 2 ) {

min = nd_gbblock[i]; max = nd_gbblock[i+1];

if ( t0 >= min && t0 <= max && t >= min && t <= max )

break;

}

if ( nd_gbblock[i] >= 0 )

continue;

}

NEXTND_pairs(r0,r);

r->i1 = (long)BDY(h);

r->i2 = t;

Line 2786 NODE postprocess_algcoef(VL av,NODE alist,NODE r)

Line 2802 NODE postprocess_algcoef(VL av,NODE alist,NODE r)

return u0;

}

void nd_gr(LIST f,LIST v,int m,int f4,struct order_spec *ord,LIST *rp)

void nd_gr(LIST f,LIST v,int m,int homo,int f4,struct order_spec *ord,LIST *rp)

{

VL tv,fv,vv,vc,av;

NODE fd,fd0,r,r0,t,x,s,xx,alist;

int e,max,nvar,i;

NDV b;

int ishomo,nalg,mrank,trank;

int ishomo,nalg,mrank,trank,wmax,len;

NMV a;

Alg alpha,dp;

P p,zp;

Q dmy;

Line 2879 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe

Line 2896 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe

if ( b ) { NEXTNODE(fd0,fd); BDY(fd) = (pointer)b; }

}

if ( fd0 ) NEXT(fd) = 0;

ndv_setup(m,0,fd0,0,0);

if ( !ishomo && homo ) {

for ( t = fd0, wmax = max; t; t = NEXT(t) ) {

b = (NDV)BDY(t); len = LEN(b);

for ( a = BDY(b), i = 0; i < len; i++, NMV_ADV(a) )

wmax = MAX(TD(DL(a)),wmax);

}

homogenize_order(ord,nvar,&ord1);

nd_init_ord(ord1);

nd_setup_parameters(nvar+1,wmax);

for ( t = fd0; t; t = NEXT(t) )

ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos,ompos);

}

ndv_setup(m,0,fd0,nd_gbblock?1:0,0);

if ( nd_gentrace ) {

MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);

}

x = f4?nd_f4(m,&perm):nd_gb(m,ishomo,0,0,&perm);

x = f4?nd_f4(m,&perm):nd_gb(m,ishomo || homo,0,0,&perm);

if ( !x ) {

*rp = 0; return;

}

if ( !ishomo && homo ) {

/* dehomogenization */

for ( t = x; t; t = NEXT(t) ) ndv_dehomogenize((NDV)BDY(t),ord);

nd_init_ord(ord);

nd_setup_parameters(nvar,0);

}

nd_demand = 0;

x = ndv_reducebase(x,perm);

if ( nd_gentrace ) { tl1 = nd_alltracelist; nd_alltracelist = 0; }

Line 2926 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe

Line 2966 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe

}

MKLIST(l1,tl1); MKLIST(l2,tl2); MKLIST(l3,t); MKLIST(l4,tl3);

MKLIST(l5,tl4);

tr = mknode(7,*rp,0,l1,l2,l3,l4,l5); MKLIST(*rp,tr);

tr = mknode(7,*rp,(!ishomo&&homo)?ONE:0,l1,l2,l3,l4,l5); MKLIST(*rp,tr);

}

#if 0

fprintf(asir_out,"ndv_alloc=%d\n",ndv_alloc);

Line 3140 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int

Line 3180 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int

tl1 = tl2 = tl3 = tl4 = 0;

if ( Demand )

nd_demand = 1;

ret = ndv_setup(m,1,fd0,0,0);

ret = ndv_setup(m,1,fd0,nd_gbblock?1:0,0);

if ( nd_gentrace ) {

MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);

}

Line 4896 EPOS nd_create_epos(struct order_spec *ord)

Line 4936 EPOS nd_create_epos(struct order_spec *ord)

/* external interface */

void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec *ord,P *rp)

void nd_nf_p(Obj f,LIST g,LIST v,int m,struct order_spec *ord,Obj *rp)

{

NODE t,in0,in;

ND nd,nf;

ND ndf,nf;

NDV ndv;

NDV ndvf;

VL vv,tv;

int stat,nvar,max,e;

int stat,nvar,max,mrank;

union oNDC dn;

Q cont;

P pp;

LIST ppl;

if ( !f ) {

*rp = 0;

Line 4914 void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec

Line 4955 void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec

pltovl(v,&vv);

for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );

/* get the degree bound */

/* max=65536 implies nd_bpe=32 */

for ( t = BDY(g), max = 1; t; t = NEXT(t) )

max = 65536;

for ( tv = vv; tv; tv = NEXT(tv) ) {

e = getdeg(tv->v,(P)BDY(t));

max = MAX(e,max);

}

for ( tv = vv; tv; tv = NEXT(tv) ) {

e = getdeg(tv->v,f);

max = MAX(e,max);

}

nd_module = 0;

/* nd_module will be set if ord is a module ordering */

nd_init_ord(ord);

nd_setup_parameters(nvar,max);

if ( nd_module && OID(f) != O_LIST )

error("nd_nf_p : the first argument must be a list");

if ( nd_module ) mrank = length(BDY((LIST)f));

/* conversion to ndv */

for ( in0 = 0, t = BDY(g); t; t = NEXT(t) ) {

NEXTNODE(in0,in);

ptozp((P)BDY(t),1,&cont,&pp);

if ( nd_module ) {

BDY(in) = (pointer)ptondv(CO,vv,pp);

if ( !BDY(t) || OID(BDY(t)) != O_LIST

|| length(BDY((LIST)BDY(t))) != mrank )

error("nd_nf_p : inconsistent basis element");

if ( !m ) pltozpl((LIST)BDY(t),&cont,&ppl);

else ppl = (LIST)BDY(t);

BDY(in) = (pointer)pltondv(CO,vv,ppl);

} else {

if ( !m ) ptozp((P)BDY(t),1,&cont,&pp);

else pp = (P)BDY(t);

BDY(in) = (pointer)ptondv(CO,vv,pp);

}

if ( m ) ndv_mod(m,(NDV)BDY(in));

}

NEXTNODE(in0,in);

if ( in0 ) NEXT(in) = 0;

BDY(in) = (pointer)ptondv(CO,vv,f);

if ( m ) ndv_mod(m,(NDV)BDY(in));

NEXT(in) = 0;

if ( nd_module ) ndvf = pltondv(CO,vv,(LIST)f);

else ndvf = ptondv(CO,vv,(P)f);

if ( m ) ndv_mod(m,ndvf);

ndf = (pointer)ndvtond(m,ndvf);

/* dont sort, dont removecont */

ndv_setup(m,0,in0,1,1);

nd_psn--;

nd_scale=2;

while ( 1 ) {

stat = nd_nf(m,0,ndf,nd_ps,1,0,&nf);

nd = (pointer)ndvtond(m,nd_ps[nd_psn]);

if ( !stat )

stat = nd_nf(m,0,nd,nd_ps,1,0,&nf);

error("nd_nf_p : exponent too large");

if ( !stat ) {

if ( nd_module ) *rp = (Obj)ndvtopl(m,CO,vv,ndtondv(m,nf),mrank);

nd_psn++;

else *rp = (Obj)ndvtop(m,CO,vv,ndtondv(m,nf));

nd_reconstruct(0,0);

nd_psn--;

} else

break;

}

*rp = ndvtop(m,CO,vv,ndtondv(m,nf));

}

int nd_to_vect(int mod,UINT *s0,int n,ND d,UINT *r)

Line 6907 void conv_ilist(int demand,int trace,NODE g,int **indp

Line 6949 void conv_ilist(int demand,int trace,NODE g,int **indp

void parse_nd_option(NODE opt)

{

NODE t,p;

NODE t,p,u;

int i,s;

char *key;

Obj value;

nd_gentrace = 0; nd_gensyz = 0; nd_nora = 0;

nd_gentrace = 0; nd_gensyz = 0; nd_nora = 0; nd_gbblock = 0;

nd_newelim = 0;

for ( t = opt; t; t = NEXT(t) ) {

p = BDY((LIST)BDY(t));

key = BDY((STRING)BDY(p));

Line 6922 void parse_nd_option(NODE opt)

Line 6966 void parse_nd_option(NODE opt)

nd_gensyz = value?1:0;

else if ( !strcmp(key,"nora") )

nd_nora = value?1:0;

else if ( !strcmp(key,"gbblock") ) {

if ( !value || OID(value) != O_LIST )

error("nd_* : invalid value for gbblock option");

u = BDY((LIST)value);

nd_gbblock = MALLOC((2*length(u)+1)*sizeof(int));

for ( i = 0; u; u = NEXT(u) ) {

p = BDY((LIST)BDY(u));

s = nd_gbblock[i++] = QTOS((Q)BDY(p));

nd_gbblock[i++] = s+QTOS((Q)BDY(NEXT(p)))-1;

}

nd_gbblock[i] = -1;

} else if ( !strcmp(key,"newelim") )

nd_newelim = value?1:0;

}

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