version 1.12, 2001/09/03 07:01:05 |
version 1.14, 2001/10/09 01:36:06 |
|
|
* 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/builtin/gf.c,v 1.11 2001/08/02 03:59:15 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/gf.c,v 1.13 2001/09/03 07:56:19 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
|
|
BM fl; |
BM fl; |
BM *r; |
BM *r; |
VL vl,nvl; |
VL vl,nvl; |
int i,fn,dx,dy; |
int i,fn,dx,dy,d; |
NODE t,top; |
NODE t,top; |
UM fm,hm,q; |
UM fm,hm,q; |
UM *gm; |
UM *gm; |
|
|
x = VR((P)BDY(mfl)); |
x = VR((P)BDY(mfl)); |
y = vl->v == x ? vl->next->v : vl->v; |
y = vl->v == x ? vl->next->v : vl->v; |
|
|
for ( i = 0, t = mfl; i < fn; i++, t = NEXT(t) ) { |
for ( i = 0, t = mfl, d = 0; i < fn; i++, t = NEXT(t) ) { |
gm[i] = (pointer)UMALLOC(getdeg(x,(P)BDY(t))); |
gm[i] = (pointer)UMALLOC(getdeg(x,(P)BDY(t))); |
ptosfum((P)BDY(t),gm[i]); |
ptosfum((P)BDY(t),gm[i]); |
|
d += DEG(gm[i]); |
} |
} |
|
|
/* reorder f if necessary */ |
/* reorder f if necessary */ |
|
|
vl = nvl; f = g; |
vl = nvl; f = g; |
} |
} |
dx = getdeg(x,f); |
dx = getdeg(x,f); |
|
if ( dx != d ) |
|
error("sfuhensel : product of factors has incompatible degree"); |
|
|
dy = getdeg(y,f); |
dy = getdeg(y,f); |
dy = MAX(dy,bound); |
dy = MAX(dy,bound); |
fl = BMALLOC(dx,dy); |
fl = BMALLOC(dx,dy); |
|
|
LIST *rp; |
LIST *rp; |
{ |
{ |
int i,j,n; |
int i,j,n; |
Q mod; |
|
MAT mat,r; |
MAT mat,r; |
VECT u; |
VECT u; |
Q q; |
Q q; |
|
|
int i,j,l,s; |
int i,j,l,s; |
GFSN w,w1,h,inv; |
GFSN w,w1,h,inv; |
GFSN *t,*u; |
GFSN *t,*u; |
GFSN one; |
|
|
|
bzero(ind,n*sizeof(int)); |
bzero(ind,n*sizeof(int)); |
ind[0] = 0; |
ind[0] = 0; |