| version 1.31, 2015/03/16 00:08:32 |
version 1.32, 2015/08/08 14:19:41 |
|
|
| /* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.30 2003/01/15 10:21:17 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 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) |