version 1.17, 2003/01/13 06:40:41 |
version 1.19, 2003/03/12 02:06:11 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Fgfs.c,v 1.16 2003/01/06 09:23:27 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Fgfs.c,v 1.18 2003/01/17 00:17:10 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
|
|
Line 561 void mfctrsfmain(VL vl, P f, DCP *dcp) |
|
Line 561 void mfctrsfmain(VL vl, P f, DCP *dcp) |
|
return; |
return; |
} |
} |
for ( n = 0, tvl = vl; tvl; tvl = NEXT(tvl), n++ ); |
for ( n = 0, tvl = vl; tvl; tvl = NEXT(tvl), n++ ); |
va = (V *)ALLOCA(n*sizeof(int)); |
va = (V *)ALLOCA(n*sizeof(V)); |
da = (int *)ALLOCA(n*sizeof(int)); |
da = (int *)ALLOCA(n*sizeof(int)); |
/* find v s.t. diff(f,v) is nonzero and deg(f,v) is minimal */ |
/* find v s.t. diff(f,v) is nonzero and deg(f,v) is minimal */ |
imin = -1; |
imin = -1; |
Line 979 void mfctrsf_hensel(VL vl,VL rvl,P f,P pp0,P u0,P v0,P |
|
Line 979 void mfctrsf_hensel(VL vl,VL rvl,P f,P pp0,P u0,P v0,P |
|
addp(vl,v,wv,&t); v = t; |
addp(vl,v,wv,&t); v = t; |
addp(vl,vh[j],wv,&t); vh[j] = t; |
addp(vl,vh[j],wv,&t); vh[j] = t; |
} |
} |
|
gfsn_poly_to_poly(vl,u,vy,&t); |
|
shift_sf(vl,rvl,t,mev,-1,&s); |
|
if ( divtp(vl,fin,s,&q) ) |
|
cont_pp_mv_sf(vl,onevl,s,&cont,up); |
|
else |
|
*up = 0; |
} |
} |
|
|
void adjust_coef_sf(VL vl,VL rvl,P lcu,P u0,int *mev,P *r) |
void adjust_coef_sf(VL vl,VL rvl,P lcu,P u0,int *mev,P *r) |