version 1.30, 2003/01/15 10:21:17 |
version 1.32, 2015/08/08 14:19:41 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.29 2002/11/22 08:44:57 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.31 2015/03/16 00:08:32 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
#include "inline.h" |
#include "inline.h" |
Line 835 void sfhenmain2(BM f,UM g0,UM h0,int dy,BM *gp) |
|
Line 835 void sfhenmain2(BM f,UM g0,UM h0,int dy,BM *gp) |
|
DEG(f) = dy; |
DEG(f) = dy; |
for ( i = 0; i <= dy; i++ ) |
for ( i = 0; i <= dy; i++ ) |
cpyum(COEF(hk)[i],COEF(f)[i]); |
cpyum(COEF(hk)[i],COEF(f)[i]); |
|
#if defined(__MINGW32__) || defined(__MINGW64__) |
|
fflush(stderr); |
|
#endif |
} |
} |
|
|
/* a0*g+b0*h = 1 mod y -> a*g+b*h = 1 mod y^(dy+1) */ |
/* a0*g+b0*h = 1 mod y -> a*g+b*h = 1 mod y^(dy+1) */ |
Line 925 void sfexgcd_by_hensel(BM g,BM h,int dy,BM *ap,BM *bp) |
|
Line 928 void sfexgcd_by_hensel(BM g,BM h,int dy,BM *ap,BM *bp) |
|
DEG(b) = dy; |
DEG(b) = dy; |
*ap = a; |
*ap = a; |
*bp = b; |
*bp = b; |
|
#if defined(__MINGW32__) || defined(__MINGW64__) |
|
fflush(stderr); |
|
#endif |
} |
} |
|
|
/* fl->c[i] = coef_y(f,i) */ |
/* fl->c[i] = coef_y(f,i) */ |
Line 1005 void sfsqfr(P f,DCP *dcp) |
|
Line 1011 void sfsqfr(P f,DCP *dcp) |
|
} else if ( !NEXT(vl) ) |
} else if ( !NEXT(vl) ) |
sfusqfr(f,dcp); |
sfusqfr(f,dcp); |
else |
else |
sqfrsf(f,dcp); |
sqfrsf(vl,f,dcp); |
} |
} |
|
|
void sfusqfr(P f,DCP *dcp) |
void sfusqfr(P f,DCP *dcp) |
Line 1291 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 1297 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
fprintf(stderr,"total %d, omitted by degtest %d\n",z,z-dtok); |
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; |
|
#if defined(__MINGW32__) || defined(__MINGW64__) |
|
fflush(stderr); |
|
#endif |
} |
} |
|
|
void extractcoefbm(BM f,int dx,UM r) |
void extractcoefbm(BM f,int dx,UM r) |