version 1.5, 2001/06/25 04:11:42 |
version 1.6, 2001/06/25 06:05:16 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.4 2001/06/25 01:35:21 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.5 2001/06/25 04:11:42 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
|
|
|
|
|
|
q = W_UMALLOC(dx); |
q = W_UMALLOC(dx); |
rlist = MLALLOC(fn); rlist->n = fn; rlist->bound = bound; |
rlist = MLALLOC(fn); rlist->n = fn; rlist->bound = bound; |
for ( i = 0; i < fn-1; i++ ) { |
for ( i = fn-1; i >= 1; i-- ) { |
fprintf(stderr,"%d\n",i); |
fprintf(stderr,"%d\n",i); |
/* fl = gm[i]*hm mod y */ |
/* fl = gm[i]*hm mod y */ |
divsfum(fm,gm[i],hm); |
divsfum(fm,gm[i],hm); |
|
|
sfhenmain2(fl,gm[i],hm,bound,(BM *)&rlist->c[i]); |
sfhenmain2(fl,gm[i],hm,bound,(BM *)&rlist->c[i]); |
cpyum(hm,fm); |
cpyum(hm,fm); |
} |
} |
/* finally, fl must be the lift of gm[fn-1] */ |
/* finally, fl must be the lift of gm[0] */ |
rlist->c[i] = fl; |
rlist->c[i] = fl; |
|
|
/* y -> y-a */ |
/* y -> y-a */ |