version 1.4, 2001/07/03 01:41:25 |
version 1.5, 2001/07/04 07:19:20 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/H.c,v 1.3 2000/08/22 05:04:04 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/H.c,v 1.4 2001/07/03 01:41:25 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "inline.h" |
#include "inline.h" |
|
|
error("henmain2 : cannot happen(extgcd)"); |
error("henmain2 : cannot happen(extgcd)"); |
#endif |
#endif |
|
|
|
#if 1 |
fprintf(stderr,"bound=%d\n",bound); |
fprintf(stderr,"bound=%d\n",bound); |
|
#endif |
for ( k = 1; k < bound; k++ ) { |
for ( k = 1; k < bound; k++ ) { |
|
#if 1 |
fprintf(stderr,"."); |
fprintf(stderr,"."); |
|
#endif |
/* at this point, f = gk*hk mod y^k */ |
/* at this point, f = gk*hk mod y^k */ |
|
|
#if 0 |
#if 0 |
|
|
for ( i = 0; i <= DEG(wh1); i++ ) |
for ( i = 0; i <= DEG(wh1); i++ ) |
COEF(hk)[i][k] = COEF(wh1)[i]; |
COEF(hk)[i][k] = COEF(wh1)[i]; |
} |
} |
|
#if 1 |
fprintf(stderr,"\n"); |
fprintf(stderr,"\n"); |
|
#endif |
*gp = gk; |
*gp = gk; |
clearlum(n,bound,f); |
clearlum(n,bound,f); |
DEG(f) = dh; |
DEG(f) = dh; |
Line 849 ML bqlist,cqlist,*listp; |
|
Line 855 ML bqlist,cqlist,*listp; |
|
for ( j = DEG(b[i]), pp = COEF(l[i]), px = COEF(b[i]); j >= 0; j-- ) |
for ( j = DEG(b[i]), pp = COEF(l[i]), px = COEF(b[i]); j >= 0; j-- ) |
pp[j][0] = px[j]; |
pp[j][0] = px[j]; |
} |
} |
|
#if 0 |
fprintf(stderr,"bound=%d\n",bound); |
fprintf(stderr,"bound=%d\n",bound); |
|
#endif |
for ( i = 1; i < bound; i++ ) { |
for ( i = 1; i < bound; i++ ) { |
|
#if 0 |
fprintf(stderr,"."); |
fprintf(stderr,"."); |
|
#endif |
mullum(mod,i+1,l[0],l[1],wb0); |
mullum(mod,i+1,l[0],l[1],wb0); |
for ( j = 2; j < np; j++ ) { |
for ( j = 2; j < np; j++ ) { |
mullum(mod,i+1,l[j],wb0,wb1); |
mullum(mod,i+1,l[j],wb0,wb1); |
Line 880 ML bqlist,cqlist,*listp; |
|
Line 890 ML bqlist,cqlist,*listp; |
|
for ( j = n, px = COEF(wq0); j >= 0; j-- ) |
for ( j = n, px = COEF(wq0); j >= 0; j-- ) |
px[j] = 0; |
px[j] = 0; |
} |
} |
|
#if 0 |
fprintf(stderr,"\n"); |
fprintf(stderr,"\n"); |
|
#endif |
} |
} |
|
|
static double M; |
static double M; |