version 1.32, 2015/08/08 14:19:41 |
version 1.33, 2015/08/14 13:51:54 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.31 2015/03/16 00:08:32 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.32 2015/08/08 14:19:41 fujimoto 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 928 void sfexgcd_by_hensel(BM g,BM h,int dy,BM *ap,BM *bp) |
|
Line 925 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 1297 void sfdtest(P f,ML list,V x,V y,DCP *dcp) |
|
Line 1291 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) |