version 1.18, 2001/10/09 01:36:10 |
version 1.20, 2001/10/30 10:24:35 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.17 2001/09/03 07:01:06 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.19 2001/10/30 07:25:58 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
|
|
void lnfsf(int n,UM p0,UM p1,struct p_pair *list,UM np0,UM np1); |
void lnfsf(int n,UM p0,UM p1,struct p_pair *list,UM np0,UM np1); |
|
void extractcoefbm(BM f,int dx,UM r); |
|
|
int comp_dum(DUM a,DUM b) |
int comp_dum(DUM a,DUM b) |
{ |
{ |
Line 517 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
Line 518 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
BM fl; |
BM fl; |
VL vl,nvl; |
VL vl,nvl; |
V y; |
V y; |
int dx,dy; |
int dx,dy,bound; |
GFS ev; |
GFS ev; |
P f1,t,c,sf; |
P f1,t,c,sf; |
DCP dc; |
DCP dc; |
Line 554 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
Line 555 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
ptosfum(dc->c,gm[i]); |
ptosfum(dc->c,gm[i]); |
} |
} |
|
|
|
/* set bound */ |
|
/* g | f, lc_y(g) = lc_y(f) => deg_y(g) <= deg_y(f) */ |
|
/* so, bound = dy is sufficient, but we use slightly large value */ |
|
bound = dy+2; |
|
|
/* f(x,y) -> f(x,y+ev) */ |
/* f(x,y) -> f(x,y+ev) */ |
fl = BMALLOC(dx,dy); |
fl = BMALLOC(dx,bound); |
ptosfbm(dy,f,fl); |
ptosfbm(bound,f,fl); |
if ( ev ) shiftsfbm(fl,FTOIF(CONT(ev))); |
if ( ev ) shiftsfbm(fl,FTOIF(CONT(ev))); |
|
|
/* sf = f(x+ev) */ |
/* sf = f(x+ev) */ |
Line 568 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
Line 574 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
hm = W_UMALLOC(dx); |
hm = W_UMALLOC(dx); |
|
|
q = W_UMALLOC(dx); |
q = W_UMALLOC(dx); |
rlist = MLALLOC(fn); rlist->n = fn; rlist->bound = dy; |
rlist = MLALLOC(fn); rlist->n = fn; rlist->bound = bound; |
fprintf(asir_out,"%d candidates\n",fn); |
fprintf(asir_out,"%d candidates\n",fn); |
init_eg(&eg_hensel); |
init_eg(&eg_hensel); |
for ( i = 0; i < fn-1; i++ ) { |
for ( i = 0; i < fn-1; i++ ) { |
Line 578 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
Line 584 void sfhensel(int count,P f,V x,GFS *evp,P *sfp,ML *li |
|
get_eg(&tmp0); |
get_eg(&tmp0); |
/* fl = gm[i]*hm mod y */ |
/* fl = gm[i]*hm mod y */ |
divsfum(fm,gm[i],hm); |
divsfum(fm,gm[i],hm); |
/* fl is replaced by the cofactor of gk mod y^dy */ |
/* fl is replaced by the cofactor of gk mod y^bound */ |
/* rlist->c[i] = gk */ |
/* rlist->c[i] = gk */ |
sfhenmain2(fl,gm[i],hm,dy,(BM *)&rlist->c[i]); |
sfhenmain2(fl,gm[i],hm,bound,(BM *)&rlist->c[i]); |
cpyum(hm,fm); |
cpyum(hm,fm); |
get_eg(&tmp1); add_eg(&eg_hensel_t,&tmp0,&tmp1); |
get_eg(&tmp1); add_eg(&eg_hensel_t,&tmp0,&tmp1); |
add_eg(&eg_hensel,&tmp0,&tmp1); |
add_eg(&eg_hensel,&tmp0,&tmp1); |
Line 959 void sfbfctr(P f,V x,V y,DCP *dcp) |
|
Line 965 void sfbfctr(P f,V x,V y,DCP *dcp) |
|
void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
{ |
{ |
int np,dx,dy; |
int np,dx,dy; |
int i,j,k; |
int i,j,k,bound; |
int *win; |
int *win; |
P g,lcg,factor,cofactor,lcyx; |
P g,lcg,factor,cofactor,lcyx; |
P csum; |
P csum; |
DCP dcf,dcf0,dc; |
DCP dcf,dcf0,dc; |
BM *c; |
BM *c; |
BM lcy; |
BM lcy; |
UM lcg0; |
UM lcg0,lcy0,w; |
|
UM *d1c; |
ML wlist; |
ML wlist; |
struct oVL vl1,vl0; |
struct oVL vl1,vl0; |
VL vl; |
VL vl; |
int z; |
int z,dt,dtok; |
|
|
/* vl = [x,y] */ |
/* vl = [x,y] */ |
vl0.v = x; vl0.next = &vl1; vl1.v = y; vl1.next = 0; vl = &vl0; |
vl0.v = x; vl0.next = &vl1; vl1.v = y; vl1.next = 0; vl = &vl0; |
Line 982 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 989 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
win = W_ALLOC(np+1); |
win = W_ALLOC(np+1); |
wlist = W_MLALLOC(np); |
wlist = W_MLALLOC(np); |
wlist->n = list->n; |
wlist->n = list->n; |
wlist->bound = dy; |
bound = wlist->bound = list->bound; |
c = (BM *)COEF(wlist); |
c = (BM *)COEF(wlist); |
bcopy((char *)COEF(list),(char *)c,(int)(sizeof(BM)*np)); |
bcopy((char *)COEF(list),(char *)c,(int)(sizeof(BM)*np)); |
|
|
Line 1006 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 1013 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
NEWP(lcyx); VR(lcyx) = x; DC(lcyx) = dc; |
NEWP(lcyx); VR(lcyx) = x; DC(lcyx) = dc; |
ptosfbm(dy,lcyx,lcy); |
ptosfbm(dy,lcyx,lcy); |
|
|
|
/* initialize lcy0 by LC(f) */ |
|
lcy0 = W_UMALLOC(bound); |
|
ptosfum(COEF(DC(g)),lcy0); |
|
|
|
/* ((d-1 coefs)*lcy0 */ |
|
d1c = (UM *)W_ALLOC(np*sizeof(UM)); |
|
w = W_UMALLOC(2*bound); |
|
for ( i = 1; i < np; i++ ) { |
|
extractcoefbm(c[i],degbm(c[i])-1,w); |
|
d1c[i] = W_UMALLOC(2*bound); |
|
mulsfum(w,lcy0,d1c[i]); |
|
/* d1c[i] = d1c[i] mod y^(bound+1) */ |
|
if ( DEG(d1c[i]) > bound ) { |
|
for ( j = DEG(d1c[i]); j > bound; j-- ) |
|
COEF(d1c[i])[j] = 0; |
|
degum(d1c[i],bound); |
|
} |
|
} |
|
|
fprintf(stderr,"np = %d\n",np); |
fprintf(stderr,"np = %d\n",np); |
|
dtok = 0; |
for ( g = f, k = 1, dcf = dcf0 = 0, win[0] = 1, --np, z = 0; ; z++ ) { |
for ( g = f, k = 1, dcf = dcf0 = 0, win[0] = 1, --np, z = 0; ; z++ ) { |
if ( !(z % 1000) ) fprintf(stderr,"."); |
if ( !(z % 1000) ) fprintf(stderr,"."); |
if ( sfdtestmain(vl,lcg,lcg0,lcy,csum,wlist,k,win,&factor,&cofactor) ) { |
dt = sfdegtest(dy,bound,d1c,k,win); |
|
if ( dt ) |
|
dtok++; |
|
if ( dt && sfdtestmain(vl,lcg,lcg0,lcy,csum,wlist, |
|
k,win,&factor,&cofactor) ) { |
NEXTDC(dcf0,dcf); DEG(dcf) = ONE; COEF(dcf) = factor; |
NEXTDC(dcf0,dcf); DEG(dcf) = ONE; COEF(dcf) = factor; |
g = cofactor; |
g = cofactor; |
|
|
Line 1022 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 1053 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
/* update csum */ |
/* update csum */ |
sfcsump(vl,lcg,&csum); |
sfcsump(vl,lcg,&csum); |
|
|
|
/* update dy */ |
|
dy = getdeg(y,g); |
|
|
/* update lcy */ |
/* update lcy */ |
clearbm(0,lcy); |
clearbm(0,lcy); |
COEF(dc) = COEF(DC(g)); |
COEF(dc) = COEF(DC(g)); |
Line 1043 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 1077 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
else |
else |
for ( i = 1; i < k; i++ ) |
for ( i = 1; i < k; i++ ) |
win[i] = win[0] + i; |
win[i] = win[0] + i; |
|
|
|
|
|
/* update lcy0 */ |
|
ptosfum(COEF(DC(g)),lcy0); |
|
|
|
/* update d-1 coeffs */ |
|
for ( i = 1; i <= np; i++ ) { |
|
extractcoefbm(c[i],degbm(c[i])-1,w); |
|
mulsfum(w,lcy0,d1c[i]); |
|
/* d1c[i] = d1c[1] mod y^(bound+1) */ |
|
if ( DEG(d1c[i]) > bound ) { |
|
for ( j = DEG(d1c[i]); j > bound; j-- ) |
|
COEF(d1c[i])[j] = 0; |
|
degum(d1c[i],bound); |
|
} |
|
} |
} else if ( !ncombi(1,np,k,win) ) |
} else if ( !ncombi(1,np,k,win) ) |
if ( k == np ) |
if ( k == np ) |
break; |
break; |
Line 1050 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 1100 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
for ( i = 0, ++k; i < k; i++ ) |
for ( i = 0, ++k; i < k; i++ ) |
win[i] = i + 1; |
win[i] = i + 1; |
} |
} |
fprintf(stderr,"\n"); |
fprintf(stderr,"total %d, omitted by degtest %d\n",z,z-dtok); |
NEXTDC(dcf0,dcf); COEF(dcf) = g; |
NEXTDC(dcf0,dcf); COEF(dcf) = g; |
DEG(dcf) = ONE; NEXT(dcf) = 0; *dcp = dcf0; |
DEG(dcf) = ONE; NEXT(dcf) = 0; *dcp = dcf0; |
} |
} |
|
|
|
void extractcoefbm(BM f,int dx,UM r) |
|
{ |
|
int j; |
|
UM fj; |
|
|
|
for ( j = DEG(f); j >= 0; j-- ) { |
|
fj = COEF(f)[j]; |
|
if ( fj && DEG(fj) >= dx ) { |
|
COEF(r)[j] = COEF(fj)[dx]; |
|
} else |
|
COEF(r)[j] = 0; |
|
} |
|
degum(r,DEG(f)); |
|
} |
|
|
|
/* deg_y(prod mod y^(bound+1)) <= dy ? */ |
|
|
|
int sfdegtest(int dy,int bound,UM *d1c,int k,int *in) |
|
{ |
|
int i,j; |
|
UM w,w1,wt; |
|
BM t; |
|
|
|
w = W_UMALLOC(bound); |
|
w1 = W_UMALLOC(bound); |
|
clearum(w,bound); |
|
for ( i = 0; i < k; i++ ) { |
|
addsfum(w,d1c[in[i]],w1); wt = w; w = w1; w1 = wt; |
|
} |
|
return DEG(w) <= dy ? 1 : 0; |
|
} |
|
|
/* lcy = LC(g), lcg = lcy*g, lcg0 = const part of lcg */ |
/* lcy = LC(g), lcg = lcy*g, lcg0 = const part of lcg */ |
|
|
int sfdtestmain(VL vl,P lcg,UM lcg0,BM lcy,P csum,ML list, |
int sfdtestmain(VL vl,P lcg,UM lcg0,BM lcy,P csum,ML list, |
int k,int *in,P *fp,P *cofp) |
int k,int *in,P *fp,P *cofp) |
{ |
{ |