version 1.10, 2001/07/03 01:41:25 |
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.9 2001/06/28 08:57:20 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" |
Line 67 void nullspace_lm(LM **,int,int *); |
|
Line 67 void nullspace_lm(LM **,int,int *); |
|
void nullspace_gf2n(GF2N **,int,int *); |
void nullspace_gf2n(GF2N **,int,int *); |
void nullspace_gfpn(GFPN **,int,int *); |
void nullspace_gfpn(GFPN **,int,int *); |
void nullspace_gfs(GFS **,int,int *); |
void nullspace_gfs(GFS **,int,int *); |
|
void nullspace_gfsn(GFSN **,int,int *); |
void null_to_sol(int **,int *,int,int,UM *); |
void null_to_sol(int **,int *,int,int,UM *); |
|
|
void showgfmat(UM **,int); |
void showgfmat(UM **,int); |
|
|
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); |
ptosfbm(dy,f,fl); |
ptosfbm(dy,f,fl); |
shiftsfbm(fl,FTOIF(CONT(ev))); |
if ( ev ) shiftsfbm(fl,FTOIF(CONT(ev))); |
|
|
/* fm = fl mod y */ |
/* fm = fl mod y */ |
fm = W_UMALLOC(dx); |
fm = W_UMALLOC(dx); |
|
|
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; |
|
|
nullspace_gfpn((GFPN **)w,n,ind); break; |
nullspace_gfpn((GFPN **)w,n,ind); break; |
case FF_GFS: |
case FF_GFS: |
nullspace_gfs((GFS **)w,n,ind); break; |
nullspace_gfs((GFS **)w,n,ind); break; |
|
case FF_GFSN: |
|
nullspace_gfsn((GFSN **)w,n,ind); break; |
default: |
default: |
error("nullspace_ff : current_ff is not set"); |
error("nullspace_ff : current_ff is not set"); |
} |
} |
|
|
chsgngfs(u[j],&h); |
chsgngfs(u[j],&h); |
for ( s = j; s < n; s++ ) { |
for ( s = j; s < n; s++ ) { |
mulgfs(h,t[s],&w); addgfs(w,u[s],&w1); u[s] = w1; |
mulgfs(h,t[s],&w); addgfs(w,u[s],&w1); u[s] = w1; |
|
} |
|
} |
|
} |
|
} |
|
|
|
void nullspace_gfsn(mat,n,ind) |
|
GFSN **mat; |
|
int n; |
|
int *ind; |
|
{ |
|
int i,j,l,s; |
|
GFSN w,w1,h,inv; |
|
GFSN *t,*u; |
|
|
|
bzero(ind,n*sizeof(int)); |
|
ind[0] = 0; |
|
|
|
for ( i = j = 0; j < n; i++, j++ ) { |
|
for ( ; j < n; j++ ) { |
|
for ( l = i; l < n; l++ ) |
|
if ( mat[l][j] ) |
|
break; |
|
if ( l < n ) { |
|
t = mat[i]; mat[i] = mat[l]; mat[l] = t; break; |
|
} else |
|
ind[j] = 1; |
|
} |
|
if ( j == n ) |
|
break; |
|
invgfsn(mat[i][j],&inv); |
|
for ( s = j, t = mat[i]; s < n; s++ ) { |
|
mulgfsn(t[s],inv,&w); t[s] = w; |
|
} |
|
for ( l = 0; l < n; l++ ) { |
|
if ( l == i ) |
|
continue; |
|
u = mat[l]; |
|
chsgngfsn(u[j],&h); |
|
for ( s = j; s < n; s++ ) { |
|
mulgfsn(h,t[s],&w); addgfsn(w,u[s],&w1); u[s] = w1; |
} |
} |
} |
} |
} |
} |